diff options
| -rw-r--r-- | src/interpolate.c | 16 | ||||
| -rw-r--r-- | src/rotatingsymmetry180.c | 10 |
2 files changed, 26 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; diff --git a/src/rotatingsymmetry180.c b/src/rotatingsymmetry180.c index c12407c..bbd8dd5 100644 --- a/src/rotatingsymmetry180.c +++ b/src/rotatingsymmetry180.c @@ -189,6 +189,16 @@ int BndRot180VI (cGH const * restrict 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); + paritiess[var][0] = paritiess[var][1] = paritiess[var][2] = +1; + int const d1 = index % 3; + int const d2 = index / 3; + paritiess[var][d1] *= -1; + paritiess[var][d2] *= -1; } else if (CCTK_EQUALS (tensortypealias, "ddd_sym")) { assert (numvars == 12); paritiess[var][0] = paritiess[var][1] = paritiess[var][2] = +1; |