diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-11-19 17:24:58 +0100 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-11-19 17:24:58 +0100 |
commit | b19621c279744d4d3eff3b575f9b277419f07c7c (patch) | |
tree | d93a1d7b7841f65a52e27829bedf83a16b6584cd /ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc | |
parent | eaf08e4235c0f0d5dfcea5640feb3c51a055627a (diff) |
Regenerate code (Kranc commit 17deb8b4ff45a8a3cd1645cb626f6d33a5e17beb)
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc index ed3b2de..8c39ef8 100644 --- a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc +++ b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc @@ -831,10 +831,10 @@ static void ML_BSSN_convertFromADMBaseGamma_Body(cGH const * restrict const cctk kmadd(Gt311,gtu11,kmadd(Gt322,gtu22,kmadd(Gt333,gtu33,kmul(kmadd(Gt312,gtu12,kmadd(Gt313,gtu13,kmul(Gt323,gtu23))),ToReal(2))))); CCTK_REAL_VEC AL = IfThen(LapseACoeff != - 0,kneg(kdiv(kmul(kpow(alphaL,-harmonicN),knmsub(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff),dtalpL)),ToReal(harmonicF))),ToReal(0)); + 0,kneg(kmul(kmul(kpow(alphaL,-harmonicN),knmsub(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff),dtalpL)),ToReal(ScalarINV(harmonicF)))),ToReal(0)); CCTK_REAL_VEC theta = - kfmin(ToReal(1),kexp(ksub(ToReal(1),kdiv(rL,ToReal(SpatialShiftGammaCoeffRadius))))); + kfmin(ToReal(1),kexp(knmsub(rL,ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1)))); CCTK_REAL_VEC B1L; CCTK_REAL_VEC B2L; |