From 435cf8c8c7690786cf90fe069cdbb4a659491e2e Mon Sep 17 00:00:00 2001 From: schnetter Date: Tue, 23 Mar 2010 23:03:16 +0000 Subject: 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 --- src/apply.c | 5 +++++ src/interpolate.c | 5 +++++ 2 files changed, 10 insertions(+) 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); -- cgit v1.2.3