diff options
| author | schnetter <schnetter@c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5> | 2006-10-19 00:28:38 +0000 |
|---|---|---|
| committer | schnetter <schnetter@c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5> | 2006-10-19 00:28:38 +0000 |
| commit | b076bc5274bea9cebcb5c6b4161fd46bde32b4ce (patch) | |
| tree | 5f06d6de7cc191cf06c22748bfa7825f766451c3 | |
| parent | 858f0f9aa64d43a781fa379e0854209acf94c16a (diff) | |
Abort with an error if asked to apply symmetries to the Weyl scalars.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/RotatingSymmetry90/trunk@39 c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5
| -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); } |