diff options
Diffstat (limited to 'src/interpolate.c')
-rw-r--r-- | src/interpolate.c | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/src/interpolate.c b/src/interpolate.c index 6fb3b16..0d51359 100644 --- a/src/interpolate.c +++ b/src/interpolate.c @@ -94,6 +94,12 @@ Rot180_CheckTensorTypes (CCTK_ARGUMENTS) || CCTK_EQUALS (tensortypealias, "dd_sym")) { /* symmetric tensor */ assert (numvars == 6); + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { + /* tensor */ + assert (numvars == 9); } else if (CCTK_EQUALS (tensortypealias, "ddd_sym")) { /* 3rd rank tensor, symmetric in last 2 indices */ assert (numvars == 18); @@ -343,6 +349,16 @@ Rot180_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, case 5: break; default: assert(0); } + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { + assert (numvars == 9); + parities[0] = parities[1] = parities[2] = +1; + int const d1 = index % 3; + int const d2 = index / 3; + parities[d1] *= -1; + parities[d2] *= -1; } else if (CCTK_EQUALS (tensortypealias, "ddd_sym")) { assert (numvars == 18); parities[0] = parities[1] = parities[2] = +1; |