diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc | 40 |
1 files changed, 20 insertions, 20 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc index 660c154..9f3ef7a 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc @@ -556,11 +556,11 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(const cGH* restrict const } /* Calculate temporaries and grid functions */ - ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir1 = kisgn(beta1L); + ptrdiff_t dir1 CCTK_ATTRIBUTE_UNUSED = kisgn(beta1L); - ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir2 = kisgn(beta2L); + ptrdiff_t dir2 CCTK_ATTRIBUTE_UNUSED = kisgn(beta2L); - ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir3 = kisgn(beta3L); + ptrdiff_t dir3 CCTK_ATTRIBUTE_UNUSED = kisgn(beta3L); CCTK_REAL_VEC JacPDstandardNth1alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDstandardNth1gt11 CCTK_ATTRIBUTE_UNUSED; @@ -856,48 +856,48 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(const cGH* restrict const JacPDupwindNthSymm3beta3 = PDupwindNthSymm3beta3; } - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED detgt = ToReal(1); + CCTK_REAL_VEC detgt CCTK_ATTRIBUTE_UNUSED = ToReal(1); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu11 = + CCTK_REAL_VEC gtu11 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt22L,gt33L,kmul(gt23L,gt23L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu12 = + CCTK_REAL_VEC gtu12 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt13L,gt23L,kmul(gt12L,gt33L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu13 = + CCTK_REAL_VEC gtu13 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt12L,gt23L,kmul(gt13L,gt22L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu22 = + CCTK_REAL_VEC gtu22 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt11L,gt33L,kmul(gt13L,gt13L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu23 = + CCTK_REAL_VEC gtu23 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt12L,gt13L,kmul(gt11L,gt23L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu33 = + CCTK_REAL_VEC gtu33 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt11L,gt22L,kmul(gt12L,gt12L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED em4phi = IfThen(conformalMethod == + CCTK_REAL_VEC em4phi CCTK_ATTRIBUTE_UNUSED = IfThen(conformalMethod == 1,kmul(phiL,phiL),kexp(kmul(phiL,ToReal(-4)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED eta = + CCTK_REAL_VEC eta CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(SpatialBetaDriverRadius),kfmax(rL,ToReal(SpatialBetaDriverRadius))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED theta = + CCTK_REAL_VEC theta CCTK_ATTRIBUTE_UNUSED = kfmin(ToReal(1),kexp(knmsub(rL,ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt1 = + CCTK_REAL_VEC Ddetgt1 CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu22,JacPDstandardNth1gt22,kmadd(gtu33,JacPDstandardNth1gt33,kmul(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmul(gtu23,JacPDstandardNth1gt23))),ToReal(2))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt2 = + CCTK_REAL_VEC Ddetgt2 CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth2gt11,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmul(kmadd(gtu12,JacPDstandardNth2gt12,kmadd(gtu13,JacPDstandardNth2gt13,kmul(gtu23,JacPDstandardNth2gt23))),ToReal(2))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt3 = + CCTK_REAL_VEC Ddetgt3 CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtalpL = + CCTK_REAL_VEC dtalpL CCTK_ATTRIBUTE_UNUSED = kmsub(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff),kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(ksub(AL,trKL),ToReal(LapseACoeff),trKL)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetaxL = + CCTK_REAL_VEC dtbetaxL CCTK_ATTRIBUTE_UNUSED = kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod == 1,kdiv(ToReal(1),phiL),ToReal(-2))))))))),ToReal(-2),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2))),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod @@ -907,7 +907,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(const cGH* restrict const 1,kdiv(ToReal(1),phiL),ToReal(-2))))))))),ToReal(-2),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2))))))))),kmul(theta,kmul(kadd(Xt1L,kmadd(ksub(B1L,Xt1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1 + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetayL = + CCTK_REAL_VEC dtbetayL CCTK_ATTRIBUTE_UNUSED = kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta2,kmadd(beta2L,JacPDupwindNthAnti2beta2,kmadd(beta3L,JacPDupwindNthAnti3beta2,kmadd(JacPDupwindNthSymm1beta2,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta2,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta2,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod == 1,kdiv(ToReal(1),phiL),ToReal(-2))))))))),ToReal(-2),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2))),kmadd(gtu22,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod @@ -917,7 +917,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(const cGH* restrict const 1,kdiv(ToReal(1),phiL),ToReal(-2))))))))),ToReal(-2),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2))))))))),kmul(theta,kmul(kadd(Xt2L,kmadd(ksub(B2L,Xt2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1 + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetazL = + CCTK_REAL_VEC dtbetazL CCTK_ATTRIBUTE_UNUSED = kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta3,kmadd(beta2L,JacPDupwindNthAnti2beta3,kmadd(beta3L,JacPDupwindNthAnti3beta3,kmadd(JacPDupwindNthSymm1beta3,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta3,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta3,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod == 1,kdiv(ToReal(1),phiL),ToReal(-2))))))))),ToReal(-2),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2))),kmadd(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod |