diff options
author | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-05-19 21:26:21 +0000 |
---|---|---|
committer | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-05-19 21:26:21 +0000 |
commit | 33bdcffa851c4d7210c24e50ff8964cced12a23e (patch) | |
tree | 2325ef29cb9039064ea0129cc5aad536ded8d86d /src/dissipation.c | |
parent | 62a0750e9f80056d73ddc7256397313e8ed627ae (diff) |
Optimised 8-4 derivative operators. Experimental 8-4 dissipation operators.
Tests of these are on going.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@32 f69c4107-0314-4c4f-9ad4-17e986b73f4a
Diffstat (limited to 'src/dissipation.c')
-rw-r--r-- | src/dissipation.c | 20 |
1 files changed, 19 insertions, 1 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index 50d99a8..638a47a 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -16,6 +16,15 @@ void CCTK_FCALL CCTK_FNAME(dissipation_6_5) (const CCTK_REAL *var, const CCTK_REAL *dx, 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, + const CCTK_INT *nk, + const CCTK_INT *bbox, + const CCTK_INT *gsize, + const CCTK_REAL *dx, + const CCTK_REAL *epsdis, + CCTK_REAL *rhs); static void apply (int const varindex, char const * const optstring, void * const arg); @@ -103,7 +112,16 @@ apply (int const varindex, char const * const optstring, void * const arg) assert (rhsptr); if ( CCTK_Equals(norm_type,"Diagonal") ) { - assert(0); + switch(order) { + case 8: { + CCTK_FNAME(dissipation_8_4) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, dx, &epsdis, rhsptr); + break; + } + default: + assert(0); + } } else { switch(order) { case 6: { |