diff options
| -rw-r--r-- | src/CartGrid3D.c | 26 |
1 files changed, 14 insertions, 12 deletions
diff --git a/src/CartGrid3D.c b/src/CartGrid3D.c index 6a50332..af4f31d 100644 --- a/src/CartGrid3D.c +++ b/src/CartGrid3D.c @@ -423,21 +423,23 @@ void CartGrid3D(CCTK_ARGUMENTS) "y=>",loweryi,upperyi, "z=>",lowerzi,upperzi); CCTK_INFO(infoline); + + + + if ((domainsym[0]==GFSYM_ROTATION_Y || domainsym[2]==GFSYM_ROTATION_X) + && (lowerz + upperz > 1e-12)) + 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 > 1e-12)) + 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 > 1e-12)) + CCTK_WARN(0, "minimum x must equal maximum x for rotation symmetry"); } /* if (coarsest refinement level) */ - - - if ((domainsym[0]==GFSYM_ROTATION_Y || domainsym[2]==GFSYM_ROTATION_X) - && (lowerz + upperz > 1e-12)) - 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 > 1e-12)) - 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 > 1e-12)) - CCTK_WARN(0, "minimum x must equal maximum x for rotation symmetry"); #ifdef CCTK_DEBUG |