Added "bitant_rotate" and "quadrant_reflect_rotate" domains for

systems which have a rotational symmetry.
There are some notes describing these domains and the implementation
in the doc directory: rotating_sym.tex
Currently only works for a single processor.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/CartGrid3D/trunk@161 c78560ca-4b45-4335-b268-5f3340f3cb52

1 files changed, 14 insertions, 0 deletions

diff --git a/src/CartGrid3D.c b/src/CartGrid3D.c
index 64792f5..c4674aa 100644
--- a/src/CartGrid3D.c
+++ b/src/CartGrid3D.c
@@ -16,6 +16,7 @@
 #include "cctk.h"
 #include "cctk_Arguments.h"
 #include "cctk_Parameters.h"
+#include "Symmetry.h"
 
 static const char *rcsid = "$Header$";
 
@@ -425,6 +426,19 @@ void CartGrid3D(CCTK_ARGUMENTS)
   } /* if (coarsest refinement level) */
   
   
+  if ((domainsym[0]==GFSYM_ROTATION_Y || domainsym[2]==GFSYM_ROTATION_X)
+      && (lowerz != -upperz))
+    CCTK_WARN(0, "minimum z must equal maximum z for rotation symmetry");
+
+  if ((domainsym[0]==GFSYM_ROTATION_Z || domainsym[4]==GFSYM_ROTATION_X) 
+      && (lowery != -uppery))
+    CCTK_WARN(0, "minimum y must equal maximum y for rotation symmetry");
+
+  if ((domainsym[2]==GFSYM_ROTATION_Z || domainsym[4]==GFSYM_ROTATION_Y)
+       && (lowerx != -upperx))
+    CCTK_WARN(0, "minimum x must equal maximum x for rotation symmetry");
+
+
 #ifdef CCTK_DEBUG
   printf("\n");
   printf("CartGrid3D\n");