diff options
author | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-06-15 15:52:08 +0000 |
---|---|---|
committer | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-06-15 15:52:08 +0000 |
commit | 2be746a547a1c999bdb068e1023e4d1ae62004e0 (patch) | |
tree | bf914c315bd2a3b436ae7863e4c618576ae16218 /src/dissipation.c | |
parent | e08a6676391c01d9dce81e260c783f24c043c1d6 (diff) |
Added an optimized 6-3 operator and dissipation operators for the remaining
diagonal norm operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@38 f69c4107-0314-4c4f-9ad4-17e986b73f4a
Diffstat (limited to 'src/dissipation.c')
-rw-r--r-- | src/dissipation.c | 45 |
1 files changed, 45 insertions, 0 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index 2b28950..6426431 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -27,6 +27,33 @@ void CCTK_FCALL CCTK_FNAME(dissipation_6_5_opt) (const CCTK_REAL *var, const CCTK_REAL *epsdis, const CCTK_REAL *diss_fraction, CCTK_REAL *rhs); +void CCTK_FCALL CCTK_FNAME(dissipation_2_1) (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, + 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_6_3) (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_8_4) (const CCTK_REAL *var, const CCTK_INT *ni, const CCTK_INT *nj, @@ -124,6 +151,24 @@ 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], + bbox, gsize, dx, &epsdis, rhsptr); + break; + } + case 4: { + CCTK_FNAME(dissipation_4_2) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, dx, &epsdis, rhsptr); + break; + } + case 6: { + CCTK_FNAME(dissipation_6_3) + (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], + bbox, gsize, dx, &epsdis, rhsptr); + break; + } case 8: { CCTK_FNAME(dissipation_8_4) (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], |