diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-06-11 04:00:49 +0200 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-06-11 04:34:23 +0200 |
commit | f1bd0bd8c76cbaeab3ba911e428a00db9f352c70 (patch) | |
tree | e62b7c30c74c3a4a5f6cbd0cc10096907f056550 /ML_BSSN_O8 | |
parent | 158189c6675c8f03869e231f7cd36b4f96e4a75c (diff) |
Regenerate code after OutputForm change in Kranc
Kranc commit 216385dbca32325749f88b9dd3cda85705d800a0
Diffstat (limited to 'ML_BSSN_O8')
-rw-r--r-- | ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBaseGamma.cc | 48 |
1 files changed, 26 insertions, 22 deletions
diff --git a/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBaseGamma.cc b/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBaseGamma.cc index 54b6f5e..1594e0b 100644 --- a/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBaseGamma.cc +++ b/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBaseGamma.cc @@ -288,32 +288,36 @@ static void ML_BSSN_O8_convertFromADMBaseGamma_Body(cGH const * restrict const c CCTK_REAL Xt3L = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu12 + Gt313*gtu13 + Gt323*gtu23) + Gt333*gtu33; - CCTK_REAL AL = IfThen(LapseACoeff != 0,-3*(9*dtalpL + - LapseAdvectionCoeff*(-(beta1L*PDupwindNthAnti1alpha) - - beta2L*PDupwindNthAnti2alpha - beta3L*PDupwindNthAnti3alpha - - PDupwindNthSymm1alpha*Abs(beta1L) - PDupwindNthSymm2alpha*Abs(beta2L) - - PDupwindNthSymm3alpha*Abs(beta3L)))*INV(harmonicF)*pow(alphaL,-harmonicN),0); + CCTK_REAL AL = IfThen(LapseACoeff != + 0,-3*INV(ToReal(harmonicF))*pow(alphaL,-ToReal(harmonicN))*(9*dtalpL - + (beta1L*PDupwindNthAnti1alpha + beta2L*PDupwindNthAnti2alpha + + beta3L*PDupwindNthAnti3alpha + PDupwindNthSymm1alpha*Abs(beta1L) + + PDupwindNthSymm2alpha*Abs(beta2L) + + PDupwindNthSymm3alpha*Abs(beta3L))*ToReal(LapseAdvectionCoeff)),0); CCTK_REAL theta = fmin(1,exp(1 - rL*INV(ToReal(SpatialShiftGammaCoeffRadius)))); - CCTK_REAL B1L = IfThen(ShiftBCoeff*ShiftGammaCoeff != 0,3*(9*dtbetaxL - + ShiftAdvectionCoeff*(-(beta1L*PDupwindNthAnti1beta1) - - beta2L*PDupwindNthAnti2beta1 - beta3L*PDupwindNthAnti3beta1 - - PDupwindNthSymm1beta1*Abs(beta1L) - PDupwindNthSymm2beta1*Abs(beta2L) - - PDupwindNthSymm3beta1*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); - - CCTK_REAL B2L = IfThen(ShiftBCoeff*ShiftGammaCoeff != 0,3*(9*dtbetayL - + ShiftAdvectionCoeff*(-(beta1L*PDupwindNthAnti1beta2) - - beta2L*PDupwindNthAnti2beta2 - beta3L*PDupwindNthAnti3beta2 - - PDupwindNthSymm1beta2*Abs(beta1L) - PDupwindNthSymm2beta2*Abs(beta2L) - - PDupwindNthSymm3beta2*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); - - CCTK_REAL B3L = IfThen(ShiftBCoeff*ShiftGammaCoeff != 0,3*(9*dtbetazL - + ShiftAdvectionCoeff*(-(beta1L*PDupwindNthAnti1beta3) - - beta2L*PDupwindNthAnti2beta3 - beta3L*PDupwindNthAnti3beta3 - - PDupwindNthSymm1beta3*Abs(beta1L) - PDupwindNthSymm2beta3*Abs(beta2L) - - PDupwindNthSymm3beta3*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); + CCTK_REAL B1L = IfThen(ShiftBCoeff*ShiftGammaCoeff != + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetaxL - + 3*(beta1L*PDupwindNthAnti1beta1 + beta2L*PDupwindNthAnti2beta1 + + beta3L*PDupwindNthAnti3beta1 + PDupwindNthSymm1beta1*Abs(beta1L) + + PDupwindNthSymm2beta1*Abs(beta2L) + + PDupwindNthSymm3beta1*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); + + CCTK_REAL B2L = IfThen(ShiftBCoeff*ShiftGammaCoeff != + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetayL - + 3*(beta1L*PDupwindNthAnti1beta2 + beta2L*PDupwindNthAnti2beta2 + + beta3L*PDupwindNthAnti3beta2 + PDupwindNthSymm1beta2*Abs(beta1L) + + PDupwindNthSymm2beta2*Abs(beta2L) + + PDupwindNthSymm3beta2*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); + + CCTK_REAL B3L = IfThen(ShiftBCoeff*ShiftGammaCoeff != + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetazL - + 3*(beta1L*PDupwindNthAnti1beta3 + beta2L*PDupwindNthAnti2beta3 + + beta3L*PDupwindNthAnti3beta3 + PDupwindNthSymm1beta3*Abs(beta1L) + + PDupwindNthSymm2beta3*Abs(beta2L) + + PDupwindNthSymm3beta3*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); /* Copy local copies back to grid functions */ A[index] = AL; |