/* File produced by Kranc */ #define KRANC_C #include #include #include #include #include #include "cctk.h" #include "cctk_Arguments.h" #include "cctk_Parameters.h" #include "GenericFD.h" #include "Differencing.h" #include "cctk_Loop.h" #include "loopcontrol.h" #include "vectors.h" /* Define macros used in calculations */ #define INITVALUE (42) #define QAD(x) (SQR(SQR(x))) #define INV(x) (kdiv(ToReal(1.0),x)) #define SQR(x) (kmul(x,x)) #define CUB(x) (kmul(x,SQR(x))) extern "C" void psi4_calc_2nd_SelectBCs(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; CCTK_INT ierr = 0; ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4r_group","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4r_group."); return; } static void psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const imin[3], int const imax[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; /* Include user-supplied include files */ /* Initialise finite differencing variables */ ptrdiff_t const di = 1; ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); ptrdiff_t const cdi = sizeof(CCTK_REAL) * di; ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj; ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk; CCTK_REAL_VEC const dx = ToReal(CCTK_DELTA_SPACE(0)); CCTK_REAL_VEC const dy = ToReal(CCTK_DELTA_SPACE(1)); CCTK_REAL_VEC const dz = ToReal(CCTK_DELTA_SPACE(2)); CCTK_REAL_VEC const dt = ToReal(CCTK_DELTA_TIME); CCTK_REAL_VEC const t = ToReal(cctk_time); CCTK_REAL_VEC const dxi = INV(dx); CCTK_REAL_VEC const dyi = INV(dy); CCTK_REAL_VEC const dzi = INV(dz); CCTK_REAL_VEC const khalf = ToReal(0.5); CCTK_REAL_VEC const kthird = ToReal(1.0/3.0); CCTK_REAL_VEC const ktwothird = ToReal(2.0/3.0); CCTK_REAL_VEC const kfourthird = ToReal(4.0/3.0); CCTK_REAL_VEC const keightthird = ToReal(8.0/3.0); CCTK_REAL_VEC const hdxi = kmul(ToReal(0.5), dxi); CCTK_REAL_VEC const hdyi = kmul(ToReal(0.5), dyi); 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)); /* Jacobian variable pointers */ bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) && strlen(jacobian_group) > 0; if (use_jacobian && strlen(jacobian_derivative_group) == 0) { CCTK_WARN (1, "GenericFD::jacobian_group and GenericFD::jacobian_derivative_group must both be set to valid group names"); } CCTK_REAL const *restrict jacobian_ptrs[9]; if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_group, 9, jacobian_ptrs); CCTK_REAL const *restrict const J11 = use_jacobian ? jacobian_ptrs[0] : 0; CCTK_REAL const *restrict const J12 = use_jacobian ? jacobian_ptrs[1] : 0; CCTK_REAL const *restrict const J13 = use_jacobian ? jacobian_ptrs[2] : 0; CCTK_REAL const *restrict const J21 = use_jacobian ? jacobian_ptrs[3] : 0; CCTK_REAL const *restrict const J22 = use_jacobian ? jacobian_ptrs[4] : 0; CCTK_REAL const *restrict const J23 = use_jacobian ? jacobian_ptrs[5] : 0; CCTK_REAL const *restrict const J31 = use_jacobian ? jacobian_ptrs[6] : 0; CCTK_REAL const *restrict const J32 = use_jacobian ? jacobian_ptrs[7] : 0; CCTK_REAL const *restrict const J33 = use_jacobian ? jacobian_ptrs[8] : 0; CCTK_REAL const *restrict jacobian_derivative_ptrs[18]; if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_derivative_group, 18, jacobian_derivative_ptrs); CCTK_REAL const *restrict const dJ111 = use_jacobian ? jacobian_derivative_ptrs[0] : 0; CCTK_REAL const *restrict const dJ112 = use_jacobian ? jacobian_derivative_ptrs[1] : 0; CCTK_REAL const *restrict const dJ113 = use_jacobian ? jacobian_derivative_ptrs[2] : 0; CCTK_REAL const *restrict const dJ122 = use_jacobian ? jacobian_derivative_ptrs[3] : 0; CCTK_REAL const *restrict const dJ123 = use_jacobian ? jacobian_derivative_ptrs[4] : 0; CCTK_REAL const *restrict const dJ133 = use_jacobian ? jacobian_derivative_ptrs[5] : 0; CCTK_REAL const *restrict const dJ211 = use_jacobian ? jacobian_derivative_ptrs[6] : 0; CCTK_REAL const *restrict const dJ212 = use_jacobian ? jacobian_derivative_ptrs[7] : 0; CCTK_REAL const *restrict const dJ213 = use_jacobian ? jacobian_derivative_ptrs[8] : 0; CCTK_REAL const *restrict const dJ222 = use_jacobian ? jacobian_derivative_ptrs[9] : 0; CCTK_REAL const *restrict const dJ223 = use_jacobian ? jacobian_derivative_ptrs[10] : 0; CCTK_REAL const *restrict const dJ233 = use_jacobian ? jacobian_derivative_ptrs[11] : 0; CCTK_REAL const *restrict const dJ311 = use_jacobian ? jacobian_derivative_ptrs[12] : 0; CCTK_REAL const *restrict const dJ312 = use_jacobian ? jacobian_derivative_ptrs[13] : 0; CCTK_REAL const *restrict const dJ313 = use_jacobian ? jacobian_derivative_ptrs[14] : 0; CCTK_REAL const *restrict const dJ322 = use_jacobian ? jacobian_derivative_ptrs[15] : 0; CCTK_REAL const *restrict const dJ323 = use_jacobian ? jacobian_derivative_ptrs[16] : 0; CCTK_REAL const *restrict const dJ333 = use_jacobian ? jacobian_derivative_ptrs[17] : 0; /* Assign local copies of arrays functions */ /* Calculate temporaries and arrays functions */ /* Copy local copies back to grid functions */ /* Loop over the grid points */ #pragma omp parallel LC_LOOP3VEC(psi4_calc_2nd, 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_REAL_VEC_SIZE) { ptrdiff_t const index = di*i + dj*j + dk*k; /* Assign local copies of grid functions */ CCTK_REAL_VEC gxxL = vec_load(gxx[index]); CCTK_REAL_VEC gxyL = vec_load(gxy[index]); CCTK_REAL_VEC gxzL = vec_load(gxz[index]); CCTK_REAL_VEC gyyL = vec_load(gyy[index]); CCTK_REAL_VEC gyzL = vec_load(gyz[index]); CCTK_REAL_VEC gzzL = vec_load(gzz[index]); CCTK_REAL_VEC kxxL = vec_load(kxx[index]); CCTK_REAL_VEC kxyL = vec_load(kxy[index]); CCTK_REAL_VEC kxzL = vec_load(kxz[index]); CCTK_REAL_VEC kyyL = vec_load(kyy[index]); CCTK_REAL_VEC kyzL = vec_load(kyz[index]); CCTK_REAL_VEC kzzL = vec_load(kzz[index]); CCTK_REAL_VEC xL = vec_load(x[index]); CCTK_REAL_VEC yL = vec_load(y[index]); CCTK_REAL_VEC zL = vec_load(z[index]); CCTK_REAL_VEC dJ111L, dJ112L, dJ113L, dJ122L, dJ123L, dJ133L, dJ211L, dJ212L, dJ213L, dJ222L, dJ223L, dJ233L, dJ311L, dJ312L, dJ313L, dJ322L, dJ323L, dJ333L, J11L, J12L, J13L, J21L, J22L, J23L, J31L, J32L, J33L; if (use_jacobian) { dJ111L = vec_load(dJ111[index]); dJ112L = vec_load(dJ112[index]); dJ113L = vec_load(dJ113[index]); dJ122L = vec_load(dJ122[index]); dJ123L = vec_load(dJ123[index]); dJ133L = vec_load(dJ133[index]); dJ211L = vec_load(dJ211[index]); dJ212L = vec_load(dJ212[index]); dJ213L = vec_load(dJ213[index]); dJ222L = vec_load(dJ222[index]); dJ223L = vec_load(dJ223[index]); dJ233L = vec_load(dJ233[index]); dJ311L = vec_load(dJ311[index]); dJ312L = vec_load(dJ312[index]); dJ313L = vec_load(dJ313[index]); dJ322L = vec_load(dJ322[index]); dJ323L = vec_load(dJ323[index]); dJ333L = vec_load(dJ333[index]); J11L = vec_load(J11[index]); J12L = vec_load(J12[index]); J13L = vec_load(J13[index]); J21L = vec_load(J21[index]); J22L = vec_load(J22[index]); J23L = vec_load(J23[index]); J31L = vec_load(J31[index]); J32L = vec_load(J32[index]); J33L = vec_load(J33[index]); } /* Include user supplied include files */ /* Precompute derivatives */ CCTK_REAL_VEC PDstandard2nd1gxx; CCTK_REAL_VEC PDstandard2nd2gxx; CCTK_REAL_VEC PDstandard2nd3gxx; CCTK_REAL_VEC PDstandard2nd11gxx; CCTK_REAL_VEC PDstandard2nd22gxx; CCTK_REAL_VEC PDstandard2nd33gxx; CCTK_REAL_VEC PDstandard2nd12gxx; CCTK_REAL_VEC PDstandard2nd13gxx; CCTK_REAL_VEC PDstandard2nd23gxx; CCTK_REAL_VEC PDstandard2nd1gxy; CCTK_REAL_VEC PDstandard2nd2gxy; CCTK_REAL_VEC PDstandard2nd3gxy; CCTK_REAL_VEC PDstandard2nd11gxy; CCTK_REAL_VEC PDstandard2nd22gxy; CCTK_REAL_VEC PDstandard2nd33gxy; CCTK_REAL_VEC PDstandard2nd12gxy; CCTK_REAL_VEC PDstandard2nd13gxy; CCTK_REAL_VEC PDstandard2nd23gxy; CCTK_REAL_VEC PDstandard2nd1gxz; CCTK_REAL_VEC PDstandard2nd2gxz; CCTK_REAL_VEC PDstandard2nd3gxz; CCTK_REAL_VEC PDstandard2nd11gxz; CCTK_REAL_VEC PDstandard2nd22gxz; CCTK_REAL_VEC PDstandard2nd33gxz; CCTK_REAL_VEC PDstandard2nd12gxz; CCTK_REAL_VEC PDstandard2nd13gxz; CCTK_REAL_VEC PDstandard2nd23gxz; CCTK_REAL_VEC PDstandard2nd1gyy; CCTK_REAL_VEC PDstandard2nd2gyy; CCTK_REAL_VEC PDstandard2nd3gyy; CCTK_REAL_VEC PDstandard2nd11gyy; CCTK_REAL_VEC PDstandard2nd22gyy; CCTK_REAL_VEC PDstandard2nd33gyy; CCTK_REAL_VEC PDstandard2nd12gyy; CCTK_REAL_VEC PDstandard2nd13gyy; CCTK_REAL_VEC PDstandard2nd23gyy; CCTK_REAL_VEC PDstandard2nd1gyz; CCTK_REAL_VEC PDstandard2nd2gyz; CCTK_REAL_VEC PDstandard2nd3gyz; CCTK_REAL_VEC PDstandard2nd11gyz; CCTK_REAL_VEC PDstandard2nd22gyz; CCTK_REAL_VEC PDstandard2nd33gyz; CCTK_REAL_VEC PDstandard2nd12gyz; CCTK_REAL_VEC PDstandard2nd13gyz; CCTK_REAL_VEC PDstandard2nd23gyz; CCTK_REAL_VEC PDstandard2nd1gzz; CCTK_REAL_VEC PDstandard2nd2gzz; CCTK_REAL_VEC PDstandard2nd3gzz; CCTK_REAL_VEC PDstandard2nd11gzz; CCTK_REAL_VEC PDstandard2nd22gzz; CCTK_REAL_VEC PDstandard2nd33gzz; CCTK_REAL_VEC PDstandard2nd12gzz; CCTK_REAL_VEC PDstandard2nd13gzz; CCTK_REAL_VEC PDstandard2nd23gzz; CCTK_REAL_VEC PDstandard2nd1kxx; CCTK_REAL_VEC PDstandard2nd2kxx; CCTK_REAL_VEC PDstandard2nd3kxx; CCTK_REAL_VEC PDstandard2nd1kxy; CCTK_REAL_VEC PDstandard2nd2kxy; CCTK_REAL_VEC PDstandard2nd3kxy; CCTK_REAL_VEC PDstandard2nd1kxz; CCTK_REAL_VEC PDstandard2nd2kxz; CCTK_REAL_VEC PDstandard2nd3kxz; CCTK_REAL_VEC PDstandard2nd1kyy; CCTK_REAL_VEC PDstandard2nd2kyy; CCTK_REAL_VEC PDstandard2nd3kyy; CCTK_REAL_VEC PDstandard2nd1kyz; CCTK_REAL_VEC PDstandard2nd2kyz; CCTK_REAL_VEC PDstandard2nd3kyz; CCTK_REAL_VEC PDstandard2nd1kzz; CCTK_REAL_VEC PDstandard2nd2kzz; CCTK_REAL_VEC PDstandard2nd3kzz; switch(fdOrder) { case 2: PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); PDstandard2nd11gxx = PDstandard2nd11(&gxx[index]); PDstandard2nd22gxx = PDstandard2nd22(&gxx[index]); PDstandard2nd33gxx = PDstandard2nd33(&gxx[index]); PDstandard2nd12gxx = PDstandard2nd12(&gxx[index]); PDstandard2nd13gxx = PDstandard2nd13(&gxx[index]); PDstandard2nd23gxx = PDstandard2nd23(&gxx[index]); PDstandard2nd1gxy = PDstandard2nd1(&gxy[index]); PDstandard2nd2gxy = PDstandard2nd2(&gxy[index]); PDstandard2nd3gxy = PDstandard2nd3(&gxy[index]); PDstandard2nd11gxy = PDstandard2nd11(&gxy[index]); PDstandard2nd22gxy = PDstandard2nd22(&gxy[index]); PDstandard2nd33gxy = PDstandard2nd33(&gxy[index]); PDstandard2nd12gxy = PDstandard2nd12(&gxy[index]); PDstandard2nd13gxy = PDstandard2nd13(&gxy[index]); PDstandard2nd23gxy = PDstandard2nd23(&gxy[index]); PDstandard2nd1gxz = PDstandard2nd1(&gxz[index]); PDstandard2nd2gxz = PDstandard2nd2(&gxz[index]); PDstandard2nd3gxz = PDstandard2nd3(&gxz[index]); PDstandard2nd11gxz = PDstandard2nd11(&gxz[index]); PDstandard2nd22gxz = PDstandard2nd22(&gxz[index]); PDstandard2nd33gxz = PDstandard2nd33(&gxz[index]); PDstandard2nd12gxz = PDstandard2nd12(&gxz[index]); PDstandard2nd13gxz = PDstandard2nd13(&gxz[index]); PDstandard2nd23gxz = PDstandard2nd23(&gxz[index]); PDstandard2nd1gyy = PDstandard2nd1(&gyy[index]); PDstandard2nd2gyy = PDstandard2nd2(&gyy[index]); PDstandard2nd3gyy = PDstandard2nd3(&gyy[index]); PDstandard2nd11gyy = PDstandard2nd11(&gyy[index]); PDstandard2nd22gyy = PDstandard2nd22(&gyy[index]); PDstandard2nd33gyy = PDstandard2nd33(&gyy[index]); PDstandard2nd12gyy = PDstandard2nd12(&gyy[index]); PDstandard2nd13gyy = PDstandard2nd13(&gyy[index]); PDstandard2nd23gyy = PDstandard2nd23(&gyy[index]); PDstandard2nd1gyz = PDstandard2nd1(&gyz[index]); PDstandard2nd2gyz = PDstandard2nd2(&gyz[index]); PDstandard2nd3gyz = PDstandard2nd3(&gyz[index]); PDstandard2nd11gyz = PDstandard2nd11(&gyz[index]); PDstandard2nd22gyz = PDstandard2nd22(&gyz[index]); PDstandard2nd33gyz = PDstandard2nd33(&gyz[index]); PDstandard2nd12gyz = PDstandard2nd12(&gyz[index]); PDstandard2nd13gyz = PDstandard2nd13(&gyz[index]); PDstandard2nd23gyz = PDstandard2nd23(&gyz[index]); PDstandard2nd1gzz = PDstandard2nd1(&gzz[index]); PDstandard2nd2gzz = PDstandard2nd2(&gzz[index]); PDstandard2nd3gzz = PDstandard2nd3(&gzz[index]); PDstandard2nd11gzz = PDstandard2nd11(&gzz[index]); PDstandard2nd22gzz = PDstandard2nd22(&gzz[index]); PDstandard2nd33gzz = PDstandard2nd33(&gzz[index]); PDstandard2nd12gzz = PDstandard2nd12(&gzz[index]); PDstandard2nd13gzz = PDstandard2nd13(&gzz[index]); PDstandard2nd23gzz = PDstandard2nd23(&gzz[index]); PDstandard2nd1kxx = PDstandard2nd1(&kxx[index]); PDstandard2nd2kxx = PDstandard2nd2(&kxx[index]); PDstandard2nd3kxx = PDstandard2nd3(&kxx[index]); PDstandard2nd1kxy = PDstandard2nd1(&kxy[index]); PDstandard2nd2kxy = PDstandard2nd2(&kxy[index]); PDstandard2nd3kxy = PDstandard2nd3(&kxy[index]); PDstandard2nd1kxz = PDstandard2nd1(&kxz[index]); PDstandard2nd2kxz = PDstandard2nd2(&kxz[index]); PDstandard2nd3kxz = PDstandard2nd3(&kxz[index]); PDstandard2nd1kyy = PDstandard2nd1(&kyy[index]); PDstandard2nd2kyy = PDstandard2nd2(&kyy[index]); PDstandard2nd3kyy = PDstandard2nd3(&kyy[index]); PDstandard2nd1kyz = PDstandard2nd1(&kyz[index]); PDstandard2nd2kyz = PDstandard2nd2(&kyz[index]); PDstandard2nd3kyz = PDstandard2nd3(&kyz[index]); PDstandard2nd1kzz = PDstandard2nd1(&kzz[index]); PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; case 4: PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); PDstandard2nd11gxx = PDstandard2nd11(&gxx[index]); PDstandard2nd22gxx = PDstandard2nd22(&gxx[index]); PDstandard2nd33gxx = PDstandard2nd33(&gxx[index]); PDstandard2nd12gxx = PDstandard2nd12(&gxx[index]); PDstandard2nd13gxx = PDstandard2nd13(&gxx[index]); PDstandard2nd23gxx = PDstandard2nd23(&gxx[index]); PDstandard2nd1gxy = PDstandard2nd1(&gxy[index]); PDstandard2nd2gxy = PDstandard2nd2(&gxy[index]); PDstandard2nd3gxy = PDstandard2nd3(&gxy[index]); PDstandard2nd11gxy = PDstandard2nd11(&gxy[index]); PDstandard2nd22gxy = PDstandard2nd22(&gxy[index]); PDstandard2nd33gxy = PDstandard2nd33(&gxy[index]); PDstandard2nd12gxy = PDstandard2nd12(&gxy[index]); PDstandard2nd13gxy = PDstandard2nd13(&gxy[index]); PDstandard2nd23gxy = PDstandard2nd23(&gxy[index]); PDstandard2nd1gxz = PDstandard2nd1(&gxz[index]); PDstandard2nd2gxz = PDstandard2nd2(&gxz[index]); PDstandard2nd3gxz = PDstandard2nd3(&gxz[index]); PDstandard2nd11gxz = PDstandard2nd11(&gxz[index]); PDstandard2nd22gxz = PDstandard2nd22(&gxz[index]); PDstandard2nd33gxz = PDstandard2nd33(&gxz[index]); PDstandard2nd12gxz = PDstandard2nd12(&gxz[index]); PDstandard2nd13gxz = PDstandard2nd13(&gxz[index]); PDstandard2nd23gxz = PDstandard2nd23(&gxz[index]); PDstandard2nd1gyy = PDstandard2nd1(&gyy[index]); PDstandard2nd2gyy = PDstandard2nd2(&gyy[index]); PDstandard2nd3gyy = PDstandard2nd3(&gyy[index]); PDstandard2nd11gyy = PDstandard2nd11(&gyy[index]); PDstandard2nd22gyy = PDstandard2nd22(&gyy[index]); PDstandard2nd33gyy = PDstandard2nd33(&gyy[index]); PDstandard2nd12gyy = PDstandard2nd12(&gyy[index]); PDstandard2nd13gyy = PDstandard2nd13(&gyy[index]); PDstandard2nd23gyy = PDstandard2nd23(&gyy[index]); PDstandard2nd1gyz = PDstandard2nd1(&gyz[index]); PDstandard2nd2gyz = PDstandard2nd2(&gyz[index]); PDstandard2nd3gyz = PDstandard2nd3(&gyz[index]); PDstandard2nd11gyz = PDstandard2nd11(&gyz[index]); PDstandard2nd22gyz = PDstandard2nd22(&gyz[index]); PDstandard2nd33gyz = PDstandard2nd33(&gyz[index]); PDstandard2nd12gyz = PDstandard2nd12(&gyz[index]); PDstandard2nd13gyz = PDstandard2nd13(&gyz[index]); PDstandard2nd23gyz = PDstandard2nd23(&gyz[index]); PDstandard2nd1gzz = PDstandard2nd1(&gzz[index]); PDstandard2nd2gzz = PDstandard2nd2(&gzz[index]); PDstandard2nd3gzz = PDstandard2nd3(&gzz[index]); PDstandard2nd11gzz = PDstandard2nd11(&gzz[index]); PDstandard2nd22gzz = PDstandard2nd22(&gzz[index]); PDstandard2nd33gzz = PDstandard2nd33(&gzz[index]); PDstandard2nd12gzz = PDstandard2nd12(&gzz[index]); PDstandard2nd13gzz = PDstandard2nd13(&gzz[index]); PDstandard2nd23gzz = PDstandard2nd23(&gzz[index]); PDstandard2nd1kxx = PDstandard2nd1(&kxx[index]); PDstandard2nd2kxx = PDstandard2nd2(&kxx[index]); PDstandard2nd3kxx = PDstandard2nd3(&kxx[index]); PDstandard2nd1kxy = PDstandard2nd1(&kxy[index]); PDstandard2nd2kxy = PDstandard2nd2(&kxy[index]); PDstandard2nd3kxy = PDstandard2nd3(&kxy[index]); PDstandard2nd1kxz = PDstandard2nd1(&kxz[index]); PDstandard2nd2kxz = PDstandard2nd2(&kxz[index]); PDstandard2nd3kxz = PDstandard2nd3(&kxz[index]); PDstandard2nd1kyy = PDstandard2nd1(&kyy[index]); PDstandard2nd2kyy = PDstandard2nd2(&kyy[index]); PDstandard2nd3kyy = PDstandard2nd3(&kyy[index]); PDstandard2nd1kyz = PDstandard2nd1(&kyz[index]); PDstandard2nd2kyz = PDstandard2nd2(&kyz[index]); PDstandard2nd3kyz = PDstandard2nd3(&kyz[index]); PDstandard2nd1kzz = PDstandard2nd1(&kzz[index]); PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; case 6: PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); PDstandard2nd11gxx = PDstandard2nd11(&gxx[index]); PDstandard2nd22gxx = PDstandard2nd22(&gxx[index]); PDstandard2nd33gxx = PDstandard2nd33(&gxx[index]); PDstandard2nd12gxx = PDstandard2nd12(&gxx[index]); PDstandard2nd13gxx = PDstandard2nd13(&gxx[index]); PDstandard2nd23gxx = PDstandard2nd23(&gxx[index]); PDstandard2nd1gxy = PDstandard2nd1(&gxy[index]); PDstandard2nd2gxy = PDstandard2nd2(&gxy[index]); PDstandard2nd3gxy = PDstandard2nd3(&gxy[index]); PDstandard2nd11gxy = PDstandard2nd11(&gxy[index]); PDstandard2nd22gxy = PDstandard2nd22(&gxy[index]); PDstandard2nd33gxy = PDstandard2nd33(&gxy[index]); PDstandard2nd12gxy = PDstandard2nd12(&gxy[index]); PDstandard2nd13gxy = PDstandard2nd13(&gxy[index]); PDstandard2nd23gxy = PDstandard2nd23(&gxy[index]); PDstandard2nd1gxz = PDstandard2nd1(&gxz[index]); PDstandard2nd2gxz = PDstandard2nd2(&gxz[index]); PDstandard2nd3gxz = PDstandard2nd3(&gxz[index]); PDstandard2nd11gxz = PDstandard2nd11(&gxz[index]); PDstandard2nd22gxz = PDstandard2nd22(&gxz[index]); PDstandard2nd33gxz = PDstandard2nd33(&gxz[index]); PDstandard2nd12gxz = PDstandard2nd12(&gxz[index]); PDstandard2nd13gxz = PDstandard2nd13(&gxz[index]); PDstandard2nd23gxz = PDstandard2nd23(&gxz[index]); PDstandard2nd1gyy = PDstandard2nd1(&gyy[index]); PDstandard2nd2gyy = PDstandard2nd2(&gyy[index]); PDstandard2nd3gyy = PDstandard2nd3(&gyy[index]); PDstandard2nd11gyy = PDstandard2nd11(&gyy[index]); PDstandard2nd22gyy = PDstandard2nd22(&gyy[index]); PDstandard2nd33gyy = PDstandard2nd33(&gyy[index]); PDstandard2nd12gyy = PDstandard2nd12(&gyy[index]); PDstandard2nd13gyy = PDstandard2nd13(&gyy[index]); PDstandard2nd23gyy = PDstandard2nd23(&gyy[index]); PDstandard2nd1gyz = PDstandard2nd1(&gyz[index]); PDstandard2nd2gyz = PDstandard2nd2(&gyz[index]); PDstandard2nd3gyz = PDstandard2nd3(&gyz[index]); PDstandard2nd11gyz = PDstandard2nd11(&gyz[index]); PDstandard2nd22gyz = PDstandard2nd22(&gyz[index]); PDstandard2nd33gyz = PDstandard2nd33(&gyz[index]); PDstandard2nd12gyz = PDstandard2nd12(&gyz[index]); PDstandard2nd13gyz = PDstandard2nd13(&gyz[index]); PDstandard2nd23gyz = PDstandard2nd23(&gyz[index]); PDstandard2nd1gzz = PDstandard2nd1(&gzz[index]); PDstandard2nd2gzz = PDstandard2nd2(&gzz[index]); PDstandard2nd3gzz = PDstandard2nd3(&gzz[index]); PDstandard2nd11gzz = PDstandard2nd11(&gzz[index]); PDstandard2nd22gzz = PDstandard2nd22(&gzz[index]); PDstandard2nd33gzz = PDstandard2nd33(&gzz[index]); PDstandard2nd12gzz = PDstandard2nd12(&gzz[index]); PDstandard2nd13gzz = PDstandard2nd13(&gzz[index]); PDstandard2nd23gzz = PDstandard2nd23(&gzz[index]); PDstandard2nd1kxx = PDstandard2nd1(&kxx[index]); PDstandard2nd2kxx = PDstandard2nd2(&kxx[index]); PDstandard2nd3kxx = PDstandard2nd3(&kxx[index]); PDstandard2nd1kxy = PDstandard2nd1(&kxy[index]); PDstandard2nd2kxy = PDstandard2nd2(&kxy[index]); PDstandard2nd3kxy = PDstandard2nd3(&kxy[index]); PDstandard2nd1kxz = PDstandard2nd1(&kxz[index]); PDstandard2nd2kxz = PDstandard2nd2(&kxz[index]); PDstandard2nd3kxz = PDstandard2nd3(&kxz[index]); PDstandard2nd1kyy = PDstandard2nd1(&kyy[index]); PDstandard2nd2kyy = PDstandard2nd2(&kyy[index]); PDstandard2nd3kyy = PDstandard2nd3(&kyy[index]); PDstandard2nd1kyz = PDstandard2nd1(&kyz[index]); PDstandard2nd2kyz = PDstandard2nd2(&kyz[index]); PDstandard2nd3kyz = PDstandard2nd3(&kyz[index]); PDstandard2nd1kzz = PDstandard2nd1(&kzz[index]); PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; case 8: PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); PDstandard2nd11gxx = PDstandard2nd11(&gxx[index]); PDstandard2nd22gxx = PDstandard2nd22(&gxx[index]); PDstandard2nd33gxx = PDstandard2nd33(&gxx[index]); PDstandard2nd12gxx = PDstandard2nd12(&gxx[index]); PDstandard2nd13gxx = PDstandard2nd13(&gxx[index]); PDstandard2nd23gxx = PDstandard2nd23(&gxx[index]); PDstandard2nd1gxy = PDstandard2nd1(&gxy[index]); PDstandard2nd2gxy = PDstandard2nd2(&gxy[index]); PDstandard2nd3gxy = PDstandard2nd3(&gxy[index]); PDstandard2nd11gxy = PDstandard2nd11(&gxy[index]); PDstandard2nd22gxy = PDstandard2nd22(&gxy[index]); PDstandard2nd33gxy = PDstandard2nd33(&gxy[index]); PDstandard2nd12gxy = PDstandard2nd12(&gxy[index]); PDstandard2nd13gxy = PDstandard2nd13(&gxy[index]); PDstandard2nd23gxy = PDstandard2nd23(&gxy[index]); PDstandard2nd1gxz = PDstandard2nd1(&gxz[index]); PDstandard2nd2gxz = PDstandard2nd2(&gxz[index]); PDstandard2nd3gxz = PDstandard2nd3(&gxz[index]); PDstandard2nd11gxz = PDstandard2nd11(&gxz[index]); PDstandard2nd22gxz = PDstandard2nd22(&gxz[index]); PDstandard2nd33gxz = PDstandard2nd33(&gxz[index]); PDstandard2nd12gxz = PDstandard2nd12(&gxz[index]); PDstandard2nd13gxz = PDstandard2nd13(&gxz[index]); PDstandard2nd23gxz = PDstandard2nd23(&gxz[index]); PDstandard2nd1gyy = PDstandard2nd1(&gyy[index]); PDstandard2nd2gyy = PDstandard2nd2(&gyy[index]); PDstandard2nd3gyy = PDstandard2nd3(&gyy[index]); PDstandard2nd11gyy = PDstandard2nd11(&gyy[index]); PDstandard2nd22gyy = PDstandard2nd22(&gyy[index]); PDstandard2nd33gyy = PDstandard2nd33(&gyy[index]); PDstandard2nd12gyy = PDstandard2nd12(&gyy[index]); PDstandard2nd13gyy = PDstandard2nd13(&gyy[index]); PDstandard2nd23gyy = PDstandard2nd23(&gyy[index]); PDstandard2nd1gyz = PDstandard2nd1(&gyz[index]); PDstandard2nd2gyz = PDstandard2nd2(&gyz[index]); PDstandard2nd3gyz = PDstandard2nd3(&gyz[index]); PDstandard2nd11gyz = PDstandard2nd11(&gyz[index]); PDstandard2nd22gyz = PDstandard2nd22(&gyz[index]); PDstandard2nd33gyz = PDstandard2nd33(&gyz[index]); PDstandard2nd12gyz = PDstandard2nd12(&gyz[index]); PDstandard2nd13gyz = PDstandard2nd13(&gyz[index]); PDstandard2nd23gyz = PDstandard2nd23(&gyz[index]); PDstandard2nd1gzz = PDstandard2nd1(&gzz[index]); PDstandard2nd2gzz = PDstandard2nd2(&gzz[index]); PDstandard2nd3gzz = PDstandard2nd3(&gzz[index]); PDstandard2nd11gzz = PDstandard2nd11(&gzz[index]); PDstandard2nd22gzz = PDstandard2nd22(&gzz[index]); PDstandard2nd33gzz = PDstandard2nd33(&gzz[index]); PDstandard2nd12gzz = PDstandard2nd12(&gzz[index]); PDstandard2nd13gzz = PDstandard2nd13(&gzz[index]); PDstandard2nd23gzz = PDstandard2nd23(&gzz[index]); PDstandard2nd1kxx = PDstandard2nd1(&kxx[index]); PDstandard2nd2kxx = PDstandard2nd2(&kxx[index]); PDstandard2nd3kxx = PDstandard2nd3(&kxx[index]); PDstandard2nd1kxy = PDstandard2nd1(&kxy[index]); PDstandard2nd2kxy = PDstandard2nd2(&kxy[index]); PDstandard2nd3kxy = PDstandard2nd3(&kxy[index]); PDstandard2nd1kxz = PDstandard2nd1(&kxz[index]); PDstandard2nd2kxz = PDstandard2nd2(&kxz[index]); PDstandard2nd3kxz = PDstandard2nd3(&kxz[index]); PDstandard2nd1kyy = PDstandard2nd1(&kyy[index]); PDstandard2nd2kyy = PDstandard2nd2(&kyy[index]); PDstandard2nd3kyy = PDstandard2nd3(&kyy[index]); PDstandard2nd1kyz = PDstandard2nd1(&kyz[index]); PDstandard2nd2kyz = PDstandard2nd2(&kyz[index]); PDstandard2nd3kyz = PDstandard2nd3(&kyz[index]); PDstandard2nd1kzz = PDstandard2nd1(&kzz[index]); PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; } /* Calculate temporaries and grid functions */ CCTK_REAL_VEC JacPDstandard2nd11gxy; CCTK_REAL_VEC JacPDstandard2nd11gxz; CCTK_REAL_VEC JacPDstandard2nd11gyy; CCTK_REAL_VEC JacPDstandard2nd11gyz; CCTK_REAL_VEC JacPDstandard2nd11gzz; CCTK_REAL_VEC JacPDstandard2nd12gxx; CCTK_REAL_VEC JacPDstandard2nd12gxy; CCTK_REAL_VEC JacPDstandard2nd12gxz; CCTK_REAL_VEC JacPDstandard2nd12gyy; CCTK_REAL_VEC JacPDstandard2nd12gyz; CCTK_REAL_VEC JacPDstandard2nd12gzz; CCTK_REAL_VEC JacPDstandard2nd13gxx; CCTK_REAL_VEC JacPDstandard2nd13gxy; CCTK_REAL_VEC JacPDstandard2nd13gxz; CCTK_REAL_VEC JacPDstandard2nd13gyy; CCTK_REAL_VEC JacPDstandard2nd13gyz; CCTK_REAL_VEC JacPDstandard2nd13gzz; CCTK_REAL_VEC JacPDstandard2nd1gxx; CCTK_REAL_VEC JacPDstandard2nd1gxy; CCTK_REAL_VEC JacPDstandard2nd1gxz; CCTK_REAL_VEC JacPDstandard2nd1gyy; CCTK_REAL_VEC JacPDstandard2nd1gyz; CCTK_REAL_VEC JacPDstandard2nd1gzz; CCTK_REAL_VEC JacPDstandard2nd1kxy; CCTK_REAL_VEC JacPDstandard2nd1kxz; CCTK_REAL_VEC JacPDstandard2nd1kyy; CCTK_REAL_VEC JacPDstandard2nd1kyz; CCTK_REAL_VEC JacPDstandard2nd1kzz; CCTK_REAL_VEC JacPDstandard2nd21gxx; CCTK_REAL_VEC JacPDstandard2nd21gxy; CCTK_REAL_VEC JacPDstandard2nd21gxz; CCTK_REAL_VEC JacPDstandard2nd21gyy; CCTK_REAL_VEC JacPDstandard2nd21gyz; CCTK_REAL_VEC JacPDstandard2nd21gzz; CCTK_REAL_VEC JacPDstandard2nd22gxx; CCTK_REAL_VEC JacPDstandard2nd22gxy; CCTK_REAL_VEC JacPDstandard2nd22gxz; CCTK_REAL_VEC JacPDstandard2nd22gyz; CCTK_REAL_VEC JacPDstandard2nd22gzz; CCTK_REAL_VEC JacPDstandard2nd23gxx; CCTK_REAL_VEC JacPDstandard2nd23gxy; CCTK_REAL_VEC JacPDstandard2nd23gxz; CCTK_REAL_VEC JacPDstandard2nd23gyy; CCTK_REAL_VEC JacPDstandard2nd23gyz; CCTK_REAL_VEC JacPDstandard2nd23gzz; CCTK_REAL_VEC JacPDstandard2nd2gxx; CCTK_REAL_VEC JacPDstandard2nd2gxy; CCTK_REAL_VEC JacPDstandard2nd2gxz; CCTK_REAL_VEC JacPDstandard2nd2gyy; CCTK_REAL_VEC JacPDstandard2nd2gyz; CCTK_REAL_VEC JacPDstandard2nd2gzz; CCTK_REAL_VEC JacPDstandard2nd2kxx; CCTK_REAL_VEC JacPDstandard2nd2kxy; CCTK_REAL_VEC JacPDstandard2nd2kxz; CCTK_REAL_VEC JacPDstandard2nd2kyz; CCTK_REAL_VEC JacPDstandard2nd2kzz; CCTK_REAL_VEC JacPDstandard2nd31gxx; CCTK_REAL_VEC JacPDstandard2nd31gxy; CCTK_REAL_VEC JacPDstandard2nd31gxz; CCTK_REAL_VEC JacPDstandard2nd31gyy; CCTK_REAL_VEC JacPDstandard2nd31gyz; CCTK_REAL_VEC JacPDstandard2nd31gzz; CCTK_REAL_VEC JacPDstandard2nd32gxx; CCTK_REAL_VEC JacPDstandard2nd32gxy; CCTK_REAL_VEC JacPDstandard2nd32gxz; CCTK_REAL_VEC JacPDstandard2nd32gyy; CCTK_REAL_VEC JacPDstandard2nd32gyz; CCTK_REAL_VEC JacPDstandard2nd32gzz; CCTK_REAL_VEC JacPDstandard2nd33gxx; CCTK_REAL_VEC JacPDstandard2nd33gxy; CCTK_REAL_VEC JacPDstandard2nd33gxz; CCTK_REAL_VEC JacPDstandard2nd33gyy; CCTK_REAL_VEC JacPDstandard2nd33gyz; CCTK_REAL_VEC JacPDstandard2nd3gxx; CCTK_REAL_VEC JacPDstandard2nd3gxy; CCTK_REAL_VEC JacPDstandard2nd3gxz; CCTK_REAL_VEC JacPDstandard2nd3gyy; CCTK_REAL_VEC JacPDstandard2nd3gyz; CCTK_REAL_VEC JacPDstandard2nd3gzz; CCTK_REAL_VEC JacPDstandard2nd3kxx; CCTK_REAL_VEC JacPDstandard2nd3kxy; CCTK_REAL_VEC JacPDstandard2nd3kxz; CCTK_REAL_VEC JacPDstandard2nd3kyy; CCTK_REAL_VEC JacPDstandard2nd3kyz; if (use_jacobian) { JacPDstandard2nd1gxx = kmadd(J11L,PDstandard2nd1gxx,kmadd(J21L,PDstandard2nd2gxx,kmul(J31L,PDstandard2nd3gxx))); JacPDstandard2nd1gxy = kmadd(J11L,PDstandard2nd1gxy,kmadd(J21L,PDstandard2nd2gxy,kmul(J31L,PDstandard2nd3gxy))); JacPDstandard2nd1gxz = kmadd(J11L,PDstandard2nd1gxz,kmadd(J21L,PDstandard2nd2gxz,kmul(J31L,PDstandard2nd3gxz))); JacPDstandard2nd1gyy = kmadd(J11L,PDstandard2nd1gyy,kmadd(J21L,PDstandard2nd2gyy,kmul(J31L,PDstandard2nd3gyy))); JacPDstandard2nd1gyz = kmadd(J11L,PDstandard2nd1gyz,kmadd(J21L,PDstandard2nd2gyz,kmul(J31L,PDstandard2nd3gyz))); JacPDstandard2nd1gzz = kmadd(J11L,PDstandard2nd1gzz,kmadd(J21L,PDstandard2nd2gzz,kmul(J31L,PDstandard2nd3gzz))); JacPDstandard2nd1kxy = kmadd(J11L,PDstandard2nd1kxy,kmadd(J21L,PDstandard2nd2kxy,kmul(J31L,PDstandard2nd3kxy))); JacPDstandard2nd1kxz = kmadd(J11L,PDstandard2nd1kxz,kmadd(J21L,PDstandard2nd2kxz,kmul(J31L,PDstandard2nd3kxz))); JacPDstandard2nd1kyy = kmadd(J11L,PDstandard2nd1kyy,kmadd(J21L,PDstandard2nd2kyy,kmul(J31L,PDstandard2nd3kyy))); JacPDstandard2nd1kyz = kmadd(J11L,PDstandard2nd1kyz,kmadd(J21L,PDstandard2nd2kyz,kmul(J31L,PDstandard2nd3kyz))); JacPDstandard2nd1kzz = kmadd(J11L,PDstandard2nd1kzz,kmadd(J21L,PDstandard2nd2kzz,kmul(J31L,PDstandard2nd3kzz))); JacPDstandard2nd2gxx = kmadd(J12L,PDstandard2nd1gxx,kmadd(J22L,PDstandard2nd2gxx,kmul(J32L,PDstandard2nd3gxx))); JacPDstandard2nd2gxy = kmadd(J12L,PDstandard2nd1gxy,kmadd(J22L,PDstandard2nd2gxy,kmul(J32L,PDstandard2nd3gxy))); JacPDstandard2nd2gxz = kmadd(J12L,PDstandard2nd1gxz,kmadd(J22L,PDstandard2nd2gxz,kmul(J32L,PDstandard2nd3gxz))); JacPDstandard2nd2gyy = kmadd(J12L,PDstandard2nd1gyy,kmadd(J22L,PDstandard2nd2gyy,kmul(J32L,PDstandard2nd3gyy))); JacPDstandard2nd2gyz = kmadd(J12L,PDstandard2nd1gyz,kmadd(J22L,PDstandard2nd2gyz,kmul(J32L,PDstandard2nd3gyz))); JacPDstandard2nd2gzz = kmadd(J12L,PDstandard2nd1gzz,kmadd(J22L,PDstandard2nd2gzz,kmul(J32L,PDstandard2nd3gzz))); JacPDstandard2nd2kxx = kmadd(J12L,PDstandard2nd1kxx,kmadd(J22L,PDstandard2nd2kxx,kmul(J32L,PDstandard2nd3kxx))); JacPDstandard2nd2kxy = kmadd(J12L,PDstandard2nd1kxy,kmadd(J22L,PDstandard2nd2kxy,kmul(J32L,PDstandard2nd3kxy))); JacPDstandard2nd2kxz = kmadd(J12L,PDstandard2nd1kxz,kmadd(J22L,PDstandard2nd2kxz,kmul(J32L,PDstandard2nd3kxz))); JacPDstandard2nd2kyz = kmadd(J12L,PDstandard2nd1kyz,kmadd(J22L,PDstandard2nd2kyz,kmul(J32L,PDstandard2nd3kyz))); JacPDstandard2nd2kzz = kmadd(J12L,PDstandard2nd1kzz,kmadd(J22L,PDstandard2nd2kzz,kmul(J32L,PDstandard2nd3kzz))); JacPDstandard2nd3gxx = kmadd(J13L,PDstandard2nd1gxx,kmadd(J23L,PDstandard2nd2gxx,kmul(J33L,PDstandard2nd3gxx))); JacPDstandard2nd3gxy = kmadd(J13L,PDstandard2nd1gxy,kmadd(J23L,PDstandard2nd2gxy,kmul(J33L,PDstandard2nd3gxy))); JacPDstandard2nd3gxz = kmadd(J13L,PDstandard2nd1gxz,kmadd(J23L,PDstandard2nd2gxz,kmul(J33L,PDstandard2nd3gxz))); JacPDstandard2nd3gyy = kmadd(J13L,PDstandard2nd1gyy,kmadd(J23L,PDstandard2nd2gyy,kmul(J33L,PDstandard2nd3gyy))); JacPDstandard2nd3gyz = kmadd(J13L,PDstandard2nd1gyz,kmadd(J23L,PDstandard2nd2gyz,kmul(J33L,PDstandard2nd3gyz))); JacPDstandard2nd3gzz = kmadd(J13L,PDstandard2nd1gzz,kmadd(J23L,PDstandard2nd2gzz,kmul(J33L,PDstandard2nd3gzz))); JacPDstandard2nd3kxx = kmadd(J13L,PDstandard2nd1kxx,kmadd(J23L,PDstandard2nd2kxx,kmul(J33L,PDstandard2nd3kxx))); JacPDstandard2nd3kxy = kmadd(J13L,PDstandard2nd1kxy,kmadd(J23L,PDstandard2nd2kxy,kmul(J33L,PDstandard2nd3kxy))); JacPDstandard2nd3kxz = kmadd(J13L,PDstandard2nd1kxz,kmadd(J23L,PDstandard2nd2kxz,kmul(J33L,PDstandard2nd3kxz))); JacPDstandard2nd3kyy = kmadd(J13L,PDstandard2nd1kyy,kmadd(J23L,PDstandard2nd2kyy,kmul(J33L,PDstandard2nd3kyy))); JacPDstandard2nd3kyz = kmadd(J13L,PDstandard2nd1kyz,kmadd(J23L,PDstandard2nd2kyz,kmul(J33L,PDstandard2nd3kyz))); JacPDstandard2nd11gxy = kmadd(dJ111L,PDstandard2nd1gxy,kmadd(dJ211L,PDstandard2nd2gxy,kmadd(dJ311L,PDstandard2nd3gxy,kmadd(PDstandard2nd11gxy,SQR(J11L),kmadd(PDstandard2nd22gxy,SQR(J21L),kmadd(PDstandard2nd33gxy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy)),kmul(J21L,kmul(J31L,PDstandard2nd23gxy))),ToReal(2)))))))); JacPDstandard2nd11gxz = kmadd(dJ111L,PDstandard2nd1gxz,kmadd(dJ211L,PDstandard2nd2gxz,kmadd(dJ311L,PDstandard2nd3gxz,kmadd(PDstandard2nd11gxz,SQR(J11L),kmadd(PDstandard2nd22gxz,SQR(J21L),kmadd(PDstandard2nd33gxz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gxz,kmul(J31L,PDstandard2nd13gxz)),kmul(J21L,kmul(J31L,PDstandard2nd23gxz))),ToReal(2)))))))); JacPDstandard2nd11gyy = kmadd(dJ111L,PDstandard2nd1gyy,kmadd(dJ211L,PDstandard2nd2gyy,kmadd(dJ311L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,SQR(J11L),kmadd(PDstandard2nd22gyy,SQR(J21L),kmadd(PDstandard2nd33gyy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy)),kmul(J21L,kmul(J31L,PDstandard2nd23gyy))),ToReal(2)))))))); JacPDstandard2nd11gyz = kmadd(dJ111L,PDstandard2nd1gyz,kmadd(dJ211L,PDstandard2nd2gyz,kmadd(dJ311L,PDstandard2nd3gyz,kmadd(PDstandard2nd11gyz,SQR(J11L),kmadd(PDstandard2nd22gyz,SQR(J21L),kmadd(PDstandard2nd33gyz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz)),kmul(J21L,kmul(J31L,PDstandard2nd23gyz))),ToReal(2)))))))); JacPDstandard2nd11gzz = kmadd(dJ111L,PDstandard2nd1gzz,kmadd(dJ211L,PDstandard2nd2gzz,kmadd(dJ311L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,SQR(J11L),kmadd(PDstandard2nd22gzz,SQR(J21L),kmadd(PDstandard2nd33gzz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz)),kmul(J21L,kmul(J31L,PDstandard2nd23gzz))),ToReal(2)))))))); JacPDstandard2nd22gxx = kmadd(dJ122L,PDstandard2nd1gxx,kmadd(dJ222L,PDstandard2nd2gxx,kmadd(dJ322L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,SQR(J12L),kmadd(PDstandard2nd22gxx,SQR(J22L),kmadd(PDstandard2nd33gxx,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx)),kmul(J22L,kmul(J32L,PDstandard2nd23gxx))),ToReal(2)))))))); JacPDstandard2nd22gxy = kmadd(dJ122L,PDstandard2nd1gxy,kmadd(dJ222L,PDstandard2nd2gxy,kmadd(dJ322L,PDstandard2nd3gxy,kmadd(PDstandard2nd11gxy,SQR(J12L),kmadd(PDstandard2nd22gxy,SQR(J22L),kmadd(PDstandard2nd33gxy,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy)),kmul(J22L,kmul(J32L,PDstandard2nd23gxy))),ToReal(2)))))))); JacPDstandard2nd22gxz = kmadd(dJ122L,PDstandard2nd1gxz,kmadd(dJ222L,PDstandard2nd2gxz,kmadd(dJ322L,PDstandard2nd3gxz,kmadd(PDstandard2nd11gxz,SQR(J12L),kmadd(PDstandard2nd22gxz,SQR(J22L),kmadd(PDstandard2nd33gxz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz)),kmul(J22L,kmul(J32L,PDstandard2nd23gxz))),ToReal(2)))))))); JacPDstandard2nd22gyz = kmadd(dJ122L,PDstandard2nd1gyz,kmadd(dJ222L,PDstandard2nd2gyz,kmadd(dJ322L,PDstandard2nd3gyz,kmadd(PDstandard2nd11gyz,SQR(J12L),kmadd(PDstandard2nd22gyz,SQR(J22L),kmadd(PDstandard2nd33gyz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gyz,kmul(J32L,PDstandard2nd13gyz)),kmul(J22L,kmul(J32L,PDstandard2nd23gyz))),ToReal(2)))))))); JacPDstandard2nd22gzz = kmadd(dJ122L,PDstandard2nd1gzz,kmadd(dJ222L,PDstandard2nd2gzz,kmadd(dJ322L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,SQR(J12L),kmadd(PDstandard2nd22gzz,SQR(J22L),kmadd(PDstandard2nd33gzz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz)),kmul(J22L,kmul(J32L,PDstandard2nd23gzz))),ToReal(2)))))))); JacPDstandard2nd33gxx = kmadd(dJ133L,PDstandard2nd1gxx,kmadd(dJ233L,PDstandard2nd2gxx,kmadd(dJ333L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,SQR(J13L),kmadd(PDstandard2nd22gxx,SQR(J23L),kmadd(PDstandard2nd33gxx,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmul(J23L,kmul(J33L,PDstandard2nd23gxx))),ToReal(2)))))))); JacPDstandard2nd33gxy = kmadd(dJ133L,PDstandard2nd1gxy,kmadd(dJ233L,PDstandard2nd2gxy,kmadd(dJ333L,PDstandard2nd3gxy,kmadd(PDstandard2nd11gxy,SQR(J13L),kmadd(PDstandard2nd22gxy,SQR(J23L),kmadd(PDstandard2nd33gxy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmul(J23L,kmul(J33L,PDstandard2nd23gxy))),ToReal(2)))))))); JacPDstandard2nd33gxz = kmadd(dJ133L,PDstandard2nd1gxz,kmadd(dJ233L,PDstandard2nd2gxz,kmadd(dJ333L,PDstandard2nd3gxz,kmadd(PDstandard2nd11gxz,SQR(J13L),kmadd(PDstandard2nd22gxz,SQR(J23L),kmadd(PDstandard2nd33gxz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxz,kmul(J33L,PDstandard2nd13gxz)),kmul(J23L,kmul(J33L,PDstandard2nd23gxz))),ToReal(2)))))))); JacPDstandard2nd33gyy = kmadd(dJ133L,PDstandard2nd1gyy,kmadd(dJ233L,PDstandard2nd2gyy,kmadd(dJ333L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,SQR(J13L),kmadd(PDstandard2nd22gyy,SQR(J23L),kmadd(PDstandard2nd33gyy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmul(J23L,kmul(J33L,PDstandard2nd23gyy))),ToReal(2)))))))); JacPDstandard2nd33gyz = kmadd(dJ133L,PDstandard2nd1gyz,kmadd(dJ233L,PDstandard2nd2gyz,kmadd(dJ333L,PDstandard2nd3gyz,kmadd(PDstandard2nd11gyz,SQR(J13L),kmadd(PDstandard2nd22gyz,SQR(J23L),kmadd(PDstandard2nd33gyz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gyz,kmul(J33L,PDstandard2nd13gyz)),kmul(J23L,kmul(J33L,PDstandard2nd23gyz))),ToReal(2)))))))); JacPDstandard2nd12gxx = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxx,kmadd(J21L,PDstandard2nd12gxx,kmul(J31L,PDstandard2nd13gxx))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx)),kmadd(dJ112L,PDstandard2nd1gxx,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxx,kmul(J31L,PDstandard2nd23gxx)),kmadd(dJ212L,PDstandard2nd2gxx,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxx,kmul(J31L,PDstandard2nd33gxx)),kmul(dJ312L,PDstandard2nd3gxx))))))); JacPDstandard2nd12gxy = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxy,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy)),kmadd(dJ112L,PDstandard2nd1gxy,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxy,kmul(J31L,PDstandard2nd23gxy)),kmadd(dJ212L,PDstandard2nd2gxy,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxy,kmul(J31L,PDstandard2nd33gxy)),kmul(dJ312L,PDstandard2nd3gxy))))))); JacPDstandard2nd12gxz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxz,kmadd(J21L,PDstandard2nd12gxz,kmul(J31L,PDstandard2nd13gxz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz)),kmadd(dJ112L,PDstandard2nd1gxz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxz,kmul(J31L,PDstandard2nd23gxz)),kmadd(dJ212L,PDstandard2nd2gxz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxz,kmul(J31L,PDstandard2nd33gxz)),kmul(dJ312L,PDstandard2nd3gxz))))))); JacPDstandard2nd12gyy = kmadd(J12L,kmadd(J11L,PDstandard2nd11gyy,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gyy,kmul(J32L,PDstandard2nd13gyy)),kmadd(dJ112L,PDstandard2nd1gyy,kmadd(J22L,kmadd(J21L,PDstandard2nd22gyy,kmul(J31L,PDstandard2nd23gyy)),kmadd(dJ212L,PDstandard2nd2gyy,kmadd(J32L,kmadd(J21L,PDstandard2nd23gyy,kmul(J31L,PDstandard2nd33gyy)),kmul(dJ312L,PDstandard2nd3gyy))))))); JacPDstandard2nd12gyz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gyz,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gyz,kmul(J32L,PDstandard2nd13gyz)),kmadd(dJ112L,PDstandard2nd1gyz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gyz,kmul(J31L,PDstandard2nd23gyz)),kmadd(dJ212L,PDstandard2nd2gyz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gyz,kmul(J31L,PDstandard2nd33gyz)),kmul(dJ312L,PDstandard2nd3gyz))))))); JacPDstandard2nd12gzz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gzz,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz)),kmadd(dJ112L,PDstandard2nd1gzz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gzz,kmul(J31L,PDstandard2nd23gzz)),kmadd(dJ212L,PDstandard2nd2gzz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gzz,kmul(J31L,PDstandard2nd33gzz)),kmul(dJ312L,PDstandard2nd3gzz))))))); JacPDstandard2nd13gxx = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxx,kmadd(J21L,PDstandard2nd12gxx,kmul(J31L,PDstandard2nd13gxx))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmadd(dJ113L,PDstandard2nd1gxx,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxx,kmul(J31L,PDstandard2nd23gxx)),kmadd(dJ213L,PDstandard2nd2gxx,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxx,kmul(J31L,PDstandard2nd33gxx)),kmul(dJ313L,PDstandard2nd3gxx))))))); JacPDstandard2nd13gxy = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxy,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmadd(dJ113L,PDstandard2nd1gxy,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxy,kmul(J31L,PDstandard2nd23gxy)),kmadd(dJ213L,PDstandard2nd2gxy,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxy,kmul(J31L,PDstandard2nd33gxy)),kmul(dJ313L,PDstandard2nd3gxy))))))); JacPDstandard2nd13gxz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxz,kmadd(J21L,PDstandard2nd12gxz,kmul(J31L,PDstandard2nd13gxz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxz,kmul(J33L,PDstandard2nd13gxz)),kmadd(dJ113L,PDstandard2nd1gxz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxz,kmul(J31L,PDstandard2nd23gxz)),kmadd(dJ213L,PDstandard2nd2gxz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxz,kmul(J31L,PDstandard2nd33gxz)),kmul(dJ313L,PDstandard2nd3gxz))))))); JacPDstandard2nd13gyy = kmadd(J13L,kmadd(J11L,PDstandard2nd11gyy,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmadd(dJ113L,PDstandard2nd1gyy,kmadd(J23L,kmadd(J21L,PDstandard2nd22gyy,kmul(J31L,PDstandard2nd23gyy)),kmadd(dJ213L,PDstandard2nd2gyy,kmadd(J33L,kmadd(J21L,PDstandard2nd23gyy,kmul(J31L,PDstandard2nd33gyy)),kmul(dJ313L,PDstandard2nd3gyy))))))); JacPDstandard2nd13gyz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gyz,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gyz,kmul(J33L,PDstandard2nd13gyz)),kmadd(dJ113L,PDstandard2nd1gyz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gyz,kmul(J31L,PDstandard2nd23gyz)),kmadd(dJ213L,PDstandard2nd2gyz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gyz,kmul(J31L,PDstandard2nd33gyz)),kmul(dJ313L,PDstandard2nd3gyz))))))); JacPDstandard2nd13gzz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gzz,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gzz,kmul(J33L,PDstandard2nd13gzz)),kmadd(dJ113L,PDstandard2nd1gzz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gzz,kmul(J31L,PDstandard2nd23gzz)),kmadd(dJ213L,PDstandard2nd2gzz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gzz,kmul(J31L,PDstandard2nd33gzz)),kmul(dJ313L,PDstandard2nd3gzz))))))); JacPDstandard2nd21gxx = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxx,kmadd(J21L,PDstandard2nd12gxx,kmul(J31L,PDstandard2nd13gxx))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx)),kmadd(dJ112L,PDstandard2nd1gxx,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxx,kmul(J31L,PDstandard2nd23gxx)),kmadd(dJ212L,PDstandard2nd2gxx,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxx,kmul(J31L,PDstandard2nd33gxx)),kmul(dJ312L,PDstandard2nd3gxx))))))); JacPDstandard2nd21gxy = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxy,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy)),kmadd(dJ112L,PDstandard2nd1gxy,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxy,kmul(J31L,PDstandard2nd23gxy)),kmadd(dJ212L,PDstandard2nd2gxy,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxy,kmul(J31L,PDstandard2nd33gxy)),kmul(dJ312L,PDstandard2nd3gxy))))))); JacPDstandard2nd21gxz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxz,kmadd(J21L,PDstandard2nd12gxz,kmul(J31L,PDstandard2nd13gxz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz)),kmadd(dJ112L,PDstandard2nd1gxz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxz,kmul(J31L,PDstandard2nd23gxz)),kmadd(dJ212L,PDstandard2nd2gxz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxz,kmul(J31L,PDstandard2nd33gxz)),kmul(dJ312L,PDstandard2nd3gxz))))))); JacPDstandard2nd21gyy = kmadd(J12L,kmadd(J11L,PDstandard2nd11gyy,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gyy,kmul(J32L,PDstandard2nd13gyy)),kmadd(dJ112L,PDstandard2nd1gyy,kmadd(J22L,kmadd(J21L,PDstandard2nd22gyy,kmul(J31L,PDstandard2nd23gyy)),kmadd(dJ212L,PDstandard2nd2gyy,kmadd(J32L,kmadd(J21L,PDstandard2nd23gyy,kmul(J31L,PDstandard2nd33gyy)),kmul(dJ312L,PDstandard2nd3gyy))))))); JacPDstandard2nd21gyz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gyz,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gyz,kmul(J32L,PDstandard2nd13gyz)),kmadd(dJ112L,PDstandard2nd1gyz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gyz,kmul(J31L,PDstandard2nd23gyz)),kmadd(dJ212L,PDstandard2nd2gyz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gyz,kmul(J31L,PDstandard2nd33gyz)),kmul(dJ312L,PDstandard2nd3gyz))))))); JacPDstandard2nd21gzz = kmadd(J12L,kmadd(J11L,PDstandard2nd11gzz,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz)),kmadd(dJ112L,PDstandard2nd1gzz,kmadd(J22L,kmadd(J21L,PDstandard2nd22gzz,kmul(J31L,PDstandard2nd23gzz)),kmadd(dJ212L,PDstandard2nd2gzz,kmadd(J32L,kmadd(J21L,PDstandard2nd23gzz,kmul(J31L,PDstandard2nd33gzz)),kmul(dJ312L,PDstandard2nd3gzz))))))); JacPDstandard2nd23gxx = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxx,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmadd(dJ123L,PDstandard2nd1gxx,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxx,kmul(J32L,PDstandard2nd23gxx)),kmadd(dJ223L,PDstandard2nd2gxx,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxx,kmul(J32L,PDstandard2nd33gxx)),kmul(dJ323L,PDstandard2nd3gxx))))))); JacPDstandard2nd23gxy = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxy,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmadd(dJ123L,PDstandard2nd1gxy,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxy,kmul(J32L,PDstandard2nd23gxy)),kmadd(dJ223L,PDstandard2nd2gxy,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxy,kmul(J32L,PDstandard2nd33gxy)),kmul(dJ323L,PDstandard2nd3gxy))))))); JacPDstandard2nd23gxz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxz,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxz,kmul(J33L,PDstandard2nd13gxz)),kmadd(dJ123L,PDstandard2nd1gxz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxz,kmul(J32L,PDstandard2nd23gxz)),kmadd(dJ223L,PDstandard2nd2gxz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxz,kmul(J32L,PDstandard2nd33gxz)),kmul(dJ323L,PDstandard2nd3gxz))))))); JacPDstandard2nd23gyy = kmadd(J13L,kmadd(J12L,PDstandard2nd11gyy,kmadd(J22L,PDstandard2nd12gyy,kmul(J32L,PDstandard2nd13gyy))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmadd(dJ123L,PDstandard2nd1gyy,kmadd(J23L,kmadd(J22L,PDstandard2nd22gyy,kmul(J32L,PDstandard2nd23gyy)),kmadd(dJ223L,PDstandard2nd2gyy,kmadd(J33L,kmadd(J22L,PDstandard2nd23gyy,kmul(J32L,PDstandard2nd33gyy)),kmul(dJ323L,PDstandard2nd3gyy))))))); JacPDstandard2nd23gyz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gyz,kmadd(J22L,PDstandard2nd12gyz,kmul(J32L,PDstandard2nd13gyz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gyz,kmul(J33L,PDstandard2nd13gyz)),kmadd(dJ123L,PDstandard2nd1gyz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gyz,kmul(J32L,PDstandard2nd23gyz)),kmadd(dJ223L,PDstandard2nd2gyz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gyz,kmul(J32L,PDstandard2nd33gyz)),kmul(dJ323L,PDstandard2nd3gyz))))))); JacPDstandard2nd23gzz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gzz,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gzz,kmul(J33L,PDstandard2nd13gzz)),kmadd(dJ123L,PDstandard2nd1gzz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gzz,kmul(J32L,PDstandard2nd23gzz)),kmadd(dJ223L,PDstandard2nd2gzz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gzz,kmul(J32L,PDstandard2nd33gzz)),kmul(dJ323L,PDstandard2nd3gzz))))))); JacPDstandard2nd31gxx = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxx,kmadd(J21L,PDstandard2nd12gxx,kmul(J31L,PDstandard2nd13gxx))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmadd(dJ113L,PDstandard2nd1gxx,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxx,kmul(J31L,PDstandard2nd23gxx)),kmadd(dJ213L,PDstandard2nd2gxx,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxx,kmul(J31L,PDstandard2nd33gxx)),kmul(dJ313L,PDstandard2nd3gxx))))))); JacPDstandard2nd31gxy = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxy,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmadd(dJ113L,PDstandard2nd1gxy,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxy,kmul(J31L,PDstandard2nd23gxy)),kmadd(dJ213L,PDstandard2nd2gxy,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxy,kmul(J31L,PDstandard2nd33gxy)),kmul(dJ313L,PDstandard2nd3gxy))))))); JacPDstandard2nd31gxz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gxz,kmadd(J21L,PDstandard2nd12gxz,kmul(J31L,PDstandard2nd13gxz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gxz,kmul(J33L,PDstandard2nd13gxz)),kmadd(dJ113L,PDstandard2nd1gxz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gxz,kmul(J31L,PDstandard2nd23gxz)),kmadd(dJ213L,PDstandard2nd2gxz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gxz,kmul(J31L,PDstandard2nd33gxz)),kmul(dJ313L,PDstandard2nd3gxz))))))); JacPDstandard2nd31gyy = kmadd(J13L,kmadd(J11L,PDstandard2nd11gyy,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmadd(dJ113L,PDstandard2nd1gyy,kmadd(J23L,kmadd(J21L,PDstandard2nd22gyy,kmul(J31L,PDstandard2nd23gyy)),kmadd(dJ213L,PDstandard2nd2gyy,kmadd(J33L,kmadd(J21L,PDstandard2nd23gyy,kmul(J31L,PDstandard2nd33gyy)),kmul(dJ313L,PDstandard2nd3gyy))))))); JacPDstandard2nd31gyz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gyz,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gyz,kmul(J33L,PDstandard2nd13gyz)),kmadd(dJ113L,PDstandard2nd1gyz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gyz,kmul(J31L,PDstandard2nd23gyz)),kmadd(dJ213L,PDstandard2nd2gyz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gyz,kmul(J31L,PDstandard2nd33gyz)),kmul(dJ313L,PDstandard2nd3gyz))))))); JacPDstandard2nd31gzz = kmadd(J13L,kmadd(J11L,PDstandard2nd11gzz,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz))),kmadd(J11L,kmadd(J23L,PDstandard2nd12gzz,kmul(J33L,PDstandard2nd13gzz)),kmadd(dJ113L,PDstandard2nd1gzz,kmadd(J23L,kmadd(J21L,PDstandard2nd22gzz,kmul(J31L,PDstandard2nd23gzz)),kmadd(dJ213L,PDstandard2nd2gzz,kmadd(J33L,kmadd(J21L,PDstandard2nd23gzz,kmul(J31L,PDstandard2nd33gzz)),kmul(dJ313L,PDstandard2nd3gzz))))))); JacPDstandard2nd32gxx = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxx,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmadd(dJ123L,PDstandard2nd1gxx,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxx,kmul(J32L,PDstandard2nd23gxx)),kmadd(dJ223L,PDstandard2nd2gxx,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxx,kmul(J32L,PDstandard2nd33gxx)),kmul(dJ323L,PDstandard2nd3gxx))))))); JacPDstandard2nd32gxy = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxy,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmadd(dJ123L,PDstandard2nd1gxy,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxy,kmul(J32L,PDstandard2nd23gxy)),kmadd(dJ223L,PDstandard2nd2gxy,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxy,kmul(J32L,PDstandard2nd33gxy)),kmul(dJ323L,PDstandard2nd3gxy))))))); JacPDstandard2nd32gxz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gxz,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gxz,kmul(J33L,PDstandard2nd13gxz)),kmadd(dJ123L,PDstandard2nd1gxz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gxz,kmul(J32L,PDstandard2nd23gxz)),kmadd(dJ223L,PDstandard2nd2gxz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gxz,kmul(J32L,PDstandard2nd33gxz)),kmul(dJ323L,PDstandard2nd3gxz))))))); JacPDstandard2nd32gyy = kmadd(J13L,kmadd(J12L,PDstandard2nd11gyy,kmadd(J22L,PDstandard2nd12gyy,kmul(J32L,PDstandard2nd13gyy))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmadd(dJ123L,PDstandard2nd1gyy,kmadd(J23L,kmadd(J22L,PDstandard2nd22gyy,kmul(J32L,PDstandard2nd23gyy)),kmadd(dJ223L,PDstandard2nd2gyy,kmadd(J33L,kmadd(J22L,PDstandard2nd23gyy,kmul(J32L,PDstandard2nd33gyy)),kmul(dJ323L,PDstandard2nd3gyy))))))); JacPDstandard2nd32gyz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gyz,kmadd(J22L,PDstandard2nd12gyz,kmul(J32L,PDstandard2nd13gyz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gyz,kmul(J33L,PDstandard2nd13gyz)),kmadd(dJ123L,PDstandard2nd1gyz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gyz,kmul(J32L,PDstandard2nd23gyz)),kmadd(dJ223L,PDstandard2nd2gyz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gyz,kmul(J32L,PDstandard2nd33gyz)),kmul(dJ323L,PDstandard2nd3gyz))))))); JacPDstandard2nd32gzz = kmadd(J13L,kmadd(J12L,PDstandard2nd11gzz,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz))),kmadd(J12L,kmadd(J23L,PDstandard2nd12gzz,kmul(J33L,PDstandard2nd13gzz)),kmadd(dJ123L,PDstandard2nd1gzz,kmadd(J23L,kmadd(J22L,PDstandard2nd22gzz,kmul(J32L,PDstandard2nd23gzz)),kmadd(dJ223L,PDstandard2nd2gzz,kmadd(J33L,kmadd(J22L,PDstandard2nd23gzz,kmul(J32L,PDstandard2nd33gzz)),kmul(dJ323L,PDstandard2nd3gzz))))))); } else { JacPDstandard2nd1gxx = PDstandard2nd1gxx; JacPDstandard2nd1gxy = PDstandard2nd1gxy; JacPDstandard2nd1gxz = PDstandard2nd1gxz; JacPDstandard2nd1gyy = PDstandard2nd1gyy; JacPDstandard2nd1gyz = PDstandard2nd1gyz; JacPDstandard2nd1gzz = PDstandard2nd1gzz; JacPDstandard2nd1kxy = PDstandard2nd1kxy; JacPDstandard2nd1kxz = PDstandard2nd1kxz; JacPDstandard2nd1kyy = PDstandard2nd1kyy; JacPDstandard2nd1kyz = PDstandard2nd1kyz; JacPDstandard2nd1kzz = PDstandard2nd1kzz; JacPDstandard2nd2gxx = PDstandard2nd2gxx; JacPDstandard2nd2gxy = PDstandard2nd2gxy; JacPDstandard2nd2gxz = PDstandard2nd2gxz; JacPDstandard2nd2gyy = PDstandard2nd2gyy; JacPDstandard2nd2gyz = PDstandard2nd2gyz; JacPDstandard2nd2gzz = PDstandard2nd2gzz; JacPDstandard2nd2kxx = PDstandard2nd2kxx; JacPDstandard2nd2kxy = PDstandard2nd2kxy; JacPDstandard2nd2kxz = PDstandard2nd2kxz; JacPDstandard2nd2kyz = PDstandard2nd2kyz; JacPDstandard2nd2kzz = PDstandard2nd2kzz; JacPDstandard2nd3gxx = PDstandard2nd3gxx; JacPDstandard2nd3gxy = PDstandard2nd3gxy; JacPDstandard2nd3gxz = PDstandard2nd3gxz; JacPDstandard2nd3gyy = PDstandard2nd3gyy; JacPDstandard2nd3gyz = PDstandard2nd3gyz; JacPDstandard2nd3gzz = PDstandard2nd3gzz; JacPDstandard2nd3kxx = PDstandard2nd3kxx; JacPDstandard2nd3kxy = PDstandard2nd3kxy; JacPDstandard2nd3kxz = PDstandard2nd3kxz; JacPDstandard2nd3kyy = PDstandard2nd3kyy; JacPDstandard2nd3kyz = PDstandard2nd3kyz; JacPDstandard2nd11gxy = PDstandard2nd11gxy; JacPDstandard2nd11gxz = PDstandard2nd11gxz; JacPDstandard2nd11gyy = PDstandard2nd11gyy; JacPDstandard2nd11gyz = PDstandard2nd11gyz; JacPDstandard2nd11gzz = PDstandard2nd11gzz; JacPDstandard2nd22gxx = PDstandard2nd22gxx; JacPDstandard2nd22gxy = PDstandard2nd22gxy; JacPDstandard2nd22gxz = PDstandard2nd22gxz; JacPDstandard2nd22gyz = PDstandard2nd22gyz; JacPDstandard2nd22gzz = PDstandard2nd22gzz; JacPDstandard2nd33gxx = PDstandard2nd33gxx; JacPDstandard2nd33gxy = PDstandard2nd33gxy; JacPDstandard2nd33gxz = PDstandard2nd33gxz; JacPDstandard2nd33gyy = PDstandard2nd33gyy; JacPDstandard2nd33gyz = PDstandard2nd33gyz; JacPDstandard2nd12gxx = PDstandard2nd12gxx; JacPDstandard2nd12gxy = PDstandard2nd12gxy; JacPDstandard2nd12gxz = PDstandard2nd12gxz; JacPDstandard2nd12gyy = PDstandard2nd12gyy; JacPDstandard2nd12gyz = PDstandard2nd12gyz; JacPDstandard2nd12gzz = PDstandard2nd12gzz; JacPDstandard2nd13gxx = PDstandard2nd13gxx; JacPDstandard2nd13gxy = PDstandard2nd13gxy; JacPDstandard2nd13gxz = PDstandard2nd13gxz; JacPDstandard2nd13gyy = PDstandard2nd13gyy; JacPDstandard2nd13gyz = PDstandard2nd13gyz; JacPDstandard2nd13gzz = PDstandard2nd13gzz; JacPDstandard2nd21gxx = PDstandard2nd12gxx; JacPDstandard2nd21gxy = PDstandard2nd12gxy; JacPDstandard2nd21gxz = PDstandard2nd12gxz; JacPDstandard2nd21gyy = PDstandard2nd12gyy; JacPDstandard2nd21gyz = PDstandard2nd12gyz; JacPDstandard2nd21gzz = PDstandard2nd12gzz; JacPDstandard2nd23gxx = PDstandard2nd23gxx; JacPDstandard2nd23gxy = PDstandard2nd23gxy; JacPDstandard2nd23gxz = PDstandard2nd23gxz; JacPDstandard2nd23gyy = PDstandard2nd23gyy; JacPDstandard2nd23gyz = PDstandard2nd23gyz; JacPDstandard2nd23gzz = PDstandard2nd23gzz; JacPDstandard2nd31gxx = PDstandard2nd13gxx; JacPDstandard2nd31gxy = PDstandard2nd13gxy; JacPDstandard2nd31gxz = PDstandard2nd13gxz; JacPDstandard2nd31gyy = PDstandard2nd13gyy; JacPDstandard2nd31gyz = PDstandard2nd13gyz; JacPDstandard2nd31gzz = PDstandard2nd13gzz; JacPDstandard2nd32gxx = PDstandard2nd23gxx; JacPDstandard2nd32gxy = PDstandard2nd23gxy; JacPDstandard2nd32gxz = PDstandard2nd23gxz; JacPDstandard2nd32gyy = PDstandard2nd23gyy; JacPDstandard2nd32gyz = PDstandard2nd23gyz; JacPDstandard2nd32gzz = PDstandard2nd23gzz; } 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))))))); CCTK_REAL_VEC invdetg = INV(detg); CCTK_REAL_VEC gInv11 = kmul(invdetg,kmsub(gyyL,gzzL,SQR(gyzL))); CCTK_REAL_VEC gInv12 = kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); CCTK_REAL_VEC gInv13 = kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); 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 gInv23 = kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); CCTK_REAL_VEC gInv31 = kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); 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 gamma111 = kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard2nd1gxx,knmsub(gInv12,JacPDstandard2nd2gxx,kmsub(kmadd(gInv12,JacPDstandard2nd1gxy,kmul(gInv13,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv13,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC gamma211 = kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard2nd1gxx,knmsub(gInv22,JacPDstandard2nd2gxx,kmsub(kmadd(gInv22,JacPDstandard2nd1gxy,kmul(gInv23,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv23,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC gamma311 = kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard2nd1gxx,knmsub(gInv32,JacPDstandard2nd2gxx,kmsub(kmadd(gInv32,JacPDstandard2nd1gxy,kmul(gInv33,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv33,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC gamma121 = kmul(kmadd(gInv12,JacPDstandard2nd1gyy,kmadd(gInv11,JacPDstandard2nd2gxx,kmul(gInv13,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC gamma221 = kmul(kmadd(gInv22,JacPDstandard2nd1gyy,kmadd(gInv21,JacPDstandard2nd2gxx,kmul(gInv23,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC gamma321 = kmul(kmadd(gInv32,JacPDstandard2nd1gyy,kmadd(gInv31,JacPDstandard2nd2gxx,kmul(gInv33,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC gamma131 = kmul(kmadd(gInv13,JacPDstandard2nd1gzz,kmadd(gInv11,JacPDstandard2nd3gxx,kmul(gInv12,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC gamma231 = kmul(kmadd(gInv23,JacPDstandard2nd1gzz,kmadd(gInv21,JacPDstandard2nd3gxx,kmul(gInv22,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC gamma331 = kmul(kmadd(gInv33,JacPDstandard2nd1gzz,kmadd(gInv31,JacPDstandard2nd3gxx,kmul(gInv32,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC gamma122 = kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2nd2gyy,kmadd(gInv11,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv13,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC gamma222 = kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2nd2gyy,kmadd(gInv21,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv23,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC gamma322 = kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2nd2gyy,kmadd(gInv31,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv33,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC gamma132 = kmul(kmadd(gInv13,JacPDstandard2nd2gzz,kmadd(gInv12,JacPDstandard2nd3gyy,kmul(gInv11,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC gamma232 = kmul(kmadd(gInv23,JacPDstandard2nd2gzz,kmadd(gInv22,JacPDstandard2nd3gyy,kmul(gInv21,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC gamma332 = kmul(kmadd(gInv33,JacPDstandard2nd2gzz,kmadd(gInv32,JacPDstandard2nd3gyy,kmul(gInv31,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC gamma133 = kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard2nd3gzz,kmadd(gInv11,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv12,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); CCTK_REAL_VEC gamma233 = kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard2nd3gzz,kmadd(gInv21,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv22,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); CCTK_REAL_VEC gamma333 = kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard2nd3gzz,kmadd(gInv31,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv32,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); CCTK_REAL_VEC xmoved = ksub(xL,ToReal(xorig)); CCTK_REAL_VEC ymoved = ksub(yL,ToReal(yorig)); CCTK_REAL_VEC zmoved = ksub(zL,ToReal(zorig)); CCTK_REAL_VEC va1 = kneg(ymoved); CCTK_REAL_VEC va2 = kadd(xmoved,ToReal(offset)); CCTK_REAL_VEC va3 = ToReal(0); CCTK_REAL_VEC vb1 = kadd(xmoved,ToReal(offset)); CCTK_REAL_VEC vb2 = ymoved; CCTK_REAL_VEC vb3 = zmoved; CCTK_REAL_VEC vc1 = kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv11,va2,kmul(gInv12,va1)),kmadd(vb1,kmsub(gInv12,va3,kmul(gInv13,va2)),kmul(vb2,kmsub(gInv13,va1,kmul(gInv11,va3)))))); CCTK_REAL_VEC vc2 = kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv21,va2,kmul(gInv22,va1)),kmadd(vb1,kmsub(gInv22,va3,kmul(gInv23,va2)),kmul(vb2,kmsub(gInv23,va1,kmul(gInv21,va3)))))); CCTK_REAL_VEC vc3 = kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv31,va2,kmul(gInv32,va1)),kmadd(vb1,kmsub(gInv32,va3,kmul(gInv33,va2)),kmul(vb2,kmsub(gInv33,va1,kmul(gInv31,va3)))))); CCTK_REAL_VEC wa1 = va1; CCTK_REAL_VEC wa2 = va2; 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))))); CCTK_REAL_VEC ea1 = kmul(wa1,INV(ksqrt(omega11))); CCTK_REAL_VEC ea2 = kmul(wa2,INV(ksqrt(omega11))); CCTK_REAL_VEC ea3 = kmul(wa3,INV(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)))))); CCTK_REAL_VEC wb1 = knmsub(ea1,omega12,vb1); CCTK_REAL_VEC wb2 = knmsub(ea2,omega12,vb2); 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))))); CCTK_REAL_VEC eb1 = kmul(wb1,INV(ksqrt(omega22))); CCTK_REAL_VEC eb2 = kmul(wb2,INV(ksqrt(omega22))); CCTK_REAL_VEC eb3 = kmul(wb3,INV(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)))))); CCTK_REAL_VEC omega23 = kmadd(eb1,kmadd(gxxL,vc1,kmadd(gxyL,vc2,kmul(gxzL,vc3))),kmadd(eb2,kmadd(gxyL,vc1,kmadd(gyyL,vc2,kmul(gyzL,vc3))),kmul(eb3,kmadd(gxzL,vc1,kmadd(gyzL,vc2,kmul(gzzL,vc3)))))); CCTK_REAL_VEC wc1 = ksub(vc1,kmadd(eb1,omega23,kmul(ea1,omega13))); CCTK_REAL_VEC wc2 = ksub(vc2,kmadd(eb2,omega23,kmul(ea2,omega13))); 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))))); CCTK_REAL_VEC ec1 = kmul(wc1,INV(ksqrt(omega33))); CCTK_REAL_VEC ec2 = kmul(wc2,INV(ksqrt(omega33))); CCTK_REAL_VEC ec3 = kmul(wc3,INV(ksqrt(omega33))); CCTK_REAL_VEC isqrt2 = ToReal(0.707106781186547524); CCTK_REAL_VEC n1 = kneg(kmul(eb1,isqrt2)); CCTK_REAL_VEC n2 = kneg(kmul(eb2,isqrt2)); CCTK_REAL_VEC n3 = kneg(kmul(eb3,isqrt2)); CCTK_REAL_VEC rm1 = kmul(ec1,isqrt2); CCTK_REAL_VEC rm2 = kmul(ec2,isqrt2); CCTK_REAL_VEC rm3 = kmul(ec3,isqrt2); CCTK_REAL_VEC im1 = kmul(ea1,isqrt2); CCTK_REAL_VEC im2 = kmul(ea2,isqrt2); CCTK_REAL_VEC im3 = kmul(ea3,isqrt2); CCTK_REAL_VEC rmbar1 = kmul(ec1,isqrt2); CCTK_REAL_VEC rmbar2 = kmul(ec2,isqrt2); CCTK_REAL_VEC rmbar3 = kmul(ec3,isqrt2); CCTK_REAL_VEC imbar1 = kneg(kmul(ea1,isqrt2)); CCTK_REAL_VEC imbar2 = kneg(kmul(ea2,isqrt2)); CCTK_REAL_VEC imbar3 = kneg(kmul(ea3,isqrt2)); CCTK_REAL_VEC nn = isqrt2; CCTK_REAL_VEC R1111 = ToReal(0); CCTK_REAL_VEC R1112 = kmul(ksub(JacPDstandard2nd21gxx,JacPDstandard2nd12gxx),ToReal(0.5)); CCTK_REAL_VEC R1113 = kmul(ksub(JacPDstandard2nd31gxx,JacPDstandard2nd13gxx),ToReal(0.5)); CCTK_REAL_VEC R1121 = kmul(ksub(JacPDstandard2nd21gxx,JacPDstandard2nd12gxx),ToReal(0.5)); CCTK_REAL_VEC R1122 = ksub(JacPDstandard2nd21gxy,JacPDstandard2nd12gxy); CCTK_REAL_VEC R1123 = kmul(kadd(JacPDstandard2nd21gxz,ksub(JacPDstandard2nd31gxy,kadd(JacPDstandard2nd13gxy,JacPDstandard2nd12gxz))),ToReal(0.5)); CCTK_REAL_VEC R1131 = kmul(ksub(JacPDstandard2nd31gxx,JacPDstandard2nd13gxx),ToReal(0.5)); CCTK_REAL_VEC R1132 = kmul(kadd(JacPDstandard2nd21gxz,ksub(JacPDstandard2nd31gxy,kadd(JacPDstandard2nd13gxy,JacPDstandard2nd12gxz))),ToReal(0.5)); CCTK_REAL_VEC R1133 = ksub(JacPDstandard2nd31gxz,JacPDstandard2nd13gxz); CCTK_REAL_VEC R1211 = kmul(ksub(JacPDstandard2nd12gxx,JacPDstandard2nd21gxx),ToReal(0.5)); CCTK_REAL_VEC R1212 = kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxy,kadd(JacPDstandard2nd21gxy,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))))),JacPDstandard2nd22gxx),JacPDstandard2nd11gyy))))); CCTK_REAL_VEC R1213 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxz,kadd(JacPDstandard2nd31gxy,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),JacPDstandard2nd23gxx),JacPDstandard2nd11gyz)))))); 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(JacPDstandard2nd11gyy,kadd(JacPDstandard2nd22gxx,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),JacPDstandard2nd21gxy),JacPDstandard2nd12gxy)))))),kmul(gxxL,SQR(gamma121))))); CCTK_REAL_VEC R1222 = kmul(ksub(JacPDstandard2nd21gyy,JacPDstandard2nd12gyy),ToReal(0.5)); CCTK_REAL_VEC R1223 = kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd22gxz,kadd(JacPDstandard2nd31gyy,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2),JacPDstandard2nd23gxy),JacPDstandard2nd12gyz)))))); CCTK_REAL_VEC R1231 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd11gyz,kadd(JacPDstandard2nd32gxx,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)))))),JacPDstandard2nd21gxz),JacPDstandard2nd13gxy)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R1232 = kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd21gyz,kadd(JacPDstandard2nd32gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2)),kmul(gamma221,kmul(gamma232,ToReal(2)))))),JacPDstandard2nd22gxz),JacPDstandard2nd13gyy)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R1233 = kmul(kadd(JacPDstandard2nd31gyz,ksub(JacPDstandard2nd32gxz,kadd(JacPDstandard2nd23gxz,JacPDstandard2nd13gyz))),ToReal(0.5)); CCTK_REAL_VEC R1311 = kmul(ksub(JacPDstandard2nd13gxx,JacPDstandard2nd31gxx),ToReal(0.5)); CCTK_REAL_VEC R1312 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxy,kadd(JacPDstandard2nd21gxz,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),JacPDstandard2nd32gxx),JacPDstandard2nd11gyz)))))); CCTK_REAL_VEC R1313 = kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxz,kadd(JacPDstandard2nd31gxz,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))))),JacPDstandard2nd33gxx),JacPDstandard2nd11gzz))))); CCTK_REAL_VEC R1321 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd11gyz,kadd(JacPDstandard2nd23gxx,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)))))),JacPDstandard2nd31gxy),JacPDstandard2nd12gxz)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R1322 = kmul(kadd(JacPDstandard2nd21gyz,ksub(JacPDstandard2nd23gxy,kadd(JacPDstandard2nd32gxy,JacPDstandard2nd12gyz))),ToReal(0.5)); CCTK_REAL_VEC R1323 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd31gyz,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),JacPDstandard2nd33gxy),JacPDstandard2nd12gzz)))))); 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(JacPDstandard2nd11gzz,kadd(JacPDstandard2nd33gxx,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),JacPDstandard2nd31gxz),JacPDstandard2nd13gxz)))))),kmul(gxxL,SQR(gamma131))))); CCTK_REAL_VEC R1332 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd21gzz,kadd(JacPDstandard2nd33gxy,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)))))),JacPDstandard2nd32gxz),JacPDstandard2nd13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); CCTK_REAL_VEC R1333 = kmul(ksub(JacPDstandard2nd31gzz,JacPDstandard2nd13gzz),ToReal(0.5)); CCTK_REAL_VEC R2111 = kmul(ksub(JacPDstandard2nd12gxx,JacPDstandard2nd21gxx),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(JacPDstandard2nd11gyy,kadd(JacPDstandard2nd22gxx,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),JacPDstandard2nd21gxy),JacPDstandard2nd12gxy)))))),kmul(gxxL,SQR(gamma121))))); CCTK_REAL_VEC R2113 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd11gyz,kadd(JacPDstandard2nd32gxx,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)))))),JacPDstandard2nd21gxz),JacPDstandard2nd13gxy)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R2121 = kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxy,kadd(JacPDstandard2nd21gxy,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))))),JacPDstandard2nd22gxx),JacPDstandard2nd11gyy))))); CCTK_REAL_VEC R2122 = kmul(ksub(JacPDstandard2nd21gyy,JacPDstandard2nd12gyy),ToReal(0.5)); CCTK_REAL_VEC R2123 = kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd21gyz,kadd(JacPDstandard2nd32gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2)),kmul(gamma221,kmul(gamma232,ToReal(2)))))),JacPDstandard2nd22gxz),JacPDstandard2nd13gyy)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R2131 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxz,kadd(JacPDstandard2nd31gxy,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),JacPDstandard2nd23gxx),JacPDstandard2nd11gyz)))))); CCTK_REAL_VEC R2132 = kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd22gxz,kadd(JacPDstandard2nd31gyy,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2),JacPDstandard2nd23gxy),JacPDstandard2nd12gyz)))))); CCTK_REAL_VEC R2133 = kmul(kadd(JacPDstandard2nd31gyz,ksub(JacPDstandard2nd32gxz,kadd(JacPDstandard2nd23gxz,JacPDstandard2nd13gyz))),ToReal(0.5)); CCTK_REAL_VEC R2211 = ksub(JacPDstandard2nd12gxy,JacPDstandard2nd21gxy); CCTK_REAL_VEC R2212 = kmul(ksub(JacPDstandard2nd12gyy,JacPDstandard2nd21gyy),ToReal(0.5)); CCTK_REAL_VEC R2213 = kmul(kadd(JacPDstandard2nd12gyz,ksub(JacPDstandard2nd32gxy,kadd(JacPDstandard2nd23gxy,JacPDstandard2nd21gyz))),ToReal(0.5)); CCTK_REAL_VEC R2221 = kmul(ksub(JacPDstandard2nd12gyy,JacPDstandard2nd21gyy),ToReal(0.5)); CCTK_REAL_VEC R2222 = ToReal(0); CCTK_REAL_VEC R2223 = kmul(ksub(JacPDstandard2nd32gyy,JacPDstandard2nd23gyy),ToReal(0.5)); CCTK_REAL_VEC R2231 = kmul(kadd(JacPDstandard2nd12gyz,ksub(JacPDstandard2nd32gxy,kadd(JacPDstandard2nd23gxy,JacPDstandard2nd21gyz))),ToReal(0.5)); CCTK_REAL_VEC R2232 = kmul(ksub(JacPDstandard2nd32gyy,JacPDstandard2nd23gyy),ToReal(0.5)); CCTK_REAL_VEC R2233 = ksub(JacPDstandard2nd32gyz,JacPDstandard2nd23gyz); CCTK_REAL_VEC R2311 = kmul(kadd(JacPDstandard2nd12gxz,ksub(JacPDstandard2nd13gxy,kadd(JacPDstandard2nd31gxy,JacPDstandard2nd21gxz))),ToReal(0.5)); CCTK_REAL_VEC R2312 = kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gyy,kadd(JacPDstandard2nd22gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2),JacPDstandard2nd32gxy),JacPDstandard2nd21gyz)))))); CCTK_REAL_VEC R2313 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gyz,kadd(JacPDstandard2nd32gxz,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),JacPDstandard2nd33gxy),JacPDstandard2nd21gzz)))))); CCTK_REAL_VEC R2321 = kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd12gyz,kadd(JacPDstandard2nd23gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2)),kmul(gamma221,kmul(gamma232,ToReal(2)))))),JacPDstandard2nd31gyy),JacPDstandard2nd22gxz)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R2322 = kmul(ksub(JacPDstandard2nd23gyy,JacPDstandard2nd32gyy),ToReal(0.5)); CCTK_REAL_VEC R2323 = kmul(ToReal(0.5),kadd(JacPDstandard2nd23gyz,kadd(JacPDstandard2nd32gyz,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))))),JacPDstandard2nd33gyy),JacPDstandard2nd22gzz))))); CCTK_REAL_VEC R2331 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd12gzz,kadd(JacPDstandard2nd33gxy,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)))))),JacPDstandard2nd31gyz),JacPDstandard2nd23gxz)))),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(JacPDstandard2nd22gzz,kadd(JacPDstandard2nd33gyy,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),JacPDstandard2nd32gyz),JacPDstandard2nd23gyz)))))),kmul(gxxL,SQR(gamma132))))); CCTK_REAL_VEC R2333 = kmul(ksub(JacPDstandard2nd32gzz,JacPDstandard2nd23gzz),ToReal(0.5)); CCTK_REAL_VEC R3111 = kmul(ksub(JacPDstandard2nd13gxx,JacPDstandard2nd31gxx),ToReal(0.5)); CCTK_REAL_VEC R3112 = kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd11gyz,kadd(JacPDstandard2nd23gxx,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)))))),JacPDstandard2nd31gxy),JacPDstandard2nd12gxz)))),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(JacPDstandard2nd11gzz,kadd(JacPDstandard2nd33gxx,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),JacPDstandard2nd31gxz),JacPDstandard2nd13gxz)))))),kmul(gxxL,SQR(gamma131))))); CCTK_REAL_VEC R3121 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxy,kadd(JacPDstandard2nd21gxz,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),JacPDstandard2nd32gxx),JacPDstandard2nd11gyz)))))); CCTK_REAL_VEC R3122 = kmul(kadd(JacPDstandard2nd21gyz,ksub(JacPDstandard2nd23gxy,kadd(JacPDstandard2nd32gxy,JacPDstandard2nd12gyz))),ToReal(0.5)); CCTK_REAL_VEC R3123 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd21gzz,kadd(JacPDstandard2nd33gxy,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)))))),JacPDstandard2nd32gxz),JacPDstandard2nd13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); CCTK_REAL_VEC R3131 = kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxz,kadd(JacPDstandard2nd31gxz,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))))),JacPDstandard2nd33gxx),JacPDstandard2nd11gzz))))); CCTK_REAL_VEC R3132 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd31gyz,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),JacPDstandard2nd33gxy),JacPDstandard2nd12gzz)))))); CCTK_REAL_VEC R3133 = kmul(ksub(JacPDstandard2nd31gzz,JacPDstandard2nd13gzz),ToReal(0.5)); CCTK_REAL_VEC R3211 = kmul(kadd(JacPDstandard2nd12gxz,ksub(JacPDstandard2nd13gxy,kadd(JacPDstandard2nd31gxy,JacPDstandard2nd21gxz))),ToReal(0.5)); CCTK_REAL_VEC R3212 = kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd12gyz,kadd(JacPDstandard2nd23gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2)),kmul(gamma221,kmul(gamma232,ToReal(2)))))),JacPDstandard2nd31gyy),JacPDstandard2nd22gxz)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); CCTK_REAL_VEC R3213 = kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard2nd12gzz,kadd(JacPDstandard2nd33gxy,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)))))),JacPDstandard2nd31gyz),JacPDstandard2nd23gxz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); CCTK_REAL_VEC R3221 = kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gyy,kadd(JacPDstandard2nd22gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2),JacPDstandard2nd32gxy),JacPDstandard2nd21gyz)))))); CCTK_REAL_VEC R3222 = kmul(ksub(JacPDstandard2nd23gyy,JacPDstandard2nd32gyy),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(JacPDstandard2nd22gzz,kadd(JacPDstandard2nd33gyy,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),JacPDstandard2nd32gyz),JacPDstandard2nd23gyz)))))),kmul(gxxL,SQR(gamma132))))); CCTK_REAL_VEC R3231 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard2nd13gyz,kadd(JacPDstandard2nd32gxz,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),JacPDstandard2nd33gxy),JacPDstandard2nd21gzz)))))); CCTK_REAL_VEC R3232 = kmul(ToReal(0.5),kadd(JacPDstandard2nd23gyz,kadd(JacPDstandard2nd32gyz,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))))),JacPDstandard2nd33gyy),JacPDstandard2nd22gzz))))); CCTK_REAL_VEC R3233 = kmul(ksub(JacPDstandard2nd32gzz,JacPDstandard2nd23gzz),ToReal(0.5)); CCTK_REAL_VEC R3311 = ksub(JacPDstandard2nd13gxz,JacPDstandard2nd31gxz); CCTK_REAL_VEC R3312 = kmul(kadd(JacPDstandard2nd13gyz,ksub(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd32gxz,JacPDstandard2nd31gyz))),ToReal(0.5)); CCTK_REAL_VEC R3313 = kmul(ksub(JacPDstandard2nd13gzz,JacPDstandard2nd31gzz),ToReal(0.5)); CCTK_REAL_VEC R3321 = kmul(kadd(JacPDstandard2nd13gyz,ksub(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd32gxz,JacPDstandard2nd31gyz))),ToReal(0.5)); CCTK_REAL_VEC R3322 = ksub(JacPDstandard2nd23gyz,JacPDstandard2nd32gyz); CCTK_REAL_VEC R3323 = kmul(ksub(JacPDstandard2nd23gzz,JacPDstandard2nd32gzz),ToReal(0.5)); CCTK_REAL_VEC R3331 = kmul(ksub(JacPDstandard2nd13gzz,JacPDstandard2nd31gzz),ToReal(0.5)); CCTK_REAL_VEC R3332 = kmul(ksub(JacPDstandard2nd23gzz,JacPDstandard2nd32gzz),ToReal(0.5)); CCTK_REAL_VEC R3333 = ToReal(0); CCTK_REAL_VEC R4p1111 = R1111; CCTK_REAL_VEC R4p1112 = R1112; CCTK_REAL_VEC R4p1113 = R1113; CCTK_REAL_VEC R4p1121 = R1121; CCTK_REAL_VEC R4p1122 = R1122; CCTK_REAL_VEC R4p1123 = R1123; CCTK_REAL_VEC R4p1131 = R1131; CCTK_REAL_VEC R4p1132 = R1132; CCTK_REAL_VEC R4p1133 = R1133; CCTK_REAL_VEC R4p1211 = R1211; CCTK_REAL_VEC R4p1212 = kmadd(kxxL,kyyL,ksub(R1212,SQR(kxyL))); 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 R4p1222 = R1222; CCTK_REAL_VEC R4p1223 = kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1223)); CCTK_REAL_VEC R4p1231 = kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1231)); CCTK_REAL_VEC R4p1232 = kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R1232)); CCTK_REAL_VEC R4p1233 = R1233; CCTK_REAL_VEC R4p1311 = R1311; 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 R4p1321 = kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1321)); CCTK_REAL_VEC R4p1322 = R1322; 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 R4p1332 = kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1332)); CCTK_REAL_VEC R4p1333 = R1333; CCTK_REAL_VEC R4p2111 = R2111; CCTK_REAL_VEC R4p2112 = kadd(R2112,knmsub(kxxL,kyyL,SQR(kxyL))); 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 R4p2122 = R2122; CCTK_REAL_VEC R4p2123 = kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2123)); CCTK_REAL_VEC R4p2131 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R2131)); CCTK_REAL_VEC R4p2132 = kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2132)); CCTK_REAL_VEC R4p2133 = R2133; CCTK_REAL_VEC R4p2211 = R2211; CCTK_REAL_VEC R4p2212 = R2212; CCTK_REAL_VEC R4p2213 = R2213; CCTK_REAL_VEC R4p2221 = R2221; CCTK_REAL_VEC R4p2222 = R2222; CCTK_REAL_VEC R4p2223 = R2223; CCTK_REAL_VEC R4p2231 = R2231; CCTK_REAL_VEC R4p2232 = R2232; CCTK_REAL_VEC R4p2233 = R2233; CCTK_REAL_VEC R4p2311 = R2311; CCTK_REAL_VEC R4p2312 = kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2312)); CCTK_REAL_VEC R4p2313 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R2313)); CCTK_REAL_VEC R4p2321 = kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2321)); CCTK_REAL_VEC R4p2322 = R2322; CCTK_REAL_VEC R4p2323 = kmadd(kyyL,kzzL,ksub(R2323,SQR(kyzL))); 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 R4p2333 = R2333; CCTK_REAL_VEC R4p3111 = R3111; 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 R4p3121 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)); CCTK_REAL_VEC R4p3122 = R3122; 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 R4p3132 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)); CCTK_REAL_VEC R4p3133 = R3133; CCTK_REAL_VEC R4p3211 = R3211; CCTK_REAL_VEC R4p3212 = kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R3212)); CCTK_REAL_VEC R4p3213 = kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3213)); CCTK_REAL_VEC R4p3221 = kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R3221)); CCTK_REAL_VEC R4p3222 = R3222; CCTK_REAL_VEC R4p3223 = kadd(R3223,knmsub(kyyL,kzzL,SQR(kyzL))); 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 R4p3233 = R3233; CCTK_REAL_VEC R4p3311 = R3311; CCTK_REAL_VEC R4p3312 = R3312; CCTK_REAL_VEC R4p3313 = R3313; CCTK_REAL_VEC R4p3321 = R3321; CCTK_REAL_VEC R4p3322 = R3322; CCTK_REAL_VEC R4p3323 = R3323; CCTK_REAL_VEC R4p3331 = R3331; CCTK_REAL_VEC R4p3332 = R3332; CCTK_REAL_VEC R4p3333 = R3333; CCTK_REAL_VEC Ro111 = ToReal(0); CCTK_REAL_VEC Ro112 = kmadd(kxxL,gamma121,kmadd(kxzL,gamma321,kadd(JacPDstandard2nd1kxy,knmsub(kyyL,gamma211,knmsub(kyzL,gamma311,kmsub(kxyL,ksub(gamma221,gamma111),JacPDstandard2nd2kxx)))))); CCTK_REAL_VEC Ro113 = kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,kadd(JacPDstandard2nd1kxz,knmsub(kyzL,gamma211,knmsub(kzzL,gamma311,kmsub(kxzL,ksub(gamma331,gamma111),JacPDstandard2nd3kxx)))))); CCTK_REAL_VEC Ro121 = kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,kadd(JacPDstandard2nd2kxx,knmsub(kxxL,gamma121,knmsub(kxzL,gamma321,kmsub(kxyL,ksub(gamma111,gamma221),JacPDstandard2nd1kxy)))))); CCTK_REAL_VEC Ro122 = ToReal(0); CCTK_REAL_VEC Ro123 = kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kadd(JacPDstandard2nd2kxz,knmsub(kxzL,gamma121,knmsub(kzzL,gamma321,kmsub(kyzL,ksub(gamma331,gamma221),JacPDstandard2nd3kxy)))))); CCTK_REAL_VEC Ro131 = kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,kadd(JacPDstandard2nd3kxx,knmsub(kxxL,gamma131,knmsub(kxyL,gamma231,kmsub(kxzL,ksub(gamma111,gamma331),JacPDstandard2nd1kxz)))))); CCTK_REAL_VEC Ro132 = kmadd(kxzL,gamma121,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,knmsub(kxyL,gamma131,knmsub(kyyL,gamma231,kmsub(kyzL,ksub(gamma221,gamma331),JacPDstandard2nd2kxz)))))); CCTK_REAL_VEC Ro133 = ToReal(0); CCTK_REAL_VEC Ro211 = ToReal(0); CCTK_REAL_VEC Ro212 = kmadd(kxxL,gamma122,kmadd(kxzL,gamma322,kadd(JacPDstandard2nd1kyy,knmsub(kyyL,gamma221,knmsub(kyzL,gamma321,kmsub(kxyL,ksub(gamma222,gamma121),JacPDstandard2nd2kxy)))))); CCTK_REAL_VEC Ro213 = kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kadd(JacPDstandard2nd1kyz,knmsub(kyzL,gamma221,knmsub(kzzL,gamma321,kmsub(kxzL,ksub(gamma332,gamma121),JacPDstandard2nd3kxy)))))); CCTK_REAL_VEC Ro221 = kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,kadd(JacPDstandard2nd2kxy,knmsub(kxxL,gamma122,knmsub(kxzL,gamma322,kmsub(kxyL,ksub(gamma121,gamma222),JacPDstandard2nd1kyy)))))); CCTK_REAL_VEC Ro222 = ToReal(0); CCTK_REAL_VEC Ro223 = kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,kadd(JacPDstandard2nd2kyz,knmsub(kxzL,gamma122,knmsub(kzzL,gamma322,kmsub(kyzL,ksub(gamma332,gamma222),JacPDstandard2nd3kyy)))))); CCTK_REAL_VEC Ro231 = kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,knmsub(kxxL,gamma132,knmsub(kxyL,gamma232,kmsub(kxzL,ksub(gamma121,gamma332),JacPDstandard2nd1kyz)))))); CCTK_REAL_VEC Ro232 = kmadd(kxzL,gamma122,kmadd(kzzL,gamma322,kadd(JacPDstandard2nd3kyy,knmsub(kxyL,gamma132,knmsub(kyyL,gamma232,kmsub(kyzL,ksub(gamma222,gamma332),JacPDstandard2nd2kyz)))))); CCTK_REAL_VEC Ro233 = ToReal(0); CCTK_REAL_VEC Ro311 = ToReal(0); CCTK_REAL_VEC Ro312 = kmadd(kxxL,gamma132,kmadd(kxzL,gamma332,kadd(JacPDstandard2nd1kyz,knmsub(kyyL,gamma231,knmsub(kyzL,gamma331,kmsub(kxyL,ksub(gamma232,gamma131),JacPDstandard2nd2kxz)))))); CCTK_REAL_VEC Ro313 = kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,kadd(JacPDstandard2nd1kzz,knmsub(kyzL,gamma231,knmsub(kzzL,gamma331,kmsub(kxzL,ksub(gamma333,gamma131),JacPDstandard2nd3kxz)))))); CCTK_REAL_VEC Ro321 = kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kadd(JacPDstandard2nd2kxz,knmsub(kxxL,gamma132,knmsub(kxzL,gamma332,kmsub(kxyL,ksub(gamma131,gamma232),JacPDstandard2nd1kyz)))))); CCTK_REAL_VEC Ro322 = ToReal(0); CCTK_REAL_VEC Ro323 = kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,kadd(JacPDstandard2nd2kzz,knmsub(kxzL,gamma132,knmsub(kzzL,gamma332,kmsub(kyzL,ksub(gamma333,gamma232),JacPDstandard2nd3kyz)))))); CCTK_REAL_VEC Ro331 = kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,kadd(JacPDstandard2nd3kxz,knmsub(kxxL,gamma133,knmsub(kxyL,gamma233,kmsub(kxzL,ksub(gamma131,gamma333),JacPDstandard2nd1kzz)))))); CCTK_REAL_VEC Ro332 = kmadd(kxzL,gamma132,kmadd(kzzL,gamma332,kadd(JacPDstandard2nd3kyz,knmsub(kxyL,gamma133,knmsub(kyyL,gamma233,kmsub(kyzL,ksub(gamma232,gamma333),JacPDstandard2nd2kzz)))))); 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)))))))))))); 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))))))))))); 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))))))))))); 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))))))))))); 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)))))))))))); 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)))))))))))); 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))))))))))); 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)))))))))))); 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)))))))))))); 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)))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))); 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)))))))))))))))))))))))))); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); vec_store_nta_partial(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } LC_ENDLOOP3VEC(psi4_calc_2nd); } extern "C" void psi4_calc_2nd(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; if (verbose > 1) { CCTK_VInfo(CCTK_THORNSTRING,"Entering psi4_calc_2nd_Body"); } if (cctk_iteration % psi4_calc_2nd_calc_every != psi4_calc_2nd_calc_offset) { return; } const char *const groups[] = { "admbase::curv", "admbase::metric", "grid::coordinates", "WeylScal4::Psi4i_group", "WeylScal4::Psi4r_group"}; GenericFD_AssertGroupStorage(cctkGH, "psi4_calc_2nd", 5, groups); switch(fdOrder) { case 2: GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_2nd", 1, 1, 1); break; case 4: GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_2nd", 1, 1, 1); break; case 6: GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_2nd", 1, 1, 1); break; case 8: GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_2nd", 1, 1, 1); break; } GenericFD_LoopOverInterior(cctkGH, psi4_calc_2nd_Body); if (verbose > 1) { CCTK_VInfo(CCTK_THORNSTRING,"Leaving psi4_calc_2nd_Body"); } }