diff options
Diffstat (limited to 'src/dissipation.c')
| -rw-r--r-- | src/dissipation.c | 17 |
1 files changed, 16 insertions, 1 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index 362bfbb..fd59e31 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -13,6 +13,8 @@ CCTK_FNAME(apply_dissipation) (CCTK_REAL const * const var, int const * const ni, int const * const nj, int const * const nk, + CCTK_REAL const * const dx, + CCTK_INT const * const order, CCTK_REAL const * const epsdis); static void @@ -38,11 +40,23 @@ apply (int const varindex, char const * const optstring, void * const arg) cGroup vardata, rhsdata; CCTK_REAL const * varptr; CCTK_REAL * rhsptr; + CCTK_REAL dx[3]; int n; + int d; int ierr; assert (varindex >= 0); + for (d=0; d<cctk_dim; ++d) { + if (cctk_nghostzones[d] < order/2) { + CCTK_WARN (0, "This thorn requires at least order/2 ghost zones"); + } + } + + for (d=0; d<3; ++d) { + dx[d] = CCTK_DELTA_SPACE(d); + } + rhsindex = MoLQueryEvolvedRHS (varindex); if (rhsindex < 0) { char * const fullvarname = CCTK_FullName (varindex); @@ -89,5 +103,6 @@ apply (int const varindex, char const * const optstring, void * const arg) assert (rhsptr); CCTK_FNAME(apply_dissipation) - (varptr, rhsptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], &epsdis); + (varptr, rhsptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + dx, &order, &epsdis); } |