aboutsummaryrefslogtreecommitdiff
path: root/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc
diff options
context:
space:
mode:
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc')
-rw-r--r--ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc180
1 files changed, 94 insertions, 86 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc
index 0d06a21..61b8fb2 100644
--- a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc
+++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc
@@ -18,10 +18,14 @@
/* Define macros used in calculations */
#define INITVALUE (42)
-#define QAD(x) (SQR(SQR(x)))
+#define ScalarINV(x) ((CCTK_REAL)1.0 / (x))
+#define ScalarSQR(x) ((x) * (x))
+#define ScalarCUB(x) ((x) * ScalarSQR(x))
+#define ScalarQAD(x) (ScalarSQR(ScalarSQR(x)))
#define INV(x) (kdiv(ToReal(1.0),x))
#define SQR(x) (kmul(x,x))
#define CUB(x) (kmul(x,SQR(x)))
+#define QAD(x) (SQR(SQR(x)))
extern "C" void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_SelectBCs(CCTK_ARGUMENTS)
{
@@ -71,84 +75,84 @@ static void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
CCTK_REAL_VEC const hdzi = kmul(ToReal(0.5), dzi);
/* Initialize predefined quantities */
- CCTK_REAL_VEC const p1o1024dx = kmul(INV(dx),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o1024dy = kmul(INV(dy),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o1024dz = kmul(INV(dz),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o120dx = kmul(INV(dx),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o120dy = kmul(INV(dy),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o120dz = kmul(INV(dz),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dx = kmul(INV(dx),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dy = kmul(INV(dy),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dz = kmul(INV(dz),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o144dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o144dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o144dydz = kmul(INV(kmul(dy,dz)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o1680dx = kmul(INV(dx),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o1680dy = kmul(INV(dy),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o1680dz = kmul(INV(dz),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o16dx = kmul(INV(dx),ToReal(0.0625));
- CCTK_REAL_VEC const p1o16dy = kmul(INV(dy),ToReal(0.0625));
- CCTK_REAL_VEC const p1o16dz = kmul(INV(dz),ToReal(0.0625));
- CCTK_REAL_VEC const p1o180dx2 = kmul(INV(SQR(dx)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o180dy2 = kmul(INV(SQR(dy)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o180dz2 = kmul(INV(SQR(dz)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o24dx = kmul(INV(dx),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o24dy = kmul(INV(dy),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o24dz = kmul(INV(dz),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o256dx = kmul(INV(dx),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o256dy = kmul(INV(dy),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o256dz = kmul(INV(dz),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o2dx = kmul(INV(dx),ToReal(0.5));
- CCTK_REAL_VEC const p1o2dy = kmul(INV(dy),ToReal(0.5));
- CCTK_REAL_VEC const p1o2dz = kmul(INV(dz),ToReal(0.5));
- CCTK_REAL_VEC const p1o3600dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o3600dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o3600dydz = kmul(INV(kmul(dy,dz)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o4dx = kmul(INV(dx),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dy = kmul(INV(dy),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dydz = kmul(INV(kmul(dy,dz)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dz = kmul(INV(dz),ToReal(0.25));
- CCTK_REAL_VEC const p1o5040dx2 = kmul(INV(SQR(dx)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o5040dy2 = kmul(INV(SQR(dy)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o5040dz2 = kmul(INV(SQR(dz)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o560dx = kmul(INV(dx),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o560dy = kmul(INV(dy),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o560dz = kmul(INV(dz),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o60dx = kmul(INV(dx),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o60dy = kmul(INV(dy),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o60dz = kmul(INV(dz),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o64dx = kmul(INV(dx),ToReal(0.015625));
- CCTK_REAL_VEC const p1o64dy = kmul(INV(dy),ToReal(0.015625));
- CCTK_REAL_VEC const p1o64dz = kmul(INV(dz),ToReal(0.015625));
- CCTK_REAL_VEC const p1o705600dxdy = kmul(INV(kmul(dx,dy)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o705600dxdz = kmul(INV(kmul(dx,dz)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o705600dydz = kmul(INV(kmul(dy,dz)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o840dx = kmul(INV(dx),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1o840dy = kmul(INV(dy),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1o840dz = kmul(INV(dz),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1odx = INV(dx);
- CCTK_REAL_VEC const p1odx2 = INV(SQR(dx));
- CCTK_REAL_VEC const p1ody = INV(dy);
- CCTK_REAL_VEC const p1ody2 = INV(SQR(dy));
- CCTK_REAL_VEC const p1odz = INV(dz);
- CCTK_REAL_VEC const p1odz2 = INV(SQR(dz));
- CCTK_REAL_VEC const pm1o120dx = kmul(INV(dx),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o120dy = kmul(INV(dy),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o120dz = kmul(INV(dz),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dx2 = kmul(INV(SQR(dx)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dy2 = kmul(INV(SQR(dy)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dz2 = kmul(INV(SQR(dz)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o2dx = kmul(INV(dx),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o2dy = kmul(INV(dy),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o2dz = kmul(INV(dz),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o4dx = kmul(INV(dx),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o4dy = kmul(INV(dy),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o4dz = kmul(INV(dz),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o60dx = kmul(INV(dx),ToReal(-0.0166666666666666666666666666667));
- CCTK_REAL_VEC const pm1o60dy = kmul(INV(dy),ToReal(-0.0166666666666666666666666666667));
- CCTK_REAL_VEC const pm1o60dz = kmul(INV(dz),ToReal(-0.0166666666666666666666666666667));
+ CCTK_REAL_VEC const p1o1024dx = kdiv(ToReal(0.0009765625),dx);
+ CCTK_REAL_VEC const p1o1024dy = kdiv(ToReal(0.0009765625),dy);
+ CCTK_REAL_VEC const p1o1024dz = kdiv(ToReal(0.0009765625),dz);
+ CCTK_REAL_VEC const p1o120dx = kdiv(ToReal(0.00833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const p1o120dy = kdiv(ToReal(0.00833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const p1o120dz = kdiv(ToReal(0.00833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const p1o12dx = kdiv(ToReal(0.0833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const p1o12dy = kdiv(ToReal(0.0833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const p1o12dz = kdiv(ToReal(0.0833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const p1o144dxdy = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o144dxdz = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o144dydz = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o1680dx = kdiv(ToReal(0.000595238095238095238095238095238),dx);
+ CCTK_REAL_VEC const p1o1680dy = kdiv(ToReal(0.000595238095238095238095238095238),dy);
+ CCTK_REAL_VEC const p1o1680dz = kdiv(ToReal(0.000595238095238095238095238095238),dz);
+ CCTK_REAL_VEC const p1o16dx = kdiv(ToReal(0.0625),dx);
+ CCTK_REAL_VEC const p1o16dy = kdiv(ToReal(0.0625),dy);
+ CCTK_REAL_VEC const p1o16dz = kdiv(ToReal(0.0625),dz);
+ CCTK_REAL_VEC const p1o180dx2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dx,dx));
+ CCTK_REAL_VEC const p1o180dy2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dy,dy));
+ CCTK_REAL_VEC const p1o180dz2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dz,dz));
+ CCTK_REAL_VEC const p1o24dx = kdiv(ToReal(0.0416666666666666666666666666667),dx);
+ CCTK_REAL_VEC const p1o24dy = kdiv(ToReal(0.0416666666666666666666666666667),dy);
+ CCTK_REAL_VEC const p1o24dz = kdiv(ToReal(0.0416666666666666666666666666667),dz);
+ CCTK_REAL_VEC const p1o256dx = kdiv(ToReal(0.00390625),dx);
+ CCTK_REAL_VEC const p1o256dy = kdiv(ToReal(0.00390625),dy);
+ CCTK_REAL_VEC const p1o256dz = kdiv(ToReal(0.00390625),dz);
+ CCTK_REAL_VEC const p1o2dx = kdiv(ToReal(0.5),dx);
+ CCTK_REAL_VEC const p1o2dy = kdiv(ToReal(0.5),dy);
+ CCTK_REAL_VEC const p1o2dz = kdiv(ToReal(0.5),dz);
+ CCTK_REAL_VEC const p1o3600dxdy = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o3600dxdz = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o3600dydz = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o4dx = kdiv(ToReal(0.25),dx);
+ CCTK_REAL_VEC const p1o4dxdy = kdiv(ToReal(0.25),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o4dxdz = kdiv(ToReal(0.25),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o4dy = kdiv(ToReal(0.25),dy);
+ CCTK_REAL_VEC const p1o4dydz = kdiv(ToReal(0.25),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o4dz = kdiv(ToReal(0.25),dz);
+ CCTK_REAL_VEC const p1o5040dx2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dx,dx));
+ CCTK_REAL_VEC const p1o5040dy2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dy,dy));
+ CCTK_REAL_VEC const p1o5040dz2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dz,dz));
+ CCTK_REAL_VEC const p1o560dx = kdiv(ToReal(0.00178571428571428571428571428571),dx);
+ CCTK_REAL_VEC const p1o560dy = kdiv(ToReal(0.00178571428571428571428571428571),dy);
+ CCTK_REAL_VEC const p1o560dz = kdiv(ToReal(0.00178571428571428571428571428571),dz);
+ CCTK_REAL_VEC const p1o60dx = kdiv(ToReal(0.0166666666666666666666666666667),dx);
+ CCTK_REAL_VEC const p1o60dy = kdiv(ToReal(0.0166666666666666666666666666667),dy);
+ CCTK_REAL_VEC const p1o60dz = kdiv(ToReal(0.0166666666666666666666666666667),dz);
+ CCTK_REAL_VEC const p1o64dx = kdiv(ToReal(0.015625),dx);
+ CCTK_REAL_VEC const p1o64dy = kdiv(ToReal(0.015625),dy);
+ CCTK_REAL_VEC const p1o64dz = kdiv(ToReal(0.015625),dz);
+ CCTK_REAL_VEC const p1o705600dxdy = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o705600dxdz = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o705600dydz = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o840dx = kdiv(ToReal(0.00119047619047619047619047619048),dx);
+ CCTK_REAL_VEC const p1o840dy = kdiv(ToReal(0.00119047619047619047619047619048),dy);
+ CCTK_REAL_VEC const p1o840dz = kdiv(ToReal(0.00119047619047619047619047619048),dz);
+ CCTK_REAL_VEC const p1odx = kdiv(ToReal(1),dx);
+ CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx));
+ CCTK_REAL_VEC const p1ody = kdiv(ToReal(1),dy);
+ CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy));
+ CCTK_REAL_VEC const p1odz = kdiv(ToReal(1),dz);
+ CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz));
+ CCTK_REAL_VEC const pm1o120dx = kdiv(ToReal(-0.00833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const pm1o120dy = kdiv(ToReal(-0.00833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const pm1o120dz = kdiv(ToReal(-0.00833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const pm1o12dx2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx));
+ CCTK_REAL_VEC const pm1o12dy2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy));
+ CCTK_REAL_VEC const pm1o12dz2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz));
+ CCTK_REAL_VEC const pm1o2dx = kdiv(ToReal(-0.5),dx);
+ CCTK_REAL_VEC const pm1o2dy = kdiv(ToReal(-0.5),dy);
+ CCTK_REAL_VEC const pm1o2dz = kdiv(ToReal(-0.5),dz);
+ CCTK_REAL_VEC const pm1o4dx = kdiv(ToReal(-0.25),dx);
+ CCTK_REAL_VEC const pm1o4dy = kdiv(ToReal(-0.25),dy);
+ CCTK_REAL_VEC const pm1o4dz = kdiv(ToReal(-0.25),dz);
+ CCTK_REAL_VEC const pm1o60dx = kdiv(ToReal(-0.0166666666666666666666666666667),dx);
+ CCTK_REAL_VEC const pm1o60dy = kdiv(ToReal(-0.0166666666666666666666666666667),dy);
+ CCTK_REAL_VEC const pm1o60dz = kdiv(ToReal(-0.0166666666666666666666666666667),dz);
/* Jacobian variable pointers */
bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map)
@@ -207,7 +211,7 @@ static void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
#pragma omp parallel
LC_LOOP3VEC(ML_CCZ4_convertToADMBaseDtLapseShiftBoundary,
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
- cctk_lsh[0],cctk_lsh[1],cctk_lsh[2],
+ cctk_ash[0],cctk_ash[1],cctk_ash[2],
CCTK_REAL_VEC_SIZE)
{
ptrdiff_t const index = di*i + dj*j + dk*k;
@@ -252,13 +256,14 @@ static void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
/* Calculate temporaries and grid functions */
CCTK_REAL_VEC eta =
- kmul(INV(kfmax(rL,ToReal(SpatialBetaDriverRadius))),ToReal(SpatialBetaDriverRadius));
+ kdiv(ToReal(SpatialBetaDriverRadius),kfmax(rL,ToReal(SpatialBetaDriverRadius)));
CCTK_REAL_VEC theta =
- kfmin(ToReal(1),kexp(knmsub(rL,INV(ToReal(SpatialShiftGammaCoeffRadius)),ToReal(1))));
+ kfmin(ToReal(1),kexp(knmsub(rL,ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1))));
CCTK_REAL_VEC dtalpL =
- kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(AL,ToReal(LapseACoeff),kmul(kmsub(ThetaL,ToReal(2),trKL),kadd(ToReal(-1),ToReal(LapseACoeff)))))));
+ kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(AL,ToReal(LapseACoeff),kmul(kmsub(ThetaL,ToReal(2),trKL),ToReal(-1
+ + LapseACoeff))))));
CCTK_REAL_VEC dtbetaxL;
CCTK_REAL_VEC dtbetayL;
@@ -275,13 +280,16 @@ static void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
else
{
dtbetaxL =
- kmul(theta,kmul(kadd(Xt1L,kmadd(ksub(B1L,Xt1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,kmul(ToReal(BetaDriver),kadd(ToReal(-1),ToReal(ShiftBCoeff))))))),ToReal(ShiftGammaCoeff)));
+ kmul(theta,kmul(kadd(Xt1L,kmadd(ksub(B1L,Xt1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1
+ + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)));
dtbetayL =
- kmul(theta,kmul(kadd(Xt2L,kmadd(ksub(B2L,Xt2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,kmul(ToReal(BetaDriver),kadd(ToReal(-1),ToReal(ShiftBCoeff))))))),ToReal(ShiftGammaCoeff)));
+ kmul(theta,kmul(kadd(Xt2L,kmadd(ksub(B2L,Xt2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1
+ + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)));
dtbetazL =
- kmul(theta,kmul(kadd(Xt3L,kmadd(ksub(B3L,Xt3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,kmul(ToReal(BetaDriver),kadd(ToReal(-1),ToReal(ShiftBCoeff))))))),ToReal(ShiftGammaCoeff)));
+ kmul(theta,kmul(kadd(Xt3L,kmadd(ksub(B3L,Xt3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,ToReal(BetaDriver*(-1
+ + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)));
}
/* Copy local copies back to grid functions */
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-29 06:09:27 +0000