From 119825afea3a659098d5a6f08c433d3feb99fcdf Mon Sep 17 00:00:00 2001 From: eschnett Date: Mon, 26 Sep 2011 01:39:09 +0000 Subject: Implement the "dd" tensor type Implement the "dd" tensor type, a full 3x3 tensor without symmetries. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/RotatingSymmetry180/trunk@65 20f44201-0f4f-0410-9130-e5fc2714a787 --- src/interpolate.c | 16 ++++++++++++++++ src/rotatingsymmetry180.c | 10 ++++++++++ 2 files changed, 26 insertions(+) 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; -- cgit v1.2.3