diff options
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); } |