diff options
Diffstat (limited to 'ML_WaveToy')
| -rw-r--r-- | ML_WaveToy/src/WT_Gaussian.cc | 2 | ||||
| -rw-r--r-- | ML_WaveToy/src/WT_Standing.cc | 2 |
2 files changed, 2 insertions, 2 deletions
diff --git a/ML_WaveToy/src/WT_Gaussian.cc b/ML_WaveToy/src/WT_Gaussian.cc index 2c8f46a..b8862a5 100644 --- a/ML_WaveToy/src/WT_Gaussian.cc +++ b/ML_WaveToy/src/WT_Gaussian.cc @@ -98,7 +98,7 @@ static void WT_Gaussian_Body(cGH const * restrict const cctkGH, int const dir, i /* Calculate temporaries and grid functions */ CCTK_REAL_VEC uL = - kmul(kexp(kmul(kmul(kmul(rL,rL),ToReal(-0.5)),ToReal(ScalarINV(ScalarSQR(width))))),ToReal(amplitude)); + kmul(kexp(kdiv(kmul(kmul(rL,rL),ToReal(-0.5)),ToReal(ScalarSQR(width)))),ToReal(amplitude)); CCTK_REAL_VEC rhoL = ToReal(0); diff --git a/ML_WaveToy/src/WT_Standing.cc b/ML_WaveToy/src/WT_Standing.cc index a3d0fa6..bda635f 100644 --- a/ML_WaveToy/src/WT_Standing.cc +++ b/ML_WaveToy/src/WT_Standing.cc @@ -99,7 +99,7 @@ static void WT_Standing_Body(cGH const * restrict const cctkGH, int const dir, i /* Precompute derivatives */ /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC k = ToReal(Pi*ScalarINV(width)); + CCTK_REAL_VEC k = kdiv(ToReal(Pi),ToReal(width)); CCTK_REAL_VEC omega = ksqrt(kmul(kmul(k,k),ToReal(3))); |