diff options
Diffstat (limited to 'src/setup_epsdis.c')
| -rw-r--r-- | src/setup_epsdis.c | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/src/setup_epsdis.c b/src/setup_epsdis.c index 8703c63..7afce0d 100644 --- a/src/setup_epsdis.c +++ b/src/setup_epsdis.c @@ -14,7 +14,7 @@ #include "util_Table.h" #define MAXDIM 3 -#define REFLEVEL ((int)(0.1 + log10((CCTK_REAL)(cctk_levfac[0]))/log10(2.0))) +#define REFLEVEL (ilogb(cctk_levfac[0])) void setup_epsdis (CCTK_ARGUMENTS) @@ -22,6 +22,7 @@ setup_epsdis (CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; + int ai,aj,ak; int ni,nj,nk; int i,j,k,s,l,m; int index,indexP; @@ -43,15 +44,18 @@ setup_epsdis (CCTK_ARGUMENTS) (double)cctk_time); } + ai=cctk_ash[0]; + aj=cctk_ash[1]; + ak=cctk_ash[2]; ni=cctk_lsh[0]; nj=cctk_lsh[1]; nk=cctk_lsh[2]; if (epsdis_for_level[reflvl] > 0.0) - for (i=0; i<ni*nj*nk; ++i) + for (i=0; i<ai*aj*ak; ++i) epsdisA[i] = epsdis_for_level[reflvl]; else - for (i=0;i<ni*nj*nk;i++) + for (i=0;i<ai*aj*ak;i++) epsdisA[i]=epsdis; if (extra_dissipation_at_outerbound) |