diff options
Diffstat (limited to 'src/psi4_calc_4th.cc')
-rw-r--r-- | src/psi4_calc_4th.cc | 229 |
1 files changed, 124 insertions, 105 deletions
diff --git a/src/psi4_calc_4th.cc b/src/psi4_calc_4th.cc index 284a741..3b79e7f 100644 --- a/src/psi4_calc_4th.cc +++ b/src/psi4_calc_4th.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 psi4_calc_4th_SelectBCs(CCTK_ARGUMENTS) { @@ -71,42 +75,42 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC const hdzi = kmul(ToReal(0.5), dzi); /* Initialize predefined quantities */ - 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 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 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 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 p1o4dydz = kmul(INV(kmul(dy,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 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 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 p1odx2 = INV(SQR(dx)); - CCTK_REAL_VEC const p1ody2 = INV(SQR(dy)); - CCTK_REAL_VEC const p1odz2 = INV(SQR(dz)); - 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 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 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 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 p1o4dxdy = kdiv(ToReal(0.25),kmul(dy,dx)); + CCTK_REAL_VEC const p1o4dxdz = kdiv(ToReal(0.25),kmul(dz,dx)); + CCTK_REAL_VEC const p1o4dydz = kdiv(ToReal(0.25),kmul(dz,dy)); + 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 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 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 p1odx2 = kdiv(ToReal(1),kmul(dx,dx)); + CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy)); + CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,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)); /* Jacobian variable pointers */ bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) @@ -165,7 +169,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, #pragma omp parallel LC_LOOP3VEC(psi4_calc_4th, 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; @@ -789,49 +793,49 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmadd(J13L,PDstandard4th1kyz,kmadd(J23L,PDstandard4th2kyz,kmul(J33L,PDstandard4th3kyz))); JacPDstandard4th11gxy = - kmadd(dJ111L,PDstandard4th1gxy,kmadd(dJ211L,PDstandard4th2gxy,kmadd(dJ311L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,SQR(J11L),kmadd(PDstandard4th22gxy,SQR(J21L),kmadd(PDstandard4th33gxy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy)),kmul(J21L,kmul(J31L,PDstandard4th23gxy))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard4th1gxy,kmadd(dJ211L,PDstandard4th2gxy,kmadd(dJ311L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J11L,J11L),kmadd(PDstandard4th22gxy,kmul(J21L,J21L),kmadd(PDstandard4th33gxy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy)),kmul(J21L,kmul(J31L,PDstandard4th23gxy))),ToReal(2)))))))); JacPDstandard4th11gxz = - kmadd(dJ111L,PDstandard4th1gxz,kmadd(dJ211L,PDstandard4th2gxz,kmadd(dJ311L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,SQR(J11L),kmadd(PDstandard4th22gxz,SQR(J21L),kmadd(PDstandard4th33gxz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz)),kmul(J21L,kmul(J31L,PDstandard4th23gxz))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard4th1gxz,kmadd(dJ211L,PDstandard4th2gxz,kmadd(dJ311L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J11L,J11L),kmadd(PDstandard4th22gxz,kmul(J21L,J21L),kmadd(PDstandard4th33gxz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz)),kmul(J21L,kmul(J31L,PDstandard4th23gxz))),ToReal(2)))))))); JacPDstandard4th11gyy = - kmadd(dJ111L,PDstandard4th1gyy,kmadd(dJ211L,PDstandard4th2gyy,kmadd(dJ311L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,SQR(J11L),kmadd(PDstandard4th22gyy,SQR(J21L),kmadd(PDstandard4th33gyy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy)),kmul(J21L,kmul(J31L,PDstandard4th23gyy))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard4th1gyy,kmadd(dJ211L,PDstandard4th2gyy,kmadd(dJ311L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,kmul(J11L,J11L),kmadd(PDstandard4th22gyy,kmul(J21L,J21L),kmadd(PDstandard4th33gyy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy)),kmul(J21L,kmul(J31L,PDstandard4th23gyy))),ToReal(2)))))))); JacPDstandard4th11gyz = - kmadd(dJ111L,PDstandard4th1gyz,kmadd(dJ211L,PDstandard4th2gyz,kmadd(dJ311L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,SQR(J11L),kmadd(PDstandard4th22gyz,SQR(J21L),kmadd(PDstandard4th33gyz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz)),kmul(J21L,kmul(J31L,PDstandard4th23gyz))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard4th1gyz,kmadd(dJ211L,PDstandard4th2gyz,kmadd(dJ311L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,kmul(J11L,J11L),kmadd(PDstandard4th22gyz,kmul(J21L,J21L),kmadd(PDstandard4th33gyz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz)),kmul(J21L,kmul(J31L,PDstandard4th23gyz))),ToReal(2)))))))); JacPDstandard4th11gzz = - kmadd(dJ111L,PDstandard4th1gzz,kmadd(dJ211L,PDstandard4th2gzz,kmadd(dJ311L,PDstandard4th3gzz,kmadd(PDstandard4th11gzz,SQR(J11L),kmadd(PDstandard4th22gzz,SQR(J21L),kmadd(PDstandard4th33gzz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz)),kmul(J21L,kmul(J31L,PDstandard4th23gzz))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard4th1gzz,kmadd(dJ211L,PDstandard4th2gzz,kmadd(dJ311L,PDstandard4th3gzz,kmadd(PDstandard4th11gzz,kmul(J11L,J11L),kmadd(PDstandard4th22gzz,kmul(J21L,J21L),kmadd(PDstandard4th33gzz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz)),kmul(J21L,kmul(J31L,PDstandard4th23gzz))),ToReal(2)))))))); JacPDstandard4th22gxx = - kmadd(dJ122L,PDstandard4th1gxx,kmadd(dJ222L,PDstandard4th2gxx,kmadd(dJ322L,PDstandard4th3gxx,kmadd(PDstandard4th11gxx,SQR(J12L),kmadd(PDstandard4th22gxx,SQR(J22L),kmadd(PDstandard4th33gxx,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmul(J22L,kmul(J32L,PDstandard4th23gxx))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard4th1gxx,kmadd(dJ222L,PDstandard4th2gxx,kmadd(dJ322L,PDstandard4th3gxx,kmadd(PDstandard4th11gxx,kmul(J12L,J12L),kmadd(PDstandard4th22gxx,kmul(J22L,J22L),kmadd(PDstandard4th33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmul(J22L,kmul(J32L,PDstandard4th23gxx))),ToReal(2)))))))); JacPDstandard4th22gxy = - kmadd(dJ122L,PDstandard4th1gxy,kmadd(dJ222L,PDstandard4th2gxy,kmadd(dJ322L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,SQR(J12L),kmadd(PDstandard4th22gxy,SQR(J22L),kmadd(PDstandard4th33gxy,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy)),kmul(J22L,kmul(J32L,PDstandard4th23gxy))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard4th1gxy,kmadd(dJ222L,PDstandard4th2gxy,kmadd(dJ322L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J12L,J12L),kmadd(PDstandard4th22gxy,kmul(J22L,J22L),kmadd(PDstandard4th33gxy,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy)),kmul(J22L,kmul(J32L,PDstandard4th23gxy))),ToReal(2)))))))); JacPDstandard4th22gxz = - kmadd(dJ122L,PDstandard4th1gxz,kmadd(dJ222L,PDstandard4th2gxz,kmadd(dJ322L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,SQR(J12L),kmadd(PDstandard4th22gxz,SQR(J22L),kmadd(PDstandard4th33gxz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz)),kmul(J22L,kmul(J32L,PDstandard4th23gxz))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard4th1gxz,kmadd(dJ222L,PDstandard4th2gxz,kmadd(dJ322L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J12L,J12L),kmadd(PDstandard4th22gxz,kmul(J22L,J22L),kmadd(PDstandard4th33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz)),kmul(J22L,kmul(J32L,PDstandard4th23gxz))),ToReal(2)))))))); JacPDstandard4th22gyz = - kmadd(dJ122L,PDstandard4th1gyz,kmadd(dJ222L,PDstandard4th2gyz,kmadd(dJ322L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,SQR(J12L),kmadd(PDstandard4th22gyz,SQR(J22L),kmadd(PDstandard4th33gyz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz)),kmul(J22L,kmul(J32L,PDstandard4th23gyz))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard4th1gyz,kmadd(dJ222L,PDstandard4th2gyz,kmadd(dJ322L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,kmul(J12L,J12L),kmadd(PDstandard4th22gyz,kmul(J22L,J22L),kmadd(PDstandard4th33gyz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz)),kmul(J22L,kmul(J32L,PDstandard4th23gyz))),ToReal(2)))))))); JacPDstandard4th22gzz = - kmadd(dJ122L,PDstandard4th1gzz,kmadd(dJ222L,PDstandard4th2gzz,kmadd(dJ322L,PDstandard4th3gzz,kmadd(PDstandard4th11gzz,SQR(J12L),kmadd(PDstandard4th22gzz,SQR(J22L),kmadd(PDstandard4th33gzz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz)),kmul(J22L,kmul(J32L,PDstandard4th23gzz))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard4th1gzz,kmadd(dJ222L,PDstandard4th2gzz,kmadd(dJ322L,PDstandard4th3gzz,kmadd(PDstandard4th11gzz,kmul(J12L,J12L),kmadd(PDstandard4th22gzz,kmul(J22L,J22L),kmadd(PDstandard4th33gzz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz)),kmul(J22L,kmul(J32L,PDstandard4th23gzz))),ToReal(2)))))))); JacPDstandard4th33gxx = - kmadd(dJ133L,PDstandard4th1gxx,kmadd(dJ233L,PDstandard4th2gxx,kmadd(dJ333L,PDstandard4th3gxx,kmadd(PDstandard4th11gxx,SQR(J13L),kmadd(PDstandard4th22gxx,SQR(J23L),kmadd(PDstandard4th33gxx,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmul(J23L,kmul(J33L,PDstandard4th23gxx))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard4th1gxx,kmadd(dJ233L,PDstandard4th2gxx,kmadd(dJ333L,PDstandard4th3gxx,kmadd(PDstandard4th11gxx,kmul(J13L,J13L),kmadd(PDstandard4th22gxx,kmul(J23L,J23L),kmadd(PDstandard4th33gxx,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmul(J23L,kmul(J33L,PDstandard4th23gxx))),ToReal(2)))))))); JacPDstandard4th33gxy = - kmadd(dJ133L,PDstandard4th1gxy,kmadd(dJ233L,PDstandard4th2gxy,kmadd(dJ333L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,SQR(J13L),kmadd(PDstandard4th22gxy,SQR(J23L),kmadd(PDstandard4th33gxy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmul(J23L,kmul(J33L,PDstandard4th23gxy))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard4th1gxy,kmadd(dJ233L,PDstandard4th2gxy,kmadd(dJ333L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J13L,J13L),kmadd(PDstandard4th22gxy,kmul(J23L,J23L),kmadd(PDstandard4th33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmul(J23L,kmul(J33L,PDstandard4th23gxy))),ToReal(2)))))))); JacPDstandard4th33gxz = - kmadd(dJ133L,PDstandard4th1gxz,kmadd(dJ233L,PDstandard4th2gxz,kmadd(dJ333L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,SQR(J13L),kmadd(PDstandard4th22gxz,SQR(J23L),kmadd(PDstandard4th33gxz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmul(J23L,kmul(J33L,PDstandard4th23gxz))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard4th1gxz,kmadd(dJ233L,PDstandard4th2gxz,kmadd(dJ333L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J13L,J13L),kmadd(PDstandard4th22gxz,kmul(J23L,J23L),kmadd(PDstandard4th33gxz,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmul(J23L,kmul(J33L,PDstandard4th23gxz))),ToReal(2)))))))); JacPDstandard4th33gyy = - kmadd(dJ133L,PDstandard4th1gyy,kmadd(dJ233L,PDstandard4th2gyy,kmadd(dJ333L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,SQR(J13L),kmadd(PDstandard4th22gyy,SQR(J23L),kmadd(PDstandard4th33gyy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmul(J23L,kmul(J33L,PDstandard4th23gyy))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard4th1gyy,kmadd(dJ233L,PDstandard4th2gyy,kmadd(dJ333L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,kmul(J13L,J13L),kmadd(PDstandard4th22gyy,kmul(J23L,J23L),kmadd(PDstandard4th33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmul(J23L,kmul(J33L,PDstandard4th23gyy))),ToReal(2)))))))); JacPDstandard4th33gyz = - kmadd(dJ133L,PDstandard4th1gyz,kmadd(dJ233L,PDstandard4th2gyz,kmadd(dJ333L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,SQR(J13L),kmadd(PDstandard4th22gyz,SQR(J23L),kmadd(PDstandard4th33gyz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmul(J23L,kmul(J33L,PDstandard4th23gyz))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard4th1gyz,kmadd(dJ233L,PDstandard4th2gyz,kmadd(dJ333L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,kmul(J13L,J13L),kmadd(PDstandard4th22gyz,kmul(J23L,J23L),kmadd(PDstandard4th33gyz,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmul(J23L,kmul(J33L,PDstandard4th23gyz))),ToReal(2)))))))); JacPDstandard4th12gxx = kmadd(J12L,kmadd(J11L,PDstandard4th11gxx,kmadd(J21L,PDstandard4th12gxx,kmul(J31L,PDstandard4th13gxx))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmadd(dJ112L,PDstandard4th1gxx,kmadd(J22L,kmadd(J21L,PDstandard4th22gxx,kmul(J31L,PDstandard4th23gxx)),kmadd(dJ212L,PDstandard4th2gxx,kmadd(J32L,kmadd(J21L,PDstandard4th23gxx,kmul(J31L,PDstandard4th33gxx)),kmul(dJ312L,PDstandard4th3gxx))))))); @@ -1113,11 +1117,12 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, } CCTK_REAL_VEC detg = - knmsub(gyyL,SQR(gxzL),knmsub(gxxL,SQR(gyzL),kmadd(gzzL,kmsub(gxxL,gyyL,SQR(gxyL)),kmul(gxyL,kmul(gxzL,kmul(gyzL,ToReal(2))))))); + knmsub(gyyL,kmul(gxzL,gxzL),knmsub(gxxL,kmul(gyzL,gyzL),kmadd(gzzL,kmsub(gxxL,gyyL,kmul(gxyL,gxyL)),kmul(gxyL,kmul(gxzL,kmul(gyzL,ToReal(2))))))); - CCTK_REAL_VEC invdetg = INV(detg); + CCTK_REAL_VEC invdetg = kdiv(ToReal(1),detg); - CCTK_REAL_VEC gInv11 = kmul(invdetg,kmsub(gyyL,gzzL,SQR(gyzL))); + CCTK_REAL_VEC gInv11 = + kmul(invdetg,kmsub(gyyL,gzzL,kmul(gyzL,gyzL))); CCTK_REAL_VEC gInv12 = kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); @@ -1128,7 +1133,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC gInv21 = kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); - CCTK_REAL_VEC gInv22 = kmul(invdetg,kmsub(gxxL,gzzL,SQR(gxzL))); + CCTK_REAL_VEC gInv22 = + kmul(invdetg,kmsub(gxxL,gzzL,kmul(gxzL,gxzL))); CCTK_REAL_VEC gInv23 = kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); @@ -1139,7 +1145,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC gInv32 = kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); - CCTK_REAL_VEC gInv33 = kmul(invdetg,kmsub(gxxL,gyyL,SQR(gxyL))); + CCTK_REAL_VEC gInv33 = + kmul(invdetg,kmsub(gxxL,gyyL,kmul(gxyL,gxyL))); CCTK_REAL_VEC gamma111 = kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard4th1gxx,knmsub(gInv12,JacPDstandard4th2gxx,kmsub(kmadd(gInv12,JacPDstandard4th1gxy,kmul(gInv13,JacPDstandard4th1gxz)),ToReal(2),kmul(gInv13,JacPDstandard4th3gxx))))); @@ -1229,13 +1236,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC wa3 = va3; CCTK_REAL_VEC omega11 = - kmadd(gxxL,SQR(wa1),kmadd(gyyL,SQR(wa2),kmadd(gzzL,SQR(wa3),kmul(kmadd(gyzL,kmul(wa2,wa3),kmul(wa1,kmadd(gxyL,wa2,kmul(gxzL,wa3)))),ToReal(2))))); + kmadd(gxxL,kmul(wa1,wa1),kmadd(gyyL,kmul(wa2,wa2),kmadd(gzzL,kmul(wa3,wa3),kmul(kmadd(gyzL,kmul(wa2,wa3),kmul(wa1,kmadd(gxyL,wa2,kmul(gxzL,wa3)))),ToReal(2))))); - CCTK_REAL_VEC ea1 = kmul(wa1,INV(ksqrt(omega11))); + CCTK_REAL_VEC ea1 = kdiv(wa1,ksqrt(omega11)); - CCTK_REAL_VEC ea2 = kmul(wa2,INV(ksqrt(omega11))); + CCTK_REAL_VEC ea2 = kdiv(wa2,ksqrt(omega11)); - CCTK_REAL_VEC ea3 = kmul(wa3,INV(ksqrt(omega11))); + CCTK_REAL_VEC ea3 = kdiv(wa3,ksqrt(omega11)); CCTK_REAL_VEC omega12 = kmadd(ea1,kmadd(gxxL,vb1,kmadd(gxyL,vb2,kmul(gxzL,vb3))),kmadd(ea2,kmadd(gxyL,vb1,kmadd(gyyL,vb2,kmul(gyzL,vb3))),kmul(ea3,kmadd(gxzL,vb1,kmadd(gyzL,vb2,kmul(gzzL,vb3)))))); @@ -1247,13 +1254,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC wb3 = knmsub(ea3,omega12,vb3); CCTK_REAL_VEC omega22 = - kmadd(gxxL,SQR(wb1),kmadd(gyyL,SQR(wb2),kmadd(gzzL,SQR(wb3),kmul(kmadd(gyzL,kmul(wb2,wb3),kmul(wb1,kmadd(gxyL,wb2,kmul(gxzL,wb3)))),ToReal(2))))); + kmadd(gxxL,kmul(wb1,wb1),kmadd(gyyL,kmul(wb2,wb2),kmadd(gzzL,kmul(wb3,wb3),kmul(kmadd(gyzL,kmul(wb2,wb3),kmul(wb1,kmadd(gxyL,wb2,kmul(gxzL,wb3)))),ToReal(2))))); - CCTK_REAL_VEC eb1 = kmul(wb1,INV(ksqrt(omega22))); + CCTK_REAL_VEC eb1 = kdiv(wb1,ksqrt(omega22)); - CCTK_REAL_VEC eb2 = kmul(wb2,INV(ksqrt(omega22))); + CCTK_REAL_VEC eb2 = kdiv(wb2,ksqrt(omega22)); - CCTK_REAL_VEC eb3 = kmul(wb3,INV(ksqrt(omega22))); + CCTK_REAL_VEC eb3 = kdiv(wb3,ksqrt(omega22)); CCTK_REAL_VEC omega13 = kmadd(ea1,kmadd(gxxL,vc1,kmadd(gxyL,vc2,kmul(gxzL,vc3))),kmadd(ea2,kmadd(gxyL,vc1,kmadd(gyyL,vc2,kmul(gyzL,vc3))),kmul(ea3,kmadd(gxzL,vc1,kmadd(gyzL,vc2,kmul(gzzL,vc3)))))); @@ -1268,13 +1275,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC wc3 = ksub(vc3,kmadd(eb3,omega23,kmul(ea3,omega13))); CCTK_REAL_VEC omega33 = - kmadd(gxxL,SQR(wc1),kmadd(gyyL,SQR(wc2),kmadd(gzzL,SQR(wc3),kmul(kmadd(gyzL,kmul(wc2,wc3),kmul(wc1,kmadd(gxyL,wc2,kmul(gxzL,wc3)))),ToReal(2))))); + kmadd(gxxL,kmul(wc1,wc1),kmadd(gyyL,kmul(wc2,wc2),kmadd(gzzL,kmul(wc3,wc3),kmul(kmadd(gyzL,kmul(wc2,wc3),kmul(wc1,kmadd(gxyL,wc2,kmul(gxzL,wc3)))),ToReal(2))))); - CCTK_REAL_VEC ec1 = kmul(wc1,INV(ksqrt(omega33))); + CCTK_REAL_VEC ec1 = kdiv(wc1,ksqrt(omega33)); - CCTK_REAL_VEC ec2 = kmul(wc2,INV(ksqrt(omega33))); + CCTK_REAL_VEC ec2 = kdiv(wc2,ksqrt(omega33)); - CCTK_REAL_VEC ec3 = kmul(wc3,INV(ksqrt(omega33))); + CCTK_REAL_VEC ec3 = kdiv(wc3,ksqrt(omega33)); CCTK_REAL_VEC isqrt2 = ToReal(0.707106781186547524); @@ -1340,13 +1347,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmul(ksub(JacPDstandard4th12gxx,JacPDstandard4th21gxx),ToReal(0.5)); CCTK_REAL_VEC R1212 = - kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma121),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma221),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma121,gamma121),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); CCTK_REAL_VEC R1213 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th12gxz,kadd(JacPDstandard4th31gxy,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2),JacPDstandard4th23gxx),JacPDstandard4th11gyz)))))); CCTK_REAL_VEC R1221 = - kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma221),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,SQR(gamma121))))); + kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,kmul(gamma121,gamma121))))); CCTK_REAL_VEC R1222 = kmul(ksub(JacPDstandard4th21gyy,JacPDstandard4th12gyy),ToReal(0.5)); @@ -1370,7 +1377,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gxy,kadd(JacPDstandard4th21gxz,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2),JacPDstandard4th32gxx),JacPDstandard4th11gyz)))))); CCTK_REAL_VEC R1313 = - kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma131),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma231),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma131,gamma131),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); CCTK_REAL_VEC R1321 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th23gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2)),kmul(gamma211,kmul(gamma232,ToReal(2)))))),JacPDstandard4th31gxy),JacPDstandard4th12gxz)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); @@ -1382,7 +1389,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th23gxz,kadd(JacPDstandard4th31gyz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2),JacPDstandard4th33gxy),JacPDstandard4th12gzz)))))); CCTK_REAL_VEC R1331 = - kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma231),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,SQR(gamma131))))); + kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,kmul(gamma131,gamma131))))); CCTK_REAL_VEC R1332 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2)),kmul(gamma221,kmul(gamma233,ToReal(2)))))),JacPDstandard4th32gxz),JacPDstandard4th13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); @@ -1394,13 +1401,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmul(ksub(JacPDstandard4th12gxx,JacPDstandard4th21gxx),ToReal(0.5)); CCTK_REAL_VEC R2112 = - kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma221),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,SQR(gamma121))))); + kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,kmul(gamma121,gamma121))))); CCTK_REAL_VEC R2113 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th32gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2)),kmul(gamma211,kmul(gamma232,ToReal(2)))))),JacPDstandard4th21gxz),JacPDstandard4th13gxy)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R2121 = - kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma121),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma221),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma121,gamma121),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); CCTK_REAL_VEC R2122 = kmul(ksub(JacPDstandard4th21gyy,JacPDstandard4th12gyy),ToReal(0.5)); @@ -1459,13 +1466,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmul(ksub(JacPDstandard4th23gyy,JacPDstandard4th32gyy),ToReal(0.5)); CCTK_REAL_VEC R2323 = - kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma132),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma232),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma132,gamma132),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); CCTK_REAL_VEC R2331 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th12gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2)),kmul(gamma221,kmul(gamma233,ToReal(2)))))),JacPDstandard4th31gyz),JacPDstandard4th23gxz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); CCTK_REAL_VEC R2332 = - kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma232),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,SQR(gamma132))))); + kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,kmul(gamma132,gamma132))))); CCTK_REAL_VEC R2333 = kmul(ksub(JacPDstandard4th32gzz,JacPDstandard4th23gzz),ToReal(0.5)); @@ -1477,7 +1484,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th23gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2)),kmul(gamma211,kmul(gamma232,ToReal(2)))))),JacPDstandard4th31gxy),JacPDstandard4th12gxz)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R3113 = - kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma231),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,SQR(gamma131))))); + kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,kmul(gamma131,gamma131))))); CCTK_REAL_VEC R3121 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gxy,kadd(JacPDstandard4th21gxz,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2),JacPDstandard4th32gxx),JacPDstandard4th11gyz)))))); @@ -1489,7 +1496,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2)),kmul(gamma221,kmul(gamma233,ToReal(2)))))),JacPDstandard4th32gxz),JacPDstandard4th13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); CCTK_REAL_VEC R3131 = - kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma131),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma231),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma131,gamma131),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); CCTK_REAL_VEC R3132 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th23gxz,kadd(JacPDstandard4th31gyz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2),JacPDstandard4th33gxy),JacPDstandard4th12gzz)))))); @@ -1513,13 +1520,13 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, kmul(ksub(JacPDstandard4th23gyy,JacPDstandard4th32gyy),ToReal(0.5)); CCTK_REAL_VEC R3223 = - kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4))),kmadd(gyyL,kmul(SQR(gamma232),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,SQR(gamma132))))); + kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4))),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(-2)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(-2)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,kmul(gamma132,gamma132))))); CCTK_REAL_VEC R3231 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gyz,kadd(JacPDstandard4th32gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2),JacPDstandard4th33gxy),JacPDstandard4th21gzz)))))); CCTK_REAL_VEC R3232 = - kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(SQR(gamma132),ToReal(2)),kmadd(gyyL,kmul(SQR(gamma232),ToReal(2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma132,gamma132),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); CCTK_REAL_VEC R3233 = kmul(ksub(JacPDstandard4th32gzz,JacPDstandard4th23gzz),ToReal(0.5)); @@ -1570,12 +1577,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p1211 = R1211; - CCTK_REAL_VEC R4p1212 = kmadd(kxxL,kyyL,ksub(R1212,SQR(kxyL))); + CCTK_REAL_VEC R4p1212 = + kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)); CCTK_REAL_VEC R4p1213 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)); - CCTK_REAL_VEC R4p1221 = kadd(R1221,knmsub(kxxL,kyyL,SQR(kxyL))); + CCTK_REAL_VEC R4p1221 = + kadd(R1221,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))); CCTK_REAL_VEC R4p1222 = R1222; @@ -1595,7 +1604,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p1312 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1312)); - CCTK_REAL_VEC R4p1313 = kmadd(kxxL,kzzL,ksub(R1313,SQR(kxzL))); + CCTK_REAL_VEC R4p1313 = + kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)); CCTK_REAL_VEC R4p1321 = kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1321)); @@ -1605,7 +1615,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p1323 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)); - CCTK_REAL_VEC R4p1331 = kadd(R1331,knmsub(kxxL,kzzL,SQR(kxzL))); + CCTK_REAL_VEC R4p1331 = + kadd(R1331,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))); CCTK_REAL_VEC R4p1332 = kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1332)); @@ -1614,12 +1625,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p2111 = R2111; - CCTK_REAL_VEC R4p2112 = kadd(R2112,knmsub(kxxL,kyyL,SQR(kxyL))); + CCTK_REAL_VEC R4p2112 = + kadd(R2112,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))); CCTK_REAL_VEC R4p2113 = kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R2113)); - CCTK_REAL_VEC R4p2121 = kmadd(kxxL,kyyL,ksub(R2121,SQR(kxyL))); + CCTK_REAL_VEC R4p2121 = + kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R2121)); CCTK_REAL_VEC R4p2122 = R2122; @@ -1665,12 +1678,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p2322 = R2322; - CCTK_REAL_VEC R4p2323 = kmadd(kyyL,kzzL,ksub(R2323,SQR(kyzL))); + CCTK_REAL_VEC R4p2323 = + kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323)); CCTK_REAL_VEC R4p2331 = kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R2331)); - CCTK_REAL_VEC R4p2332 = kadd(R2332,knmsub(kyyL,kzzL,SQR(kyzL))); + CCTK_REAL_VEC R4p2332 = + kadd(R2332,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))); CCTK_REAL_VEC R4p2333 = R2333; @@ -1679,7 +1694,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p3112 = kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R3112)); - CCTK_REAL_VEC R4p3113 = kadd(R3113,knmsub(kxxL,kzzL,SQR(kxzL))); + CCTK_REAL_VEC R4p3113 = + kadd(R3113,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))); CCTK_REAL_VEC R4p3121 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)); @@ -1689,7 +1705,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p3123 = kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3123)); - CCTK_REAL_VEC R4p3131 = kmadd(kxxL,kzzL,ksub(R3131,SQR(kxzL))); + CCTK_REAL_VEC R4p3131 = + kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R3131)); CCTK_REAL_VEC R4p3132 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)); @@ -1709,12 +1726,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC R4p3222 = R3222; - CCTK_REAL_VEC R4p3223 = kadd(R3223,knmsub(kyyL,kzzL,SQR(kyzL))); + CCTK_REAL_VEC R4p3223 = + kadd(R3223,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))); CCTK_REAL_VEC R4p3231 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3231)); - CCTK_REAL_VEC R4p3232 = kmadd(kyyL,kzzL,ksub(R3232,SQR(kyzL))); + CCTK_REAL_VEC R4p3232 = + kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R3232)); CCTK_REAL_VEC R4p3233 = R3233; @@ -1809,37 +1828,37 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC Ro333 = ToReal(0); CCTK_REAL_VEC Rojo11 = - kmadd(gInv11,R1111,kmadd(gInv12,R1112,kmadd(gInv13,R1113,kmadd(gInv21,R1211,kmadd(gInv31,R1311,kmadd(gInv23,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv32,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1312)),kmadd(gInv22,kmadd(kxxL,kyyL,ksub(R1212,SQR(kxyL))),kmul(gInv33,kmadd(kxxL,kzzL,ksub(R1313,SQR(kxzL)))))))))))); + kmadd(gInv11,R1111,kmadd(gInv12,R1112,kmadd(gInv13,R1113,kmadd(gInv21,R1211,kmadd(gInv31,R1311,kmadd(gInv23,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv32,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1312)),kmadd(gInv22,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313))))))))))); CCTK_REAL_VEC Rojo12 = - kmadd(gInv11,R1121,kmadd(gInv21,R1221,kmadd(gInv22,R1222,kmadd(gInv23,R1223,kmadd(gInv31,R1321,kmadd(gInv12,kadd(R1122,knmsub(kxxL,kyyL,SQR(kxyL))),kmadd(gInv32,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1322)),kmadd(gInv13,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1123)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))))); + kmadd(gInv11,R1121,kmadd(gInv21,R1221,kmadd(gInv22,R1222,kmadd(gInv23,R1223,kmadd(gInv31,R1321,kmadd(gInv12,kadd(R1122,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))),kmadd(gInv32,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1322)),kmadd(gInv13,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1123)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))))); CCTK_REAL_VEC Rojo13 = - kmadd(gInv11,R1131,kmadd(gInv21,R1231,kmadd(gInv31,R1331,kmadd(gInv32,R1332,kmadd(gInv33,R1333,kmadd(gInv12,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1132)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R1232)),kmadd(gInv13,kadd(R1133,knmsub(kxxL,kzzL,SQR(kxzL))),kmul(gInv23,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1233))))))))))); + kmadd(gInv11,R1131,kmadd(gInv21,R1231,kmadd(gInv31,R1331,kmadd(gInv32,R1332,kmadd(gInv33,R1333,kmadd(gInv12,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1132)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R1232)),kmadd(gInv13,kadd(R1133,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))),kmul(gInv23,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1233))))))))))); CCTK_REAL_VEC Rojo21 = - kmadd(gInv11,R2111,kmadd(gInv12,R2112,kmadd(gInv13,R2113,kmadd(gInv22,R2212,kmadd(gInv32,R2312,kmadd(gInv21,kadd(R2211,knmsub(kxxL,kyyL,SQR(kxyL))),kmadd(gInv23,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2213)),kmadd(gInv31,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R2311)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R2313))))))))))); + kmadd(gInv11,R2111,kmadd(gInv12,R2112,kmadd(gInv13,R2113,kmadd(gInv22,R2212,kmadd(gInv32,R2312,kmadd(gInv21,kadd(R2211,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))),kmadd(gInv23,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2213)),kmadd(gInv31,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R2311)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R2313))))))))))); CCTK_REAL_VEC Rojo22 = - kmadd(gInv12,R2122,kmadd(gInv21,R2221,kmadd(gInv22,R2222,kmadd(gInv23,R2223,kmadd(gInv32,R2322,kmadd(gInv13,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2123)),kmadd(gInv31,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2321)),kmadd(gInv11,kmadd(kxxL,kyyL,ksub(R2121,SQR(kxyL))),kmul(gInv33,kmadd(kyyL,kzzL,ksub(R2323,SQR(kyzL)))))))))))); + kmadd(gInv12,R2122,kmadd(gInv21,R2221,kmadd(gInv22,R2222,kmadd(gInv23,R2223,kmadd(gInv32,R2322,kmadd(gInv13,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2123)),kmadd(gInv31,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2321)),kmadd(gInv11,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R2121)),kmul(gInv33,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))))))))); CCTK_REAL_VEC Rojo23 = - kmadd(gInv12,R2132,kmadd(gInv22,R2232,kmadd(gInv31,R2331,kmadd(gInv32,R2332,kmadd(gInv33,R2333,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R2131)),kmadd(gInv21,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2231)),kmadd(gInv13,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R2133)),kmul(gInv23,kadd(R2233,knmsub(kyyL,kzzL,SQR(kyzL)))))))))))); + kmadd(gInv12,R2132,kmadd(gInv22,R2232,kmadd(gInv31,R2331,kmadd(gInv32,R2332,kmadd(gInv33,R2333,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R2131)),kmadd(gInv21,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2231)),kmadd(gInv13,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R2133)),kmul(gInv23,kadd(R2233,kmsub(kyzL,kyzL,kmul(kyyL,kzzL)))))))))))); CCTK_REAL_VEC Rojo31 = - kmadd(gInv11,R3111,kmadd(gInv12,R3112,kmadd(gInv13,R3113,kmadd(gInv23,R3213,kmadd(gInv33,R3313,kmadd(gInv21,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R3211)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R3212)),kmadd(gInv31,kadd(R3311,knmsub(kxxL,kzzL,SQR(kxzL))),kmul(gInv32,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3312))))))))))); + kmadd(gInv11,R3111,kmadd(gInv12,R3112,kmadd(gInv13,R3113,kmadd(gInv23,R3213,kmadd(gInv33,R3313,kmadd(gInv21,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R3211)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R3212)),kmadd(gInv31,kadd(R3311,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))),kmul(gInv32,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3312))))))))))); CCTK_REAL_VEC Rojo32 = - kmadd(gInv13,R3123,kmadd(gInv21,R3221,kmadd(gInv22,R3222,kmadd(gInv23,R3223,kmadd(gInv33,R3323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)),kmadd(gInv12,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R3122)),kmadd(gInv31,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3321)),kmul(gInv32,kadd(R3322,knmsub(kyyL,kzzL,SQR(kyzL)))))))))))); + kmadd(gInv13,R3123,kmadd(gInv21,R3221,kmadd(gInv22,R3222,kmadd(gInv23,R3223,kmadd(gInv33,R3323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)),kmadd(gInv12,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R3122)),kmadd(gInv31,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3321)),kmul(gInv32,kadd(R3322,kmsub(kyzL,kyzL,kmul(kyyL,kzzL)))))))))))); CCTK_REAL_VEC Rojo33 = - kmadd(gInv13,R3133,kmadd(gInv23,R3233,kmadd(gInv31,R3331,kmadd(gInv32,R3332,kmadd(gInv33,R3333,kmadd(gInv12,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)),kmadd(gInv21,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3231)),kmadd(gInv11,kmadd(kxxL,kzzL,ksub(R3131,SQR(kxzL))),kmul(gInv22,kmadd(kyyL,kzzL,ksub(R3232,SQR(kyzL)))))))))))); + kmadd(gInv13,R3133,kmadd(gInv23,R3233,kmadd(gInv31,R3331,kmadd(gInv32,R3332,kmadd(gInv33,R3333,kmadd(gInv12,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)),kmadd(gInv21,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3231)),kmadd(gInv11,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R3131)),kmul(gInv22,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R3232))))))))))); CCTK_REAL_VEC Psi4rL = - kmadd(R4p1112,kmul(rmbar1,kmul(rmbar2,SQR(n1))),kmadd(R4p1211,kmul(rmbar1,kmul(rmbar2,SQR(n1))),kmadd(R4p1113,kmul(rmbar1,kmul(rmbar3,SQR(n1))),kmadd(R4p1311,kmul(rmbar1,kmul(rmbar3,SQR(n1))),kmadd(R4p1213,kmul(rmbar2,kmul(rmbar3,SQR(n1))),kmadd(R4p1312,kmul(rmbar2,kmul(rmbar3,SQR(n1))),kmadd(R4p2122,kmul(rmbar1,kmul(rmbar2,SQR(n2))),kmadd(R4p2221,kmul(rmbar1,kmul(rmbar2,SQR(n2))),kmadd(R4p2123,kmul(rmbar1,kmul(rmbar3,SQR(n2))),kmadd(R4p2321,kmul(rmbar1,kmul(rmbar3,SQR(n2))),kmadd(R4p2223,kmul(rmbar2,kmul(rmbar3,SQR(n2))),kmadd(R4p2322,kmul(rmbar2,kmul(rmbar3,SQR(n2))),kmadd(R4p3132,kmul(rmbar1,kmul(rmbar2,SQR(n3))),kmadd(R4p3231,kmul(rmbar1,kmul(rmbar2,SQR(n3))),kmadd(R4p3133,kmul(rmbar1,kmul(rmbar3,SQR(n3))),kmadd(R4p3331,kmul(rmbar1,kmul(rmbar3,SQR(n3))),kmadd(R4p3233,kmul(rmbar2,kmul(rmbar3,SQR(n3))),kmadd(R4p3332,kmul(rmbar2,kmul(rmbar3,SQR(n3))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo12,SQR(nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo13,SQR(nn))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo21,SQR(nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo23,SQR(nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo31,SQR(nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo32,SQR(nn))),kmadd(n1,kmul(n2,kmul(R4p1121,SQR(rmbar1))),kmadd(n1,kmul(n3,kmul(R4p1131,SQR(rmbar1))),kmadd(n1,kmul(n2,kmul(R4p2111,SQR(rmbar1))),kmadd(n2,kmul(n3,kmul(R4p2131,SQR(rmbar1))),kmadd(n1,kmul(n3,kmul(R4p3111,SQR(rmbar1))),kmadd(n2,kmul(n3,kmul(R4p3121,SQR(rmbar1))),kmadd(R4p1111,SQR(kmul(n1,rmbar1)),kmadd(R4p2121,SQR(kmul(n2,rmbar1)),kmadd(R4p3131,SQR(kmul(n3,rmbar1)),kmadd(Rojo11,SQR(kmul(nn,rmbar1)),kmadd(n1,kmul(n2,kmul(R4p1222,SQR(rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1232,SQR(rmbar2))),kmadd(n1,kmul(n2,kmul(R4p2212,SQR(rmbar2))),kmadd(n2,kmul(n3,kmul(R4p2232,SQR(rmbar2))),kmadd(n1,kmul(n3,kmul(R4p3212,SQR(rmbar2))),kmadd(n2,kmul(n3,kmul(R4p3222,SQR(rmbar2))),kmadd(R4p1212,SQR(kmul(n1,rmbar2)),kmadd(R4p2222,SQR(kmul(n2,rmbar2)),kmadd(R4p3232,SQR(kmul(n3,rmbar2)),kmadd(Rojo22,SQR(kmul(nn,rmbar2)),kmadd(n1,kmul(n2,kmul(R4p1323,SQR(rmbar3))),kmadd(n1,kmul(n3,kmul(R4p1333,SQR(rmbar3))),kmadd(n1,kmul(n2,kmul(R4p2313,SQR(rmbar3))),kmadd(n2,kmul(n3,kmul(R4p2333,SQR(rmbar3))),kmadd(n1,kmul(n3,kmul(R4p3313,SQR(rmbar3))),kmadd(n2,kmul(n3,kmul(R4p3323,SQR(rmbar3))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(SQR(imbar3),ToReal(-2)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(SQR(imbar3),ToReal(-2)))),knmsub(n3,kmul(SQR(imbar3),kmadd(n1,R4p3313,kmul(n2,R4p2333))),kmadd(SQR(imbar3),kmsub(n1,kmul(nn,kmul(Ro313,ToReal(-2))),kmul(n2,kmul(n3,R4p3323))),kmadd(R4p1313,ksub(SQR(kmul(n1,rmbar3)),SQR(kmul(imbar3,n1))),kmadd(R4p2323,ksub(SQR(kmul(n2,rmbar3)),SQR(kmul(imbar3,n2))),kmadd(R4p3333,ksub(SQR(kmul(n3,rmbar3)),SQR(kmul(imbar3,n3))),kmadd(Rojo33,ksub(SQR(kmul(nn,rmbar3)),SQR(kmul(imbar3,nn))),kmadd(n1,kmul(nn,kmul(Ro111,kmul(SQR(rmbar1),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro121,kmul(SQR(rmbar1),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro131,kmul(SQR(rmbar1),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro212,kmul(SQR(rmbar2),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro222,kmul(SQR(rmbar2),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro232,kmul(SQR(rmbar2),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(SQR(rmbar3),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(SQR(rmbar3),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(SQR(rmbar3),ToReal(2)))),kmadd(n1,kmadd(n3,kmul(R4p1132,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1221,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p1231,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1123,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1331,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2113,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1233,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p1322,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p1332,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p2213,kmul(rmbar2,rmbar3)),kmadd(n2,kmsub(R4p1122,kmul(rmbar1,rmbar2),kmul(R4p1323,SQR(imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro112,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro113,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro211,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro213,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro311,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro312,ToReal(2))))))))))))))))))))),kmadd(n2,kmadd(n3,kmul(R4p2132,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p2211,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p2231,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1321,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p2133,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p2311,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p3321,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p1223,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p2233,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p2312,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p2112,kmul(rmbar1,rmbar2),kmul(R4p2313,SQR(imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro122,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro123,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro221,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro223,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro321,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro322,ToReal(2))))))))))))))))))))),kmsub(n3,kmadd(n2,kmul(R4p3122,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p3211,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p3221,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1133,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2331,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3113,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p3123,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3311,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2332,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3213,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3223,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3312,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3322,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p3112,kmul(rmbar1,rmbar2),kmul(R4p1333,SQR(imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro132,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro133,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro231,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro233,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro331,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro332,ToReal(2)))))))))))))))))))))))),kmadd(SQR(imbar1),kmadd(n2,kmul(n3,R4p2131),kmadd(n2,kmul(n3,R4p3121),kmadd(R4p1111,SQR(n1),kmadd(R4p2121,SQR(n2),kmadd(R4p3131,SQR(n3),kmadd(Rojo11,SQR(nn),kmadd(n2,kmul(nn,kmul(Ro121,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro131,ToReal(2))),kmul(n1,kmadd(n2,kadd(R4p1121,R4p2111),kmadd(n3,kadd(R4p1131,R4p3111),kmul(nn,kmul(Ro111,ToReal(2)))))))))))))),kmadd(SQR(imbar2),kmadd(n2,kmul(n3,R4p2232),kmadd(n2,kmul(n3,R4p3222),kmadd(R4p1212,SQR(n1),kmadd(R4p2222,SQR(n2),kmadd(R4p3232,SQR(n3),kmadd(Rojo22,SQR(nn),kmadd(n2,kmul(nn,kmul(Ro222,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro232,ToReal(2))),kmul(n1,kmadd(n2,kadd(R4p1222,R4p2212),kmadd(n3,kadd(R4p1232,R4p3212),kmul(nn,kmul(Ro212,ToReal(2)))))))))))))),kmadd(imbar2,kmul(imbar3,kmadd(kadd(R4p1213,R4p1312),SQR(n1),kmadd(kadd(R4p2223,R4p2322),SQR(n2),kmadd(R4p3233,SQR(n3),kmadd(R4p3332,SQR(n3),kmadd(Rojo23,SQR(nn),kmadd(Rojo32,SQR(nn),kmadd(n3,kmul(nn,kmul(Ro233,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro332,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1223,kadd(R4p1322,kadd(R4p2213,R4p2312))),kmadd(n3,kadd(R4p1233,kadd(R4p1332,kadd(R4p3213,R4p3312))),kmul(nn,kmul(kadd(Ro213,Ro312),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2233,kadd(R4p2332,kadd(R4p3223,R4p3322))),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2))))))))))))))),kmul(imbar1,kmadd(imbar2,kmadd(kadd(R4p1112,R4p1211),SQR(n1),kmadd(kadd(R4p2122,R4p2221),SQR(n2),kmadd(R4p3132,SQR(n3),kmadd(R4p3231,SQR(n3),kmadd(Rojo12,SQR(nn),kmadd(Rojo21,SQR(nn),kmadd(n3,kmul(nn,kmul(Ro132,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro231,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1122,kadd(R4p1221,kadd(R4p2112,R4p2211))),kmadd(n3,kadd(R4p1132,kadd(R4p1231,kadd(R4p3112,R4p3211))),kmul(nn,kmul(kadd(Ro112,Ro211),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2132,kadd(R4p2231,kadd(R4p3122,R4p3221))),kmul(nn,kmul(kadd(Ro122,Ro221),ToReal(2)))))))))))))),kmul(imbar3,kmadd(kadd(R4p1113,R4p1311),SQR(n1),kmadd(kadd(R4p2123,R4p2321),SQR(n2),kmadd(R4p3133,SQR(n3),kmadd(R4p3331,SQR(n3),kmadd(Rojo13,SQR(nn),kmadd(Rojo31,SQR(nn),kmadd(n3,kmul(nn,kmul(Ro133,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro331,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1123,kadd(R4p1321,kadd(R4p2113,R4p2311))),kmadd(n3,kadd(R4p1133,kadd(R4p1331,kadd(R4p3113,R4p3311))),kmul(nn,kmul(kadd(Ro113,Ro311),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2133,kadd(R4p2331,kadd(R4p3123,R4p3321))),kmul(nn,kmul(kadd(Ro123,Ro321),ToReal(2)))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))); + kmadd(R4p1112,kmul(rmbar1,kmul(rmbar2,kmul(n1,n1))),kmadd(R4p1211,kmul(rmbar1,kmul(rmbar2,kmul(n1,n1))),kmadd(R4p1113,kmul(rmbar1,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1311,kmul(rmbar1,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1213,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1312,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p2122,kmul(rmbar1,kmul(rmbar2,kmul(n2,n2))),kmadd(R4p2221,kmul(rmbar1,kmul(rmbar2,kmul(n2,n2))),kmadd(R4p2123,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2321,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2223,kmul(rmbar2,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2322,kmul(rmbar2,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p3132,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3))),kmadd(R4p3231,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3))),kmadd(R4p3133,kmul(rmbar1,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3331,kmul(rmbar1,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3233,kmul(rmbar2,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3332,kmul(rmbar2,kmul(rmbar3,kmul(n3,n3))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo12,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo13,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo21,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo23,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo31,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo32,kmul(nn,nn))),kmadd(n1,kmul(n2,kmul(R4p1121,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n3,kmul(R4p1131,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n2,kmul(R4p2111,kmul(rmbar1,rmbar1))),kmadd(n2,kmul(n3,kmul(R4p2131,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n3,kmul(R4p3111,kmul(rmbar1,rmbar1))),kmadd(n2,kmul(n3,kmul(R4p3121,kmul(rmbar1,rmbar1))),kmadd(R4p1111,kmul(kmul(n1,n1),kmul(rmbar1,rmbar1)),kmadd(R4p2121,kmul(kmul(n2,n2),kmul(rmbar1,rmbar1)),kmadd(R4p3131,kmul(kmul(n3,n3),kmul(rmbar1,rmbar1)),kmadd(Rojo11,kmul(kmul(nn,nn),kmul(rmbar1,rmbar1)),kmadd(n1,kmul(n2,kmul(R4p1222,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1232,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n2,kmul(R4p2212,kmul(rmbar2,rmbar2))),kmadd(n2,kmul(n3,kmul(R4p2232,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p3212,kmul(rmbar2,rmbar2))),kmadd(n2,kmul(n3,kmul(R4p3222,kmul(rmbar2,rmbar2))),kmadd(R4p1212,kmul(kmul(n1,n1),kmul(rmbar2,rmbar2)),kmadd(R4p2222,kmul(kmul(n2,n2),kmul(rmbar2,rmbar2)),kmadd(R4p3232,kmul(kmul(n3,n3),kmul(rmbar2,rmbar2)),kmadd(Rojo22,kmul(kmul(nn,nn),kmul(rmbar2,rmbar2)),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p1333,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n2,kmul(R4p2313,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p2333,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p3313,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p3323,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(imbar3,imbar3),ToReal(-2)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(imbar3,imbar3),ToReal(-2)))),knmsub(n3,kmul(kmul(imbar3,imbar3),kmadd(n1,R4p3313,kmul(n2,R4p2333))),kmadd(kmul(imbar3,imbar3),kmsub(n1,kmul(nn,kmul(Ro313,ToReal(-2))),kmul(n2,kmul(n3,R4p3323))),kmadd(R4p1313,kmsub(kmul(n1,n1),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n1,n1))),kmadd(R4p2323,kmsub(kmul(n2,n2),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n2,n2))),kmadd(R4p3333,kmsub(kmul(n3,n3),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n3,n3))),kmadd(Rojo33,kmsub(kmul(nn,nn),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(nn,nn))),kmadd(n1,kmul(nn,kmul(Ro111,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro121,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro131,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro212,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro222,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro232,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n1,kmadd(n3,kmul(R4p1132,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1221,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p1231,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1123,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1331,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2113,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1233,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p1322,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p1332,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p2213,kmul(rmbar2,rmbar3)),kmadd(n2,kmsub(R4p1122,kmul(rmbar1,rmbar2),kmul(R4p1323,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro112,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro113,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro211,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro213,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro311,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro312,ToReal(2))))))))))))))))))))),kmadd(n2,kmadd(n3,kmul(R4p2132,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p2211,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p2231,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1321,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p2133,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p2311,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p3321,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p1223,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p2233,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p2312,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p2112,kmul(rmbar1,rmbar2),kmul(R4p2313,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro122,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro123,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro221,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro223,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro321,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro322,ToReal(2))))))))))))))))))))),kmsub(n3,kmadd(n2,kmul(R4p3122,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p3211,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p3221,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1133,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2331,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3113,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p3123,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3311,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2332,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3213,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3223,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3312,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3322,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p3112,kmul(rmbar1,rmbar2),kmul(R4p1333,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro132,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro133,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro231,ToReal(2)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro233,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro331,ToReal(2)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro332,ToReal(2)))))))))))))))))))))))),kmadd(kmul(imbar1,imbar1),kmadd(n2,kmul(n3,R4p2131),kmadd(n2,kmul(n3,R4p3121),kmadd(R4p1111,kmul(n1,n1),kmadd(R4p2121,kmul(n2,n2),kmadd(R4p3131,kmul(n3,n3),kmadd(Rojo11,kmul(nn,nn),kmadd(n2,kmul(nn,kmul(Ro121,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro131,ToReal(2))),kmul(n1,kmadd(n2,kadd(R4p1121,R4p2111),kmadd(n3,kadd(R4p1131,R4p3111),kmul(nn,kmul(Ro111,ToReal(2)))))))))))))),kmadd(kmul(imbar2,imbar2),kmadd(n2,kmul(n3,R4p2232),kmadd(n2,kmul(n3,R4p3222),kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2222,kmul(n2,n2),kmadd(R4p3232,kmul(n3,n3),kmadd(Rojo22,kmul(nn,nn),kmadd(n2,kmul(nn,kmul(Ro222,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro232,ToReal(2))),kmul(n1,kmadd(n2,kadd(R4p1222,R4p2212),kmadd(n3,kadd(R4p1232,R4p3212),kmul(nn,kmul(Ro212,ToReal(2)))))))))))))),kmadd(imbar2,kmul(imbar3,kmadd(kadd(R4p1213,R4p1312),kmul(n1,n1),kmadd(kadd(R4p2223,R4p2322),kmul(n2,n2),kmadd(R4p3233,kmul(n3,n3),kmadd(R4p3332,kmul(n3,n3),kmadd(Rojo23,kmul(nn,nn),kmadd(Rojo32,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro233,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro332,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1223,kadd(R4p1322,kadd(R4p2213,R4p2312))),kmadd(n3,kadd(R4p1233,kadd(R4p1332,kadd(R4p3213,R4p3312))),kmul(nn,kmul(kadd(Ro213,Ro312),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2233,kadd(R4p2332,kadd(R4p3223,R4p3322))),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2))))))))))))))),kmul(imbar1,kmadd(imbar2,kmadd(kadd(R4p1112,R4p1211),kmul(n1,n1),kmadd(kadd(R4p2122,R4p2221),kmul(n2,n2),kmadd(R4p3132,kmul(n3,n3),kmadd(R4p3231,kmul(n3,n3),kmadd(Rojo12,kmul(nn,nn),kmadd(Rojo21,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro132,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro231,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1122,kadd(R4p1221,kadd(R4p2112,R4p2211))),kmadd(n3,kadd(R4p1132,kadd(R4p1231,kadd(R4p3112,R4p3211))),kmul(nn,kmul(kadd(Ro112,Ro211),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2132,kadd(R4p2231,kadd(R4p3122,R4p3221))),kmul(nn,kmul(kadd(Ro122,Ro221),ToReal(2)))))))))))))),kmul(imbar3,kmadd(kadd(R4p1113,R4p1311),kmul(n1,n1),kmadd(kadd(R4p2123,R4p2321),kmul(n2,n2),kmadd(R4p3133,kmul(n3,n3),kmadd(R4p3331,kmul(n3,n3),kmadd(Rojo13,kmul(nn,nn),kmadd(Rojo31,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro133,ToReal(2))),kmadd(n3,kmul(nn,kmul(Ro331,ToReal(2))),kmadd(n1,kmadd(n2,kadd(R4p1123,kadd(R4p1321,kadd(R4p2113,R4p2311))),kmadd(n3,kadd(R4p1133,kadd(R4p1331,kadd(R4p3113,R4p3311))),kmul(nn,kmul(kadd(Ro113,Ro311),ToReal(2))))),kmul(n2,kmadd(n3,kadd(R4p2133,kadd(R4p2331,kadd(R4p3123,R4p3321))),kmul(nn,kmul(kadd(Ro123,Ro321),ToReal(2)))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))); CCTK_REAL_VEC Psi4iL = - knmadd(im1,kmadd(R4p3132,kmul(rm2,SQR(n3)),kmadd(R4p3231,kmul(rm2,SQR(n3)),kmadd(R4p3133,kmul(rm3,SQR(n3)),kmadd(R4p3331,kmul(rm3,SQR(n3)),kmadd(rm2,kmul(Rojo12,SQR(nn)),kmadd(rm3,kmul(Rojo13,SQR(nn)),kmadd(rm2,kmul(Rojo21,SQR(nn)),kmadd(rm3,kmul(Rojo31,SQR(nn)),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro132,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro133,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro231,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro331,ToReal(2)))),kmadd(R4p3131,kmul(rm1,kmul(SQR(n3),ToReal(2))),kmadd(rm1,kmul(Rojo11,kmul(SQR(nn),ToReal(2))),kmadd(SQR(n1),kmadd(R4p1112,rm2,kmadd(R4p1211,rm2,kmadd(R4p1113,rm3,kmadd(R4p1311,rm3,kmul(R4p1111,kmul(rm1,ToReal(2))))))),kmadd(SQR(n2),kmadd(R4p2122,rm2,kmadd(R4p2221,rm2,kmadd(R4p2123,rm3,kmadd(R4p2321,rm3,kmul(R4p2121,kmul(rm1,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm2,kmadd(R4p1221,rm2,kmadd(R4p2112,rm2,kmadd(R4p2211,rm2,kmadd(R4p1123,rm3,kmadd(R4p1321,rm3,kmadd(R4p2113,rm3,kmadd(R4p2311,rm3,kmadd(R4p1121,kmul(rm1,ToReal(2)),kmul(R4p2111,kmul(rm1,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1132,rm2,kmadd(R4p1231,rm2,kmadd(R4p3112,rm2,kmadd(R4p3211,rm2,kmadd(R4p1133,rm3,kmadd(R4p1331,rm3,kmadd(R4p3113,rm3,kmadd(R4p3311,rm3,kmadd(R4p1131,kmul(rm1,ToReal(2)),kmul(R4p3111,kmul(rm1,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm2,kmadd(R4p2231,rm2,kmadd(R4p3122,rm2,kmadd(R4p3221,rm2,kmadd(R4p2133,rm3,kmadd(R4p2331,rm3,kmadd(R4p3123,rm3,kmadd(R4p3321,rm3,kmadd(R4p2131,kmul(rm1,ToReal(2)),kmul(R4p3121,kmul(rm1,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm1,kmul(Ro131,ToReal(4))))))))))))))))))))))),kmadd(im3,kmadd(R4p3133,kmul(rm1,SQR(n3)),kmadd(R4p3331,kmul(rm1,SQR(n3)),kmadd(R4p3233,kmul(rm2,SQR(n3)),kmadd(R4p3332,kmul(rm2,SQR(n3)),kmadd(rm1,kmul(Rojo13,SQR(nn)),kmadd(rm2,kmul(Rojo23,SQR(nn)),kmadd(rm1,kmul(Rojo31,SQR(nn)),kmadd(rm2,kmul(Rojo32,SQR(nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro133,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro233,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro331,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro332,ToReal(2)))),kmadd(R4p3333,kmul(rm3,kmul(SQR(n3),ToReal(2))),kmadd(rm3,kmul(Rojo33,kmul(SQR(nn),ToReal(2))),kmadd(SQR(n1),kmadd(R4p1113,rm1,kmadd(R4p1311,rm1,kmadd(R4p1213,rm2,kmadd(R4p1312,rm2,kmul(R4p1313,kmul(rm3,ToReal(2))))))),kmadd(SQR(n2),kmadd(R4p2123,rm1,kmadd(R4p2321,rm1,kmadd(R4p2223,rm2,kmadd(R4p2322,rm2,kmul(R4p2323,kmul(rm3,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1123,rm1,kmadd(R4p1321,rm1,kmadd(R4p2113,rm1,kmadd(R4p2311,rm1,kmadd(R4p1223,rm2,kmadd(R4p1322,rm2,kmadd(R4p2213,rm2,kmadd(R4p2312,rm2,kmadd(R4p1323,kmul(rm3,ToReal(2)),kmul(R4p2313,kmul(rm3,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1133,rm1,kmadd(R4p1331,rm1,kmadd(R4p3113,rm1,kmadd(R4p3311,rm1,kmadd(R4p1233,rm2,kmadd(R4p1332,rm2,kmadd(R4p3213,rm2,kmadd(R4p3312,rm2,kmadd(R4p1333,kmul(rm3,ToReal(2)),kmul(R4p3313,kmul(rm3,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2133,rm1,kmadd(R4p2331,rm1,kmadd(R4p3123,rm1,kmadd(R4p3321,rm1,kmadd(R4p2233,rm2,kmadd(R4p2332,rm2,kmadd(R4p3223,rm2,kmadd(R4p3322,rm2,kmadd(R4p2333,kmul(rm3,ToReal(2)),kmul(R4p3323,kmul(rm3,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm3,kmul(Ro333,ToReal(4))))))))))))))))))))))),kmul(im2,kmadd(R4p3132,kmul(rm1,SQR(n3)),kmadd(R4p3231,kmul(rm1,SQR(n3)),kmadd(R4p3233,kmul(rm3,SQR(n3)),kmadd(R4p3332,kmul(rm3,SQR(n3)),kmadd(rm1,kmul(Rojo12,SQR(nn)),kmadd(rm1,kmul(Rojo21,SQR(nn)),kmadd(rm3,kmul(Rojo23,SQR(nn)),kmadd(rm3,kmul(Rojo32,SQR(nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro132,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro231,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro233,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro332,ToReal(2)))),kmadd(R4p3232,kmul(rm2,kmul(SQR(n3),ToReal(2))),kmadd(rm2,kmul(Rojo22,kmul(SQR(nn),ToReal(2))),kmadd(SQR(n1),kmadd(R4p1112,rm1,kmadd(R4p1211,rm1,kmadd(R4p1213,rm3,kmadd(R4p1312,rm3,kmul(R4p1212,kmul(rm2,ToReal(2))))))),kmadd(SQR(n2),kmadd(R4p2122,rm1,kmadd(R4p2221,rm1,kmadd(R4p2223,rm3,kmadd(R4p2322,rm3,kmul(R4p2222,kmul(rm2,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm1,kmadd(R4p1221,rm1,kmadd(R4p2112,rm1,kmadd(R4p2211,rm1,kmadd(R4p1223,rm3,kmadd(R4p1322,rm3,kmadd(R4p2213,rm3,kmadd(R4p2312,rm3,kmadd(R4p1222,kmul(rm2,ToReal(2)),kmul(R4p2212,kmul(rm2,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1132,rm1,kmadd(R4p1231,rm1,kmadd(R4p3112,rm1,kmadd(R4p3211,rm1,kmadd(R4p1233,rm3,kmadd(R4p1332,rm3,kmadd(R4p3213,rm3,kmadd(R4p3312,rm3,kmadd(R4p1232,kmul(rm2,ToReal(2)),kmul(R4p3212,kmul(rm2,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm1,kmadd(R4p2231,rm1,kmadd(R4p3122,rm1,kmadd(R4p3221,rm1,kmadd(R4p2233,rm3,kmadd(R4p2332,rm3,kmadd(R4p3223,rm3,kmadd(R4p3322,rm3,kmadd(R4p2232,kmul(rm2,ToReal(2)),kmul(R4p3222,kmul(rm2,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro223,Ro322),kmul(rm2,kmul(Ro222,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm2,kmul(Ro232,ToReal(4)))))))))))))))))))))))))); + knmadd(im1,kmadd(R4p3132,kmul(rm2,kmul(n3,n3)),kmadd(R4p3231,kmul(rm2,kmul(n3,n3)),kmadd(R4p3133,kmul(rm3,kmul(n3,n3)),kmadd(R4p3331,kmul(rm3,kmul(n3,n3)),kmadd(rm2,kmul(Rojo12,kmul(nn,nn)),kmadd(rm3,kmul(Rojo13,kmul(nn,nn)),kmadd(rm2,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo31,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro132,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro133,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro231,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro331,ToReal(2)))),kmadd(R4p3131,kmul(rm1,kmul(kmul(n3,n3),ToReal(2))),kmadd(rm1,kmul(Rojo11,kmul(kmul(nn,nn),ToReal(2))),kmadd(kmul(n1,n1),kmadd(R4p1112,rm2,kmadd(R4p1211,rm2,kmadd(R4p1113,rm3,kmadd(R4p1311,rm3,kmul(R4p1111,kmul(rm1,ToReal(2))))))),kmadd(kmul(n2,n2),kmadd(R4p2122,rm2,kmadd(R4p2221,rm2,kmadd(R4p2123,rm3,kmadd(R4p2321,rm3,kmul(R4p2121,kmul(rm1,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm2,kmadd(R4p1221,rm2,kmadd(R4p2112,rm2,kmadd(R4p2211,rm2,kmadd(R4p1123,rm3,kmadd(R4p1321,rm3,kmadd(R4p2113,rm3,kmadd(R4p2311,rm3,kmadd(R4p1121,kmul(rm1,ToReal(2)),kmul(R4p2111,kmul(rm1,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1132,rm2,kmadd(R4p1231,rm2,kmadd(R4p3112,rm2,kmadd(R4p3211,rm2,kmadd(R4p1133,rm3,kmadd(R4p1331,rm3,kmadd(R4p3113,rm3,kmadd(R4p3311,rm3,kmadd(R4p1131,kmul(rm1,ToReal(2)),kmul(R4p3111,kmul(rm1,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm2,kmadd(R4p2231,rm2,kmadd(R4p3122,rm2,kmadd(R4p3221,rm2,kmadd(R4p2133,rm3,kmadd(R4p2331,rm3,kmadd(R4p3123,rm3,kmadd(R4p3321,rm3,kmadd(R4p2131,kmul(rm1,ToReal(2)),kmul(R4p3121,kmul(rm1,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm1,kmul(Ro131,ToReal(4))))))))))))))))))))))),kmadd(im3,kmadd(R4p3133,kmul(rm1,kmul(n3,n3)),kmadd(R4p3331,kmul(rm1,kmul(n3,n3)),kmadd(R4p3233,kmul(rm2,kmul(n3,n3)),kmadd(R4p3332,kmul(rm2,kmul(n3,n3)),kmadd(rm1,kmul(Rojo13,kmul(nn,nn)),kmadd(rm2,kmul(Rojo23,kmul(nn,nn)),kmadd(rm1,kmul(Rojo31,kmul(nn,nn)),kmadd(rm2,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro133,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro233,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro331,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro332,ToReal(2)))),kmadd(R4p3333,kmul(rm3,kmul(kmul(n3,n3),ToReal(2))),kmadd(rm3,kmul(Rojo33,kmul(kmul(nn,nn),ToReal(2))),kmadd(kmul(n1,n1),kmadd(R4p1113,rm1,kmadd(R4p1311,rm1,kmadd(R4p1213,rm2,kmadd(R4p1312,rm2,kmul(R4p1313,kmul(rm3,ToReal(2))))))),kmadd(kmul(n2,n2),kmadd(R4p2123,rm1,kmadd(R4p2321,rm1,kmadd(R4p2223,rm2,kmadd(R4p2322,rm2,kmul(R4p2323,kmul(rm3,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1123,rm1,kmadd(R4p1321,rm1,kmadd(R4p2113,rm1,kmadd(R4p2311,rm1,kmadd(R4p1223,rm2,kmadd(R4p1322,rm2,kmadd(R4p2213,rm2,kmadd(R4p2312,rm2,kmadd(R4p1323,kmul(rm3,ToReal(2)),kmul(R4p2313,kmul(rm3,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1133,rm1,kmadd(R4p1331,rm1,kmadd(R4p3113,rm1,kmadd(R4p3311,rm1,kmadd(R4p1233,rm2,kmadd(R4p1332,rm2,kmadd(R4p3213,rm2,kmadd(R4p3312,rm2,kmadd(R4p1333,kmul(rm3,ToReal(2)),kmul(R4p3313,kmul(rm3,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2133,rm1,kmadd(R4p2331,rm1,kmadd(R4p3123,rm1,kmadd(R4p3321,rm1,kmadd(R4p2233,rm2,kmadd(R4p2332,rm2,kmadd(R4p3223,rm2,kmadd(R4p3322,rm2,kmadd(R4p2333,kmul(rm3,ToReal(2)),kmul(R4p3323,kmul(rm3,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm3,kmul(Ro333,ToReal(4))))))))))))))))))))))),kmul(im2,kmadd(R4p3132,kmul(rm1,kmul(n3,n3)),kmadd(R4p3231,kmul(rm1,kmul(n3,n3)),kmadd(R4p3233,kmul(rm3,kmul(n3,n3)),kmadd(R4p3332,kmul(rm3,kmul(n3,n3)),kmadd(rm1,kmul(Rojo12,kmul(nn,nn)),kmadd(rm1,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo23,kmul(nn,nn)),kmadd(rm3,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro132,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro231,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro233,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro332,ToReal(2)))),kmadd(R4p3232,kmul(rm2,kmul(kmul(n3,n3),ToReal(2))),kmadd(rm2,kmul(Rojo22,kmul(kmul(nn,nn),ToReal(2))),kmadd(kmul(n1,n1),kmadd(R4p1112,rm1,kmadd(R4p1211,rm1,kmadd(R4p1213,rm3,kmadd(R4p1312,rm3,kmul(R4p1212,kmul(rm2,ToReal(2))))))),kmadd(kmul(n2,n2),kmadd(R4p2122,rm1,kmadd(R4p2221,rm1,kmadd(R4p2223,rm3,kmadd(R4p2322,rm3,kmul(R4p2222,kmul(rm2,ToReal(2))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm1,kmadd(R4p1221,rm1,kmadd(R4p2112,rm1,kmadd(R4p2211,rm1,kmadd(R4p1223,rm3,kmadd(R4p1322,rm3,kmadd(R4p2213,rm3,kmadd(R4p2312,rm3,kmadd(R4p1222,kmul(rm2,ToReal(2)),kmul(R4p2212,kmul(rm2,ToReal(2)))))))))))),kmadd(n3,kmadd(R4p1132,rm1,kmadd(R4p1231,rm1,kmadd(R4p3112,rm1,kmadd(R4p3211,rm1,kmadd(R4p1233,rm3,kmadd(R4p1332,rm3,kmadd(R4p3213,rm3,kmadd(R4p3312,rm3,kmadd(R4p1232,kmul(rm2,ToReal(2)),kmul(R4p3212,kmul(rm2,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm1,kmadd(R4p2231,rm1,kmadd(R4p3122,rm1,kmadd(R4p3221,rm1,kmadd(R4p2233,rm3,kmadd(R4p2332,rm3,kmadd(R4p3223,rm3,kmadd(R4p3322,rm3,kmadd(R4p2232,kmul(rm2,ToReal(2)),kmul(R4p3222,kmul(rm2,ToReal(2)))))))))))),kmul(nn,kmul(ToReal(2),kmadd(rm1,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro223,Ro322),kmul(rm2,kmul(Ro222,ToReal(2)))))))),kmul(n3,kmul(nn,kmul(rm2,kmul(Ro232,ToReal(4)))))))))))))))))))))))))); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); |