diff options
author | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-14 20:26:02 +0000 |
---|---|---|
committer | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-14 20:26:02 +0000 |
commit | 006083638b7fcad16e5e2c6eed71f72f2c1a2c9c (patch) | |
tree | 5b65477b75adf2b24b5d21473ef1a730c06f288a /src/dissipation.c | |
parent | efcd06a3138fbffba457960d8992111983695e82 (diff) |
Add Frank Herrmann's extra dissipation switches.
Two new parameters, Dissipation::extra_dissipation_in_horizons and
Dissipation::extra_dissipation_at_outerbound, increase the disspation
coefficient inside horizons and near the outer boundary.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@10 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9
Diffstat (limited to 'src/dissipation.c')
-rw-r--r-- | src/dissipation.c | 14 |
1 files changed, 9 insertions, 5 deletions
diff --git a/src/dissipation.c b/src/dissipation.c index 6191cc0..e04c8d5 100644 --- a/src/dissipation.c +++ b/src/dissipation.c @@ -2,6 +2,7 @@ #include <assert.h> #include <stdlib.h> +#include <stdio.h> #include "cctk.h" #include "cctk_Arguments.h" @@ -45,6 +46,9 @@ apply (int const varindex, char const * const optstring, void * const arg) int n; int d; int ierr; + int ni,nj,nk; + int i; + CCTK_REAL minrad,maxrad; assert (varindex >= 0); @@ -86,25 +90,25 @@ apply (int const varindex, char const * const optstring, void * const arg) assert (vargroup >= 0); rhsgroup = CCTK_GroupIndexFromVarI (rhsindex); assert (rhsgroup >= 0); - + ierr = CCTK_GroupData (vargroup, &vardata); assert (!ierr); ierr = CCTK_GroupData (rhsgroup, &rhsdata); assert (!ierr); - + assert (vardata.grouptype == CCTK_GF); assert (vardata.vartype == CCTK_VARIABLE_REAL); assert (vardata.dim == cctk_dim); assert (rhsdata.grouptype == CCTK_GF); assert (rhsdata.vartype == CCTK_VARIABLE_REAL); assert (rhsdata.dim == cctk_dim); - + varptr = CCTK_VarDataPtrI (cctkGH, 0, varindex); assert (varptr); rhsptr = CCTK_VarDataPtrI (cctkGH, 0, rhsindex); assert (rhsptr); - + CCTK_FNAME(apply_dissipation) (varptr, rhsptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], - dx, &dt, &order, &epsdis); + dx, &dt, &order, epsdisA); } |