diff options
Diffstat (limited to 'src/setup.c')
-rw-r--r-- | src/setup.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/src/setup.c b/src/setup.c index 0744cf8..982bbe9 100644 --- a/src/setup.c +++ b/src/setup.c @@ -1,6 +1,7 @@ /* $Header$ */ #include <assert.h> +#include <stdio.h> #include "cctk.h" #include "cctk_Arguments.h" @@ -14,13 +15,36 @@ void SphericalSurfaceInfo_Setup (CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS; CCTK_REAL const pi = 3.1415926535897932384626433832795028841971693993751; + + int group; + cGroup groupinfo; + cGroupDynamicData groupdata; + int n; + int ierr; + + + + if (nsurfaces == 0) return; + + + + group = CCTK_GroupIndex ("SphericalSurface::sf_radius"); + assert (group>=0); + ierr = CCTK_GroupData (group, &groupinfo); + assert (!ierr); + ierr = CCTK_GroupDynamicData (cctkGH, group, &groupdata); + assert (!ierr); for (n=0; n<nsurfaces; ++n) { /* internal consistency checks */ + assert (groupdata.dim == 2); + assert (groupdata.gsh[0] >= ntheta[n]); + assert (groupdata.gsh[1] >= nphi[n]); + assert (ntheta[n] >= 3*nghoststheta[n] && ntheta[n] <= maxntheta); assert (nphi[n] >= 3*nghostsphi[n] && nphi[n] <= maxnphi); |