From cfba10399922ad0bf2cbb270e41e8fe4df3f092b Mon Sep 17 00:00:00 2001 From: eschnett Date: Mon, 26 Sep 2011 01:38:46 +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/RotatingSymmetry90/trunk@65 c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5 --- src/interpolate.c | 10 ++++++++-- src/rotatingsymmetry90.c | 10 ++++++++-- 2 files changed, 16 insertions(+), 4 deletions(-) diff --git a/src/interpolate.c b/src/interpolate.c index 50e3b70..3c73174 100644 --- a/src/interpolate.c +++ b/src/interpolate.c @@ -355,7 +355,10 @@ Rot90_CheckTensorTypes (CCTK_ARGUMENTS) { /* 4-vector */ assert (numvars == DIM+1); - } else if (CCTK_EQUALS (tensortypealias, "du")) { + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { /* tensor */ assert (numvars == DIM*DIM); } else if (CCTK_EQUALS (tensortypealias, "uu_sym") @@ -718,7 +721,10 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, basevar = firstvar + off; var = vi - basevar; } - } else if (CCTK_EQUALS (tensortypealias, "du")) { + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { /* tensor */ assert (numvars == DIM*DIM); tensortype = &tensor; diff --git a/src/rotatingsymmetry90.c b/src/rotatingsymmetry90.c index cd59d46..2e4d8ad 100644 --- a/src/rotatingsymmetry90.c +++ b/src/rotatingsymmetry90.c @@ -293,7 +293,10 @@ int BndRot90VI (cGH const * restrict const cctkGH, || CCTK_EQUALS (tensortypealias, "4d")) { /* 4-vector */ assert (numvars == DIM+1); - } else if (CCTK_EQUALS (tensortypealias, "du")) { + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { /* tensor */ assert (numvars == DIM*DIM); } else if (CCTK_EQUALS (tensortypealias, "uu_sym") @@ -351,7 +354,10 @@ int BndRot90VI (cGH const * restrict const cctkGH, = off + convert_index (step, index-off, alldirs, &parities[var]); srcvi = firstvar + srcindex; } - } else if (CCTK_EQUALS (tensortypealias, "du")) { + } else if (CCTK_EQUALS (tensortypealias, "uu") + || CCTK_EQUALS (tensortypealias, "ud") + || CCTK_EQUALS (tensortypealias, "du") + || CCTK_EQUALS (tensortypealias, "dd")) { /* tensor */ assert (numvars == DIM*DIM); int index1, index2; -- cgit v1.2.3