diff options
Diffstat (limited to 'src/dissipation.c')
| -rw-r--r-- | src/dissipation.c | 19 |
1 files changed, 17 insertions, 2 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index a1557cc..29dff9f 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -46,6 +46,15 @@ void CCTK_FCALL CCTK_FNAME(dissipation_2_1) (const CCTK_REAL *var, const CCTK_REAL *dx, const CCTK_REAL *epsdis, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_2_1_alt) (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_REAL *dx, + 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, @@ -171,9 +180,15 @@ apply (int const varindex, char const * const optstring, void * const arg) if ( CCTK_Equals(norm_type,"Diagonal") ) { switch(order) { case 2: { - CCTK_FNAME(dissipation_2_1) - (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + if ( CCTK_Equals(dissipation_type,"Mattson-Svard-Nordstrom") ) { + CCTK_FNAME(dissipation_2_1) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], bbox, gsize, dx, &epsdis, rhsptr); + } else { + CCTK_FNAME(dissipation_2_1_alt) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, dx, &epsdis, rhsptr); + } break; } case 4: { |