diff options
Diffstat (limited to 'src/dissipation.c')
-rw-r--r-- | src/dissipation.c | 100 |
1 files changed, 88 insertions, 12 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index 548f522..690fce7 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -85,6 +85,18 @@ void CCTK_FCALL CCTK_FNAME(dissipation_2_1_alt) (const CCTK_REAL *var, const CCTK_REAL *dx, const CCTK_REAL *epsdis, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_2_1_delta) (const CCTK_REAL *var, + const CCTK_INT *ni, + const CCTK_INT *nj, + const CCTK_INT *nk, + const CCTK_INT *bbox, + const CCTK_INT *gsize, + const CCTK_INT *offset, + const CCTK_REAL *dx, + const CCTK_REAL *dy, + const CCTK_REAL *dz, + const CCTK_REAL *epsdis, + CCTK_REAL *rhs); void CCTK_FCALL CCTK_FNAME(dissipation_4_2) (const CCTK_REAL *var, const CCTK_INT *ni, const CCTK_INT *nj, @@ -105,6 +117,18 @@ void CCTK_FCALL CCTK_FNAME(dissipation_4_2_alt) (const CCTK_REAL *var, const CCTK_REAL *dx, const CCTK_REAL *epsdis, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_4_2_delta) (const CCTK_REAL *var, + const CCTK_INT *ni, + const CCTK_INT *nj, + const CCTK_INT *nk, + const CCTK_INT *bbox, + const CCTK_INT *gsize, + const CCTK_INT *offset, + const CCTK_REAL *dx, + const CCTK_REAL *dy, + const CCTK_REAL *dz, + const CCTK_REAL *epsdis, + CCTK_REAL *rhs); void CCTK_FCALL CCTK_FNAME(dissipation_6_3) (const CCTK_REAL *var, const CCTK_INT *ni, const CCTK_INT *nj, @@ -125,6 +149,18 @@ void CCTK_FCALL CCTK_FNAME(dissipation_6_3_alt) (const CCTK_REAL *var, const CCTK_REAL *dx, const CCTK_REAL *epsdis, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_6_3_delta) (const CCTK_REAL *var, + const CCTK_INT *ni, + const CCTK_INT *nj, + const CCTK_INT *nk, + const CCTK_INT *bbox, + const CCTK_INT *gsize, + const CCTK_INT *offset, + const CCTK_REAL *dx, + const CCTK_REAL *dy, + const CCTK_REAL *dz, + const CCTK_REAL *epsdis, + CCTK_REAL *rhs); void CCTK_FCALL CCTK_FNAME(dissipation_8_4) (const CCTK_REAL *var, const CCTK_INT *ni, const CCTK_INT *nj, @@ -145,6 +181,18 @@ void CCTK_FCALL CCTK_FNAME(dissipation_8_4_alt) (const CCTK_REAL *var, const CCTK_REAL *dx, const CCTK_REAL *epsdis, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_8_4_delta) (const CCTK_REAL *var, + const CCTK_INT *ni, + const CCTK_INT *nj, + const CCTK_INT *nk, + const CCTK_INT *bbox, + const CCTK_INT *gsize, + const CCTK_INT *offset, + const CCTK_REAL *dx, + const CCTK_REAL *dy, + const CCTK_REAL *dz, + const CCTK_REAL *epsdis, + CCTK_REAL *rhs); void get_shiftout ( const CCTK_POINTER_TO_CONST cctkGH_, CCTK_INT *offset); void CCTK_FCALL CCTK_FNAME(SBP_Poisoning) ( const CCTK_INT *ni, @@ -264,9 +312,16 @@ apply (int const varindex, char const * const optstring, void * const arg) (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], bbox, gsize, offset, dx, &epsdis, rhsptr); } else { - CCTK_FNAME(dissipation_2_1_alt) - (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], - bbox, gsize, offset, dx, &epsdis, rhsptr); + if (!use_variable_deltas) { + CCTK_FNAME(dissipation_2_1_alt) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, dx, &epsdis, rhsptr); + } else { + CCTK_FNAME(dissipation_2_1_delta) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, + sbp_dx, sbp_dy, sbp_dz, &epsdis, rhsptr); + } } break; } @@ -276,9 +331,16 @@ apply (int const varindex, char const * const optstring, void * const arg) (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], bbox, gsize, offset, dx, &epsdis, rhsptr); } else { - CCTK_FNAME(dissipation_4_2_alt) - (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], - bbox, gsize, offset, dx, &epsdis, rhsptr); + if (!use_variable_deltas) { + CCTK_FNAME(dissipation_4_2_alt) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, dx, &epsdis, rhsptr); + } else { + CCTK_FNAME(dissipation_4_2_delta) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, + sbp_dx, sbp_dy, sbp_dz, &epsdis, rhsptr); + } } break; } @@ -288,9 +350,16 @@ apply (int const varindex, char const * const optstring, void * const arg) (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], bbox, gsize, offset, dx, &epsdis, rhsptr); } else { - CCTK_FNAME(dissipation_6_3_alt) - (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], - bbox, gsize, offset, dx, &epsdis, rhsptr); + if (!use_variable_deltas) { + CCTK_FNAME(dissipation_6_3_alt) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, dx, &epsdis, rhsptr); + } else { + CCTK_FNAME(dissipation_6_3_delta) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, + sbp_dx, sbp_dy, sbp_dz, &epsdis, rhsptr); + } } break; } @@ -300,9 +369,16 @@ apply (int const varindex, char const * const optstring, void * const arg) (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], bbox, gsize, offset, dx, &epsdis, rhsptr); } else { - CCTK_FNAME(dissipation_8_4_alt) - (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], - bbox, gsize, offset, dx, &epsdis, rhsptr); + if (!use_variable_deltas) { + CCTK_FNAME(dissipation_8_4_alt) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, dx, &epsdis, rhsptr); + } else { + CCTK_FNAME(dissipation_8_4_delta) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, offset, + sbp_dx, sbp_dy, sbp_dz, &epsdis, rhsptr); + } } break; } |