diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc')
| -rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc index 9a438cf..97ed004 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc @@ -289,7 +289,7 @@ static void ML_CCZ4_convertToADMBase_Body(const cGH* restrict const cctkGH, cons } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED e4phi = IfThen(conformalMethod == + CCTK_REAL_VEC e4phi CCTK_ATTRIBUTE_UNUSED = IfThen(conformalMethod == 1,kdiv(ToReal(1),kmul(phiL,phiL)),kexp(kmul(phiL,ToReal(4)))); gxxL = kmul(gt11L,e4phi); @@ -304,31 +304,31 @@ static void ML_CCZ4_convertToADMBase_Body(const cGH* restrict const cctkGH, cons gzzL = kmul(gt33L,e4phi); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxxL = + CCTK_REAL_VEC kxxL CCTK_ATTRIBUTE_UNUSED = kmadd(At11L,e4phi,kmul(gxxL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxyL = + CCTK_REAL_VEC kxyL CCTK_ATTRIBUTE_UNUSED = kmadd(At12L,e4phi,kmul(gxyL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxzL = + CCTK_REAL_VEC kxzL CCTK_ATTRIBUTE_UNUSED = kmadd(At13L,e4phi,kmul(gxzL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kyyL = + CCTK_REAL_VEC kyyL CCTK_ATTRIBUTE_UNUSED = kmadd(At22L,e4phi,kmul(gyyL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kyzL = + CCTK_REAL_VEC kyzL CCTK_ATTRIBUTE_UNUSED = kmadd(At23L,e4phi,kmul(gyzL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kzzL = + CCTK_REAL_VEC kzzL CCTK_ATTRIBUTE_UNUSED = kmadd(At33L,e4phi,kmul(gzzL,kmul(trKL,ToReal(0.333333333333333333333333333333)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED alpL = alphaL; + CCTK_REAL_VEC alpL CCTK_ATTRIBUTE_UNUSED = alphaL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betaxL = beta1L; + CCTK_REAL_VEC betaxL CCTK_ATTRIBUTE_UNUSED = beta1L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betayL = beta2L; + CCTK_REAL_VEC betayL CCTK_ATTRIBUTE_UNUSED = beta2L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betazL = beta3L; + CCTK_REAL_VEC betazL CCTK_ATTRIBUTE_UNUSED = beta3L; /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); |