From b076bc5274bea9cebcb5c6b4161fd46bde32b4ce Mon Sep 17 00:00:00 2001
From: schnetter <schnetter@c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5>
Date: Thu, 19 Oct 2006 00:28:38 +0000
Subject: 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
---
 src/rotatingsymmetry90.c | 19 +++++++++++--------
 1 file 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);
       }
-- 
cgit v1.2.3