diff options
-rw-r--r-- | src/rotatingsymmetry90.c | 19 |
1 files changed, 11 insertions, 8 deletions
diff --git a/src/rotatingsymmetry90.c b/src/rotatingsymmetry90.c index 51150c4..5099169 100644 --- a/src/rotatingsymmetry90.c +++ b/src/rotatingsymmetry90.c @@ -285,13 +285,11 @@ int BndRot90VI (cGH const * restrict const cctkGH, } else if (CCTK_EQUALS (tensortypealias, "4scalar")) { /* 4-scalar */ } else if (CCTK_EQUALS (tensortypealias, "u") - || CCTK_EQUALS (tensortypealias, "d")) - { + || CCTK_EQUALS (tensortypealias, "d")) { /* vector */ assert (numvars == DIM); } else if (CCTK_EQUALS (tensortypealias, "4u") - || CCTK_EQUALS (tensortypealias, "4d")) - { + || CCTK_EQUALS (tensortypealias, "4d")) { /* 4-vector */ assert (numvars == DIM+1); } else if (CCTK_EQUALS (tensortypealias, "uu_sym") @@ -302,6 +300,9 @@ int BndRot90VI (cGH const * restrict const cctkGH, || CCTK_EQUALS (tensortypealias, "4dd_sym")) { /* symmetric 4-tensor */ assert (numvars == (DIM+1)*(DIM+2)/2); + } else if (CCTK_EQUALS (tensortypealias, "weylscalars_real")) { + /* Weyl scalars, stored as 10 real numbers */ + assert (numvars == 10); } else { char * groupname = CCTK_GroupName(gis[var]); assert (groupname); @@ -322,16 +323,14 @@ int BndRot90VI (cGH const * restrict const cctkGH, srcvi = vis[var]; /* do nothing */ } else if (CCTK_EQUALS (tensortypealias, "u") - || CCTK_EQUALS (tensortypealias, "d")) - { + || CCTK_EQUALS (tensortypealias, "d")) { /* vector */ int srcindex; assert (index>=0 && index<DIM); srcindex = convert_index (step, index, alldirs, &parities[var]); srcvi = firstvar + srcindex; } else if (CCTK_EQUALS (tensortypealias, "4u") - || CCTK_EQUALS (tensortypealias, "4d")) - { + || CCTK_EQUALS (tensortypealias, "4d")) { /* 4-vector */ assert (index>=0 && index<DIM+1); if (index==0) { @@ -399,6 +398,10 @@ int BndRot90VI (cGH const * restrict const cctkGH, } srcvi = firstvar + srcindex; } + } else if (CCTK_EQUALS (tensortypealias, "weylscalars_real")) { + /* Weyl scalars, stored as 10 reals */ + assert (index>=0 && index<10); + srcvi = vis[var]; } else { assert (0); } |