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_MP | |
parent | 158189c6675c8f03869e231f7cd36b4f96e4a75c (diff) |
Regenerate code after OutputForm change in Kranc
Kranc commit 216385dbca32325749f88b9dd3cda85705d800a0
Diffstat (limited to 'ML_BSSN_MP')
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.cc | 77 |
1 files changed, 37 insertions, 40 deletions
diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.cc b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.cc index e9ed908..5862ff9 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.cc +++ b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.cc @@ -378,55 +378,52 @@ static void ML_BSSN_MP_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*J11L + beta2L*J12L + - beta3L*J13L)*PDupwindNthAnti1alpha + (beta1L*J21L + beta2L*J22L + - beta3L*J23L)*PDupwindNthAnti2alpha + (beta1L*J31L + beta2L*J32L + - beta3L*J33L)*PDupwindNthAnti3alpha + (J11L*PDupwindNthSymm1alpha + - J21L*PDupwindNthSymm2alpha + J31L*PDupwindNthSymm3alpha)*Abs(beta1L) + - (J12L*PDupwindNthSymm1alpha + J22L*PDupwindNthSymm2alpha + - J32L*PDupwindNthSymm3alpha)*Abs(beta2L) + (J13L*PDupwindNthSymm1alpha + - J23L*PDupwindNthSymm2alpha + - J33L*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*J11L + beta2L*J12L + beta3L*J13L)*PDupwindNthAnti1alpha + + (beta1L*J21L + beta2L*J22L + beta3L*J23L)*PDupwindNthAnti2alpha + + (beta1L*J31L + beta2L*J32L + beta3L*J33L)*PDupwindNthAnti3alpha + + (J11L*PDupwindNthSymm1alpha + J21L*PDupwindNthSymm2alpha + + J31L*PDupwindNthSymm3alpha)*Abs(beta1L) + (J12L*PDupwindNthSymm1alpha + + J22L*PDupwindNthSymm2alpha + J32L*PDupwindNthSymm3alpha)*Abs(beta2L) + + (J13L*PDupwindNthSymm1alpha + J23L*PDupwindNthSymm2alpha + + J33L*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*J11L + beta2L*J12L + - beta3L*J13L)*PDupwindNthAnti1beta1 + (beta2L*J22L + - beta3L*J23L)*PDupwindNthAnti2beta1 + (beta2L*J32L + - beta3L*J33L)*PDupwindNthAnti3beta1 + beta1L*(J21L*PDupwindNthAnti2beta1 - + J31L*PDupwindNthAnti3beta1) + (J11L*PDupwindNthSymm1beta1 + - J21L*PDupwindNthSymm2beta1 + J31L*PDupwindNthSymm3beta1)*Abs(beta1L) + - (J12L*PDupwindNthSymm1beta1 + J22L*PDupwindNthSymm2beta1 + - J32L*PDupwindNthSymm3beta1)*Abs(beta2L) + (J13L*PDupwindNthSymm1beta1 + - J23L*PDupwindNthSymm2beta1 + - J33L*PDupwindNthSymm3beta1)*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetaxL - + 3*((beta1L*J11L + beta2L*J12L + beta3L*J13L)*PDupwindNthAnti1beta1 + + (beta1L*J21L + beta2L*J22L + beta3L*J23L)*PDupwindNthAnti2beta1 + + (beta1L*J31L + beta2L*J32L + beta3L*J33L)*PDupwindNthAnti3beta1 + + (J11L*PDupwindNthSymm1beta1 + J21L*PDupwindNthSymm2beta1 + + J31L*PDupwindNthSymm3beta1)*Abs(beta1L) + (J12L*PDupwindNthSymm1beta1 + + J22L*PDupwindNthSymm2beta1 + J32L*PDupwindNthSymm3beta1)*Abs(beta2L) + + (J13L*PDupwindNthSymm1beta1 + J23L*PDupwindNthSymm2beta1 + + J33L*PDupwindNthSymm3beta1)*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); CCTK_REAL B2L = IfThen(ShiftBCoeff*ShiftGammaCoeff != - 0,-3*(-9*dtbetayL + ShiftAdvectionCoeff*((beta1L*J11L + beta2L*J12L + - beta3L*J13L)*PDupwindNthAnti1beta2 + (beta2L*J22L + - beta3L*J23L)*PDupwindNthAnti2beta2 + (beta2L*J32L + - beta3L*J33L)*PDupwindNthAnti3beta2 + beta1L*(J21L*PDupwindNthAnti2beta2 - + J31L*PDupwindNthAnti3beta2) + (J11L*PDupwindNthSymm1beta2 + - J21L*PDupwindNthSymm2beta2 + J31L*PDupwindNthSymm3beta2)*Abs(beta1L) + - (J12L*PDupwindNthSymm1beta2 + J22L*PDupwindNthSymm2beta2 + - J32L*PDupwindNthSymm3beta2)*Abs(beta2L) + (J13L*PDupwindNthSymm1beta2 + - J23L*PDupwindNthSymm2beta2 + - J33L*PDupwindNthSymm3beta2)*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetayL - + 3*((beta1L*J11L + beta2L*J12L + beta3L*J13L)*PDupwindNthAnti1beta2 + + (beta1L*J21L + beta2L*J22L + beta3L*J23L)*PDupwindNthAnti2beta2 + + (beta1L*J31L + beta2L*J32L + beta3L*J33L)*PDupwindNthAnti3beta2 + + (J11L*PDupwindNthSymm1beta2 + J21L*PDupwindNthSymm2beta2 + + J31L*PDupwindNthSymm3beta2)*Abs(beta1L) + (J12L*PDupwindNthSymm1beta2 + + J22L*PDupwindNthSymm2beta2 + J32L*PDupwindNthSymm3beta2)*Abs(beta2L) + + (J13L*PDupwindNthSymm1beta2 + J23L*PDupwindNthSymm2beta2 + + J33L*PDupwindNthSymm3beta2)*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); CCTK_REAL B3L = IfThen(ShiftBCoeff*ShiftGammaCoeff != - 0,-3*(-9*dtbetazL + ShiftAdvectionCoeff*((beta1L*J11L + beta2L*J12L + - beta3L*J13L)*PDupwindNthAnti1beta3 + (beta2L*J22L + - beta3L*J23L)*PDupwindNthAnti2beta3 + (beta2L*J32L + - beta3L*J33L)*PDupwindNthAnti3beta3 + beta1L*(J21L*PDupwindNthAnti2beta3 - + J31L*PDupwindNthAnti3beta3) + (J11L*PDupwindNthSymm1beta3 + - J21L*PDupwindNthSymm2beta3 + J31L*PDupwindNthSymm3beta3)*Abs(beta1L) + - (J12L*PDupwindNthSymm1beta3 + J22L*PDupwindNthSymm2beta3 + - J32L*PDupwindNthSymm3beta3)*Abs(beta2L) + (J13L*PDupwindNthSymm1beta3 + - J23L*PDupwindNthSymm2beta3 + - J33L*PDupwindNthSymm3beta3)*Abs(beta3L)))*INV(ShiftGammaCoeff)*INV(theta),0); + 0,INV(theta)*INV(ToReal(ShiftGammaCoeff))*(27*dtbetazL - + 3*((beta1L*J11L + beta2L*J12L + beta3L*J13L)*PDupwindNthAnti1beta3 + + (beta1L*J21L + beta2L*J22L + beta3L*J23L)*PDupwindNthAnti2beta3 + + (beta1L*J31L + beta2L*J32L + beta3L*J33L)*PDupwindNthAnti3beta3 + + (J11L*PDupwindNthSymm1beta3 + J21L*PDupwindNthSymm2beta3 + + J31L*PDupwindNthSymm3beta3)*Abs(beta1L) + (J12L*PDupwindNthSymm1beta3 + + J22L*PDupwindNthSymm2beta3 + J32L*PDupwindNthSymm3beta3)*Abs(beta2L) + + (J13L*PDupwindNthSymm1beta3 + J23L*PDupwindNthSymm2beta3 + + J33L*PDupwindNthSymm3beta3)*Abs(beta3L))*ToReal(ShiftAdvectionCoeff)),0); /* Copy local copies back to grid functions */ A[index] = AL; |