diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc')
| -rw-r--r-- | ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc | 64 |
1 files changed, 32 insertions, 32 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc b/ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc index 8afd097..de9fe5d 100644 --- a/ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc +++ b/ML_CCZ4/src/ML_CCZ4_RHSStaticBoundary.cc @@ -156,12 +156,12 @@ static void ML_CCZ4_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in 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 p1odx = kdiv(ToReal(1),dx); - CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx)); - CCTK_REAL_VEC const p1ody = kdiv(ToReal(1),dy); - CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy)); - CCTK_REAL_VEC const p1odz = kdiv(ToReal(1),dz); - CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz)); + CCTK_REAL_VEC const p1odx = kdiv(ToReal(1.),dx); + CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1.),kmul(dx,dx)); + CCTK_REAL_VEC const p1ody = kdiv(ToReal(1.),dy); + CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1.),kmul(dy,dy)); + CCTK_REAL_VEC const p1odz = kdiv(ToReal(1.),dz); + CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1.),kmul(dz,dz)); CCTK_REAL_VEC const pm1o120dx = kdiv(ToReal(-0.00833333333333333333333333333333),dx); CCTK_REAL_VEC const pm1o120dy = kdiv(ToReal(-0.00833333333333333333333333333333),dy); CCTK_REAL_VEC const pm1o120dz = kdiv(ToReal(-0.00833333333333333333333333333333),dz); @@ -265,57 +265,57 @@ static void ML_CCZ4_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC phirhsL = ToReal(0); + CCTK_REAL_VEC phirhsL = ToReal(0.); - CCTK_REAL_VEC gt11rhsL = ToReal(0); + CCTK_REAL_VEC gt11rhsL = ToReal(0.); - CCTK_REAL_VEC gt12rhsL = ToReal(0); + CCTK_REAL_VEC gt12rhsL = ToReal(0.); - CCTK_REAL_VEC gt13rhsL = ToReal(0); + CCTK_REAL_VEC gt13rhsL = ToReal(0.); - CCTK_REAL_VEC gt22rhsL = ToReal(0); + CCTK_REAL_VEC gt22rhsL = ToReal(0.); - CCTK_REAL_VEC gt23rhsL = ToReal(0); + CCTK_REAL_VEC gt23rhsL = ToReal(0.); - CCTK_REAL_VEC gt33rhsL = ToReal(0); + CCTK_REAL_VEC gt33rhsL = ToReal(0.); - CCTK_REAL_VEC trKrhsL = ToReal(0); + CCTK_REAL_VEC trKrhsL = ToReal(0.); - CCTK_REAL_VEC At11rhsL = ToReal(0); + CCTK_REAL_VEC At11rhsL = ToReal(0.); - CCTK_REAL_VEC At12rhsL = ToReal(0); + CCTK_REAL_VEC At12rhsL = ToReal(0.); - CCTK_REAL_VEC At13rhsL = ToReal(0); + CCTK_REAL_VEC At13rhsL = ToReal(0.); - CCTK_REAL_VEC At22rhsL = ToReal(0); + CCTK_REAL_VEC At22rhsL = ToReal(0.); - CCTK_REAL_VEC At23rhsL = ToReal(0); + CCTK_REAL_VEC At23rhsL = ToReal(0.); - CCTK_REAL_VEC At33rhsL = ToReal(0); + CCTK_REAL_VEC At33rhsL = ToReal(0.); - CCTK_REAL_VEC Xt1rhsL = ToReal(0); + CCTK_REAL_VEC Xt1rhsL = ToReal(0.); - CCTK_REAL_VEC Xt2rhsL = ToReal(0); + CCTK_REAL_VEC Xt2rhsL = ToReal(0.); - CCTK_REAL_VEC Xt3rhsL = ToReal(0); + CCTK_REAL_VEC Xt3rhsL = ToReal(0.); - CCTK_REAL_VEC alpharhsL = ToReal(0); + CCTK_REAL_VEC alpharhsL = ToReal(0.); - CCTK_REAL_VEC ArhsL = ToReal(0); + CCTK_REAL_VEC ArhsL = ToReal(0.); - CCTK_REAL_VEC beta1rhsL = ToReal(0); + CCTK_REAL_VEC beta1rhsL = ToReal(0.); - CCTK_REAL_VEC beta2rhsL = ToReal(0); + CCTK_REAL_VEC beta2rhsL = ToReal(0.); - CCTK_REAL_VEC beta3rhsL = ToReal(0); + CCTK_REAL_VEC beta3rhsL = ToReal(0.); - CCTK_REAL_VEC B1rhsL = ToReal(0); + CCTK_REAL_VEC B1rhsL = ToReal(0.); - CCTK_REAL_VEC B2rhsL = ToReal(0); + CCTK_REAL_VEC B2rhsL = ToReal(0.); - CCTK_REAL_VEC B3rhsL = ToReal(0); + CCTK_REAL_VEC B3rhsL = ToReal(0.); - CCTK_REAL_VEC ThetarhsL = ToReal(0); + CCTK_REAL_VEC ThetarhsL = ToReal(0.); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); |