aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN_O8
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2011-06-11 04:00:49 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2011-06-11 04:34:23 +0200
commitf1bd0bd8c76cbaeab3ba911e428a00db9f352c70 (patch)
treee62b7c30c74c3a4a5f6cbd0cc10096907f056550 /ML_BSSN_O8
parent158189c6675c8f03869e231f7cd36b4f96e4a75c (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.cc48
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;