diff options
| author | eschnett <eschnett@c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5> | 2011-09-26 01:38:46 +0000 |
|---|---|---|
| committer | eschnett <eschnett@c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5> | 2011-09-26 01:38:46 +0000 |
| commit | cfba10399922ad0bf2cbb270e41e8fe4df3f092b (patch) | |
| tree | 12ea12653a1f4de6944d7d4bdf3af42389573a27 | |
| parent | ed6e2d41b0a069023dbfbee722b2c138b72700d7 (diff) | |
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
| -rw-r--r-- | src/interpolate.c | 10 | ||||
| -rw-r--r-- | 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; |