diff options
Diffstat (limited to 'ML_ADM/src/ML_ADM_Minkowski.cc')
-rw-r--r-- | ML_ADM/src/ML_ADM_Minkowski.cc | 38 |
1 files changed, 19 insertions, 19 deletions
diff --git a/ML_ADM/src/ML_ADM_Minkowski.cc b/ML_ADM/src/ML_ADM_Minkowski.cc index 7d2462a..38943f5 100644 --- a/ML_ADM/src/ML_ADM_Minkowski.cc +++ b/ML_ADM/src/ML_ADM_Minkowski.cc @@ -90,9 +90,9 @@ static void ML_ADM_Minkowski_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC const p1o840dx = kdiv(ToReal(0.00119047619047619047619047619048),dx); CCTK_REAL_VEC const p1o840dy = kdiv(ToReal(0.00119047619047619047619047619048),dy); CCTK_REAL_VEC const p1o840dz = kdiv(ToReal(0.00119047619047619047619047619048),dz); - CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1.),kmul(dx,dx)); - CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1.),kmul(dy,dy)); - CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1.),kmul(dz,dz)); + CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx)); + CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy)); + CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz)); CCTK_REAL_VEC const pm1o12dx2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx)); CCTK_REAL_VEC const pm1o12dy2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy)); CCTK_REAL_VEC const pm1o12dz2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz)); @@ -184,37 +184,37 @@ static void ML_ADM_Minkowski_Body(cGH const * restrict const cctkGH, int const d } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC g11L = ToReal(1.); + CCTK_REAL_VEC g11L = ToReal(1); - CCTK_REAL_VEC g12L = ToReal(0.); + CCTK_REAL_VEC g12L = ToReal(0); - CCTK_REAL_VEC g13L = ToReal(0.); + CCTK_REAL_VEC g13L = ToReal(0); - CCTK_REAL_VEC g22L = ToReal(1.); + CCTK_REAL_VEC g22L = ToReal(1); - CCTK_REAL_VEC g23L = ToReal(0.); + CCTK_REAL_VEC g23L = ToReal(0); - CCTK_REAL_VEC g33L = ToReal(1.); + CCTK_REAL_VEC g33L = ToReal(1); - CCTK_REAL_VEC K11L = ToReal(0.); + CCTK_REAL_VEC K11L = ToReal(0); - CCTK_REAL_VEC K12L = ToReal(0.); + CCTK_REAL_VEC K12L = ToReal(0); - CCTK_REAL_VEC K13L = ToReal(0.); + CCTK_REAL_VEC K13L = ToReal(0); - CCTK_REAL_VEC K22L = ToReal(0.); + CCTK_REAL_VEC K22L = ToReal(0); - CCTK_REAL_VEC K23L = ToReal(0.); + CCTK_REAL_VEC K23L = ToReal(0); - CCTK_REAL_VEC K33L = ToReal(0.); + CCTK_REAL_VEC K33L = ToReal(0); - CCTK_REAL_VEC alphaL = ToReal(1.); + CCTK_REAL_VEC alphaL = ToReal(1); - CCTK_REAL_VEC beta1L = ToReal(0.); + CCTK_REAL_VEC beta1L = ToReal(0); - CCTK_REAL_VEC beta2L = ToReal(0.); + CCTK_REAL_VEC beta2L = ToReal(0); - CCTK_REAL_VEC beta3L = ToReal(0.); + CCTK_REAL_VEC beta3L = ToReal(0); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); |