diff options
Diffstat (limited to 'src/rotatingsymmetry90.c')
-rw-r--r-- | src/rotatingsymmetry90.c | 64 |
1 files changed, 42 insertions, 22 deletions
diff --git a/src/rotatingsymmetry90.c b/src/rotatingsymmetry90.c index 2590749..42776d8 100644 --- a/src/rotatingsymmetry90.c +++ b/src/rotatingsymmetry90.c @@ -459,16 +459,18 @@ int BndRot90VI (cGH const * restrict const cctkGH, case CCTK_VARIABLE_INT: /* do nothing */ break; - case CCTK_VARIABLE_REAL: + case CCTK_VARIABLE_REAL: { + CCTK_REAL * restrict const varptr = varptrs[var]; for (k=imin[0]; k<imax[2]; ++k) { for (j=imin[1]; j<imax[1]; ++j) { for (i=imin[2]; i<imax[0]; ++i) { const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); - ((CCTK_REAL * restrict) varptrs[var]) [ind] = poison_value; + memset (&varptr[ind], poison_value, sizeof varptr[ind]); } } } break; + } case CCTK_VARIABLE_COMPLEX: /* do nothing */ break; @@ -575,20 +577,35 @@ int BndRot90VI (cGH const * restrict const cctkGH, case CCTK_VARIABLE_INT: /* do nothing */ break; - case CCTK_VARIABLE_REAL: + case CCTK_VARIABLE_REAL: { + int poison_found = 0; + CCTK_REAL const * restrict const varptr = varptrs[var]; + CCTK_REAL poison; + memset (&poison, poison_value, sizeof poison); for (k=imin[0]; k<imax[2]; ++k) { for (j=imin[1]; j<imax[1]; ++j) { for (i=imin[2]; i<imax[0]; ++i) { const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); - CCTK_REAL const val = - ((CCTK_REAL const * restrict) varptrs[var]) [ind]; - if (val == poison_value) { - CCTK_WARN (CCTK_WARN_ABORT, "Poison found in symmetry regions -- there is an error in this thorn"); + if (memcmp (&varptr[ind], &poison, sizeof poison) == 0) { + printf (" ijk=[%d,%d,%d]\n", i, j, k); + poison_found = 1; } } } } + if (poison_found) { + printf ("Poison found:\n"); + printf (" levfac=[%d,%d,%d]\n", cctk_levfac[0], cctk_levfac[1], cctk_levfac[2]); + printf (" origin_space=[%g,%g,%g]\n", cctk_origin_space[0], cctk_origin_space[1], cctk_origin_space[2]); + printf (" delta_space=[%g,%g,%g]\n", cctk_delta_space[0], cctk_delta_space[1], cctk_delta_space[2]); + printf (" lbnd=[%d,%d,%d]\n", cctk_lbnd[0], cctk_lbnd[1], cctk_lbnd[2]); + printf (" lsh=[%d,%d,%d]\n", cctk_lsh[0], cctk_lsh[1], cctk_lsh[2]); + printf (" gsh=[%d,%d,%d]\n", cctk_gsh[0], cctk_gsh[1], cctk_gsh[2]); + printf (" bbox=[%d,%d,%d,%d,%d,%d]\n", cctk_bbox[0], cctk_bbox[1], cctk_bbox[2], cctk_bbox[3], cctk_bbox[4], cctk_bbox[5]); + CCTK_WARN (CCTK_WARN_ABORT, "Poison found in symmetry regions -- there is an error in this thorn"); + } break; + } case CCTK_VARIABLE_COMPLEX: /* do nothing */ break; @@ -619,41 +636,44 @@ int BndRot90VI (cGH const * restrict const cctkGH, } assert (group.dim == DIM); switch (group.vartype) { - case CCTK_VARIABLE_INT: + case CCTK_VARIABLE_INT: { + CCTK_INT * restrict const varptr = varptrs[var]; for (k=imin[2]; k<imax[2]; ++k) { for (j=imin[1]; j<imax[1]; ++j) { for (i=imin[0]; i<imax[0]; ++i) { const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); - ((CCTK_INT * restrict) varptrs[var]) [ind] *= -1; + varptr[ind] *= -1; } } } break; - case CCTK_VARIABLE_REAL: + } + case CCTK_VARIABLE_REAL: { + CCTK_REAL * restrict const varptr = varptrs[var]; for (k=imin[2]; k<imax[2]; ++k) { for (j=imin[1]; j<imax[1]; ++j) { for (i=imin[0]; i<imax[0]; ++i) { const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); - ((CCTK_REAL * restrict) varptrs[var]) [ind] *= -1; + varptr[ind] *= -1; } } } break; - case CCTK_VARIABLE_COMPLEX: - { - CCTK_COMPLEX const czero = CCTK_Cmplx (0.0, 0.0); - for (k=imin[2]; k<imax[2]; ++k) { - for (j=imin[1]; j<imax[1]; ++j) { - for (i=imin[0]; i<imax[0]; ++i) { - const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); - CCTK_COMPLEX * restrict const ptr = - & ((CCTK_COMPLEX * restrict) varptrs[var]) [ind]; - * ptr = CCTK_CmplxSub (czero, * ptr); - } + } + case CCTK_VARIABLE_COMPLEX: { + CCTK_COMPLEX const czero = CCTK_Cmplx (0.0, 0.0); + CCTK_COMPLEX * restrict const varptr = varptrs[var]; + for (k=imin[2]; k<imax[2]; ++k) { + for (j=imin[1]; j<imax[1]; ++j) { + for (i=imin[0]; i<imax[0]; ++i) { + const int ind = CCTK_GFINDEX3D(cctkGH, i, j, k); + CCTK_COMPLEX * restrict const ptr = & varptr[ind]; + * ptr = CCTK_CmplxSub (czero, * ptr); } } } break; + } default: assert (0); } /* switch grouptype */ |