aboutsummaryrefslogtreecommitdiff
path: root/src/dissipation.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/dissipation.c')
-rw-r--r--src/dissipation.c14
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);
}