diff options
| author | eschnett <eschnett@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2013-01-20 22:08:39 +0000 |
|---|---|---|
| committer | eschnett <eschnett@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2013-01-20 22:08:39 +0000 |
| commit | f097b0d44bae95b86e1c5640792b066d3fbbfde2 (patch) | |
| tree | 67c421a85c2eee40fddfb6e095b45f0da8eb6229 /src/setup_epsdis.c | |
| parent | fc0fb477a7550b6c07d3f60bfae700171441dde5 (diff) | |
Use cctk_ash instead of cctk_lsh where necessary
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@55 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9
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) |