diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-08-12 11:33:28 +0200 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-08-12 11:33:28 +0200 |
commit | 222a08bbeedce3eb4fef78bfa91c917279ad5f90 (patch) | |
tree | 98b7149713bb3e6ac0658e4905db94df275f2f65 /ML_BSSN/src/ML_BSSN_Minkowski.cc | |
parent | 9301a1b665598cafd094b7a653419de4cef77640 (diff) |
Regenerate code
Kranc commit 4707e0c3f03cc52d82f88735e602c52cf74f0cbe
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_Minkowski.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_Minkowski.cc | 63 |
1 files changed, 32 insertions, 31 deletions
diff --git a/ML_BSSN/src/ML_BSSN_Minkowski.cc b/ML_BSSN/src/ML_BSSN_Minkowski.cc index a41f5c2..05e4343 100644 --- a/ML_BSSN/src/ML_BSSN_Minkowski.cc +++ b/ML_BSSN/src/ML_BSSN_Minkowski.cc @@ -117,12 +117,12 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const 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); @@ -226,55 +226,56 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC phiL = IfThen(conformalMethod == 1,ToReal(1),ToReal(0)); + CCTK_REAL_VEC phiL = IfThen(conformalMethod == + 1.,ToReal(1.),ToReal(0.)); - CCTK_REAL_VEC gt11L = ToReal(1); + CCTK_REAL_VEC gt11L = ToReal(1.); - CCTK_REAL_VEC gt12L = ToReal(0); + CCTK_REAL_VEC gt12L = ToReal(0.); - CCTK_REAL_VEC gt13L = ToReal(0); + CCTK_REAL_VEC gt13L = ToReal(0.); - CCTK_REAL_VEC gt22L = ToReal(1); + CCTK_REAL_VEC gt22L = ToReal(1.); - CCTK_REAL_VEC gt23L = ToReal(0); + CCTK_REAL_VEC gt23L = ToReal(0.); - CCTK_REAL_VEC gt33L = ToReal(1); + CCTK_REAL_VEC gt33L = ToReal(1.); - CCTK_REAL_VEC trKL = ToReal(0); + CCTK_REAL_VEC trKL = ToReal(0.); - CCTK_REAL_VEC At11L = ToReal(0); + CCTK_REAL_VEC At11L = ToReal(0.); - CCTK_REAL_VEC At12L = ToReal(0); + CCTK_REAL_VEC At12L = ToReal(0.); - CCTK_REAL_VEC At13L = ToReal(0); + CCTK_REAL_VEC At13L = ToReal(0.); - CCTK_REAL_VEC At22L = ToReal(0); + CCTK_REAL_VEC At22L = ToReal(0.); - CCTK_REAL_VEC At23L = ToReal(0); + CCTK_REAL_VEC At23L = ToReal(0.); - CCTK_REAL_VEC At33L = ToReal(0); + CCTK_REAL_VEC At33L = ToReal(0.); - CCTK_REAL_VEC Xt1L = ToReal(0); + CCTK_REAL_VEC Xt1L = ToReal(0.); - CCTK_REAL_VEC Xt2L = ToReal(0); + CCTK_REAL_VEC Xt2L = ToReal(0.); - CCTK_REAL_VEC Xt3L = ToReal(0); + CCTK_REAL_VEC Xt3L = ToReal(0.); - CCTK_REAL_VEC alphaL = ToReal(1); + CCTK_REAL_VEC alphaL = ToReal(1.); - CCTK_REAL_VEC AL = ToReal(0); + CCTK_REAL_VEC AL = ToReal(0.); - 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.); - CCTK_REAL_VEC B1L = ToReal(0); + CCTK_REAL_VEC B1L = ToReal(0.); - CCTK_REAL_VEC B2L = ToReal(0); + CCTK_REAL_VEC B2L = ToReal(0.); - CCTK_REAL_VEC B3L = ToReal(0); + CCTK_REAL_VEC B3L = ToReal(0.); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); |