diff options
Diffstat (limited to 'src/rotatingsymmetry180.c')
-rw-r--r-- | src/rotatingsymmetry180.c | 24 |
1 files changed, 20 insertions, 4 deletions
diff --git a/src/rotatingsymmetry180.c b/src/rotatingsymmetry180.c index 9cbde1d..cf20492 100644 --- a/src/rotatingsymmetry180.c +++ b/src/rotatingsymmetry180.c @@ -164,14 +164,12 @@ int BndRot180VI (cGH const * restrict const cctkGH, } else if (CCTK_EQUALS (tensortypealias, "4scalar")) { paritiess[var][0] = paritiess[var][1] = paritiess[var][2] = +1; } else if (CCTK_EQUALS (tensortypealias, "u") - || CCTK_EQUALS (tensortypealias, "d")) - { + || CCTK_EQUALS (tensortypealias, "d")) { assert (numvars == 3); paritiess[var][0] = paritiess[var][1] = paritiess[var][2] = +1; paritiess[var][index] = -1; } else if (CCTK_EQUALS (tensortypealias, "4u") - || CCTK_EQUALS (tensortypealias, "4d")) - { + || CCTK_EQUALS (tensortypealias, "4d")) { assert (numvars == 4); if (index == 0) { paritiess[var][0] = paritiess[var][1] = paritiess[var][2] = +1; @@ -212,6 +210,24 @@ int BndRot180VI (cGH const * restrict const cctkGH, default: assert(0); } } + } else if (CCTK_EQUALS (tensortypealias, "weylscalars_real")) { + assert (numvars == 10); + { + static int const weylparities[10][3] = + {{+1,+1,+1}, + {-1,-1,-1}, + {+1,+1,-1}, + {-1,-1,+1}, + {+1,+1,+1}, + {-1,-1,-1}, + {+1,+1,-1}, + {-1,-1,+1}, + {+1,+1,+1}, + {-1,-1,-1}}; + for (d=0; d<3; ++d) { + paritiess[var][d] = weylparities[index][d]; + } + } } else { char * groupname = CCTK_GroupName(gis[var]); assert (groupname); |