diff options
Diffstat (limited to 'ML_ADMQuantities')
| -rw-r--r-- | ML_ADMQuantities/src/ML_ADMQuantities.cc | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/ML_ADMQuantities/src/ML_ADMQuantities.cc b/ML_ADMQuantities/src/ML_ADMQuantities.cc index 49cf80c..4a65bba 100644 --- a/ML_ADMQuantities/src/ML_ADMQuantities.cc +++ b/ML_ADMQuantities/src/ML_ADMQuantities.cc @@ -17,10 +17,10 @@ /* Define macros used in calculations */ #define INITVALUE (42) -#define QAD(x) (SQR(SQR(x))) -#define INV(x) ((1.0) / (x)) +#define INV(x) ((CCTK_REAL)1.0 / (x)) #define SQR(x) ((x) * (x)) -#define CUB(x) ((x) * (x) * (x)) +#define CUB(x) ((x) * SQR(x)) +#define QAD(x) (SQR(SQR(x))) extern "C" void ML_ADMQuantities_SelectBCs(CCTK_ARGUMENTS) { @@ -140,7 +140,7 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d #pragma omp parallel CCTK_LOOP3(ML_ADMQuantities, i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_ash[0],cctk_ash[1],cctk_ash[2]) { ptrdiff_t const index = di*i + dj*j + dk*k; |