Support tensor type ddd_sym, corresponding to T_i(jk).

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/ReflectionSymmetry/trunk@33 082bdb00-0f4f-0410-b49e-b1835e5f2039

2 files changed, 10 insertions, 0 deletions

diff --git a/src/apply.c b/src/apply.c
index 278903a..c070012 100644
--- a/src/apply.c
+++ b/src/apply.c
@@ -340,6 +340,11 @@ BndReflectVI (cGH const * restrict const cctkGH,
       ttype = SYMTENSOR;
       tcomponent = vi - firstvar - 4;
     }
+  } else if (CCTK_EQUALS (tensortypealias, "ddd_sym")) {
+    /* 3rd rank tensor, symmetric in last 2 indices */
+    assert (numvars == 18);
+    ttype = SYMTENSOR;
+    tcomponent = (vi - firstvar) % 6;
   } else if (CCTK_EQUALS (tensortypealias, "weylscalars_real")) {
     /* Weyl scalars, stored as 10 real values.  NOTE: This assumes
        that Psi_0 comes first, which is NOT the default with
diff --git a/src/interpolate.c b/src/interpolate.c
index b719f41..6868f5a 100644
--- a/src/interpolate.c
+++ b/src/interpolate.c
@@ -277,6 +277,11 @@ ReflectionSymmetry_Interpolate (CCTK_POINTER_TO_CONST restrict const cctkGH_,
           ttype = SYMTENSOR;
           tcomponent = vi - firstvar - 4;
         }
+      } else if (CCTK_EQUALS (tensortypealias, "ddd_sym")) {
+        /* 3rd rank tensor, symmetric in last 2 indices */
+        assert (numvars == 18);
+        ttype = SYMTENSOR;
+        tcomponent = (vi - firstvar) % 6;
       } else if (CCTK_EQUALS (tensortypealias, "weylscalars_real")) {
         /* Weyl scalars, stored as 10 real numbers */
         assert (numvars == 10);