diff options
Diffstat (limited to 'src/psi4_calc_Nth.cc')
| -rw-r--r-- | src/psi4_calc_Nth.cc | 1889 |
1 files changed, 1889 insertions, 0 deletions
diff --git a/src/psi4_calc_Nth.cc b/src/psi4_calc_Nth.cc new file mode 100644 index 0000000..3d795d2 --- /dev/null +++ b/src/psi4_calc_Nth.cc @@ -0,0 +1,1889 @@ +/* File produced by Kranc */ + +#define KRANC_C + +#include <assert.h> +#include <math.h> +#include <stdio.h> +#include <stdlib.h> +#include <string.h> +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "GenericFD.h" +#include "Differencing.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_Nth_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_Nth_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; + + + /* Declare finite differencing variables */ + + /* 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(dx),kmul(INV(dy),ToReal(0.00694444444444444444444444444444))); + CCTK_REAL_VEC const p1o144dxdz = kmul(INV(dx),kmul(INV(dz),ToReal(0.00694444444444444444444444444444))); + CCTK_REAL_VEC const p1o144dydz = kmul(INV(dy),kmul(INV(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(dx),kmul(INV(dy),ToReal(0.000277777777777777777777777777778))); + CCTK_REAL_VEC const p1o3600dxdz = kmul(INV(dx),kmul(INV(dz),ToReal(0.000277777777777777777777777777778))); + CCTK_REAL_VEC const p1o3600dydz = kmul(INV(dy),kmul(INV(dz),ToReal(0.000277777777777777777777777777778))); + CCTK_REAL_VEC const p1o4dxdy = kmul(INV(dx),kmul(INV(dy),ToReal(0.25))); + CCTK_REAL_VEC const p1o4dxdz = kmul(INV(dx),kmul(INV(dz),ToReal(0.25))); + CCTK_REAL_VEC const p1o4dydz = kmul(INV(dy),kmul(INV(dz),ToReal(0.25))); + CCTK_REAL_VEC const p1o5040dx2 = kmul(INV(SQR(dx)),ToReal(0.000198412698412698412698412698413)); + CCTK_REAL_VEC const p1o5040dy2 = kmul(INV(SQR(dy)),ToReal(0.000198412698412698412698412698413)); + CCTK_REAL_VEC const p1o5040dz2 = kmul(INV(SQR(dz)),ToReal(0.000198412698412698412698412698413)); + CCTK_REAL_VEC const 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(dx),kmul(INV(dy),ToReal(1.41723356009070294784580498866e-6))); + CCTK_REAL_VEC const p1o705600dxdz = kmul(INV(dx),kmul(INV(dz),ToReal(1.41723356009070294784580498866e-6))); + CCTK_REAL_VEC const p1o705600dydz = kmul(INV(dy),kmul(INV(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_Nth, + 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 PDstandard1gxx; + CCTK_REAL_VEC PDstandard2gxx; + CCTK_REAL_VEC PDstandard3gxx; + CCTK_REAL_VEC PDstandard11gxx; + CCTK_REAL_VEC PDstandard22gxx; + CCTK_REAL_VEC PDstandard33gxx; + CCTK_REAL_VEC PDstandard12gxx; + CCTK_REAL_VEC PDstandard13gxx; + CCTK_REAL_VEC PDstandard23gxx; + CCTK_REAL_VEC PDstandard1gxy; + CCTK_REAL_VEC PDstandard2gxy; + CCTK_REAL_VEC PDstandard3gxy; + CCTK_REAL_VEC PDstandard11gxy; + CCTK_REAL_VEC PDstandard22gxy; + CCTK_REAL_VEC PDstandard33gxy; + CCTK_REAL_VEC PDstandard12gxy; + CCTK_REAL_VEC PDstandard13gxy; + CCTK_REAL_VEC PDstandard23gxy; + CCTK_REAL_VEC PDstandard1gxz; + CCTK_REAL_VEC PDstandard2gxz; + CCTK_REAL_VEC PDstandard3gxz; + CCTK_REAL_VEC PDstandard11gxz; + CCTK_REAL_VEC PDstandard22gxz; + CCTK_REAL_VEC PDstandard33gxz; + CCTK_REAL_VEC PDstandard12gxz; + CCTK_REAL_VEC PDstandard13gxz; + CCTK_REAL_VEC PDstandard23gxz; + CCTK_REAL_VEC PDstandard1gyy; + CCTK_REAL_VEC PDstandard2gyy; + CCTK_REAL_VEC PDstandard3gyy; + CCTK_REAL_VEC PDstandard11gyy; + CCTK_REAL_VEC PDstandard22gyy; + CCTK_REAL_VEC PDstandard33gyy; + CCTK_REAL_VEC PDstandard12gyy; + CCTK_REAL_VEC PDstandard13gyy; + CCTK_REAL_VEC PDstandard23gyy; + CCTK_REAL_VEC PDstandard1gyz; + CCTK_REAL_VEC PDstandard2gyz; + CCTK_REAL_VEC PDstandard3gyz; + CCTK_REAL_VEC PDstandard11gyz; + CCTK_REAL_VEC PDstandard22gyz; + CCTK_REAL_VEC PDstandard33gyz; + CCTK_REAL_VEC PDstandard12gyz; + CCTK_REAL_VEC PDstandard13gyz; + CCTK_REAL_VEC PDstandard23gyz; + CCTK_REAL_VEC PDstandard1gzz; + CCTK_REAL_VEC PDstandard2gzz; + CCTK_REAL_VEC PDstandard3gzz; + CCTK_REAL_VEC PDstandard11gzz; + CCTK_REAL_VEC PDstandard22gzz; + CCTK_REAL_VEC PDstandard33gzz; + CCTK_REAL_VEC PDstandard12gzz; + CCTK_REAL_VEC PDstandard13gzz; + CCTK_REAL_VEC PDstandard23gzz; + CCTK_REAL_VEC PDstandard1kxx; + CCTK_REAL_VEC PDstandard2kxx; + CCTK_REAL_VEC PDstandard3kxx; + CCTK_REAL_VEC PDstandard1kxy; + CCTK_REAL_VEC PDstandard2kxy; + CCTK_REAL_VEC PDstandard3kxy; + CCTK_REAL_VEC PDstandard1kxz; + CCTK_REAL_VEC PDstandard2kxz; + CCTK_REAL_VEC PDstandard3kxz; + CCTK_REAL_VEC PDstandard1kyy; + CCTK_REAL_VEC PDstandard2kyy; + CCTK_REAL_VEC PDstandard3kyy; + CCTK_REAL_VEC PDstandard1kyz; + CCTK_REAL_VEC PDstandard2kyz; + CCTK_REAL_VEC PDstandard3kyz; + CCTK_REAL_VEC PDstandard1kzz; + CCTK_REAL_VEC PDstandard2kzz; + CCTK_REAL_VEC PDstandard3kzz; + + switch(fdOrder) + { + case 2: + PDstandard1gxx = PDstandardfdOrder21(&gxx[index]); + PDstandard2gxx = PDstandardfdOrder22(&gxx[index]); + PDstandard3gxx = PDstandardfdOrder23(&gxx[index]); + PDstandard11gxx = PDstandardfdOrder211(&gxx[index]); + PDstandard22gxx = PDstandardfdOrder222(&gxx[index]); + PDstandard33gxx = PDstandardfdOrder233(&gxx[index]); + PDstandard12gxx = PDstandardfdOrder212(&gxx[index]); + PDstandard13gxx = PDstandardfdOrder213(&gxx[index]); + PDstandard23gxx = PDstandardfdOrder223(&gxx[index]); + PDstandard1gxy = PDstandardfdOrder21(&gxy[index]); + PDstandard2gxy = PDstandardfdOrder22(&gxy[index]); + PDstandard3gxy = PDstandardfdOrder23(&gxy[index]); + PDstandard11gxy = PDstandardfdOrder211(&gxy[index]); + PDstandard22gxy = PDstandardfdOrder222(&gxy[index]); + PDstandard33gxy = PDstandardfdOrder233(&gxy[index]); + PDstandard12gxy = PDstandardfdOrder212(&gxy[index]); + PDstandard13gxy = PDstandardfdOrder213(&gxy[index]); + PDstandard23gxy = PDstandardfdOrder223(&gxy[index]); + PDstandard1gxz = PDstandardfdOrder21(&gxz[index]); + PDstandard2gxz = PDstandardfdOrder22(&gxz[index]); + PDstandard3gxz = PDstandardfdOrder23(&gxz[index]); + PDstandard11gxz = PDstandardfdOrder211(&gxz[index]); + PDstandard22gxz = PDstandardfdOrder222(&gxz[index]); + PDstandard33gxz = PDstandardfdOrder233(&gxz[index]); + PDstandard12gxz = PDstandardfdOrder212(&gxz[index]); + PDstandard13gxz = PDstandardfdOrder213(&gxz[index]); + PDstandard23gxz = PDstandardfdOrder223(&gxz[index]); + PDstandard1gyy = PDstandardfdOrder21(&gyy[index]); + PDstandard2gyy = PDstandardfdOrder22(&gyy[index]); + PDstandard3gyy = PDstandardfdOrder23(&gyy[index]); + PDstandard11gyy = PDstandardfdOrder211(&gyy[index]); + PDstandard22gyy = PDstandardfdOrder222(&gyy[index]); + PDstandard33gyy = PDstandardfdOrder233(&gyy[index]); + PDstandard12gyy = PDstandardfdOrder212(&gyy[index]); + PDstandard13gyy = PDstandardfdOrder213(&gyy[index]); + PDstandard23gyy = PDstandardfdOrder223(&gyy[index]); + PDstandard1gyz = PDstandardfdOrder21(&gyz[index]); + PDstandard2gyz = PDstandardfdOrder22(&gyz[index]); + PDstandard3gyz = PDstandardfdOrder23(&gyz[index]); + PDstandard11gyz = PDstandardfdOrder211(&gyz[index]); + PDstandard22gyz = PDstandardfdOrder222(&gyz[index]); + PDstandard33gyz = PDstandardfdOrder233(&gyz[index]); + PDstandard12gyz = PDstandardfdOrder212(&gyz[index]); + PDstandard13gyz = PDstandardfdOrder213(&gyz[index]); + PDstandard23gyz = PDstandardfdOrder223(&gyz[index]); + PDstandard1gzz = PDstandardfdOrder21(&gzz[index]); + PDstandard2gzz = PDstandardfdOrder22(&gzz[index]); + PDstandard3gzz = PDstandardfdOrder23(&gzz[index]); + PDstandard11gzz = PDstandardfdOrder211(&gzz[index]); + PDstandard22gzz = PDstandardfdOrder222(&gzz[index]); + PDstandard33gzz = PDstandardfdOrder233(&gzz[index]); + PDstandard12gzz = PDstandardfdOrder212(&gzz[index]); + PDstandard13gzz = PDstandardfdOrder213(&gzz[index]); + PDstandard23gzz = PDstandardfdOrder223(&gzz[index]); + PDstandard1kxx = PDstandardfdOrder21(&kxx[index]); + PDstandard2kxx = PDstandardfdOrder22(&kxx[index]); + PDstandard3kxx = PDstandardfdOrder23(&kxx[index]); + PDstandard1kxy = PDstandardfdOrder21(&kxy[index]); + PDstandard2kxy = PDstandardfdOrder22(&kxy[index]); + PDstandard3kxy = PDstandardfdOrder23(&kxy[index]); + PDstandard1kxz = PDstandardfdOrder21(&kxz[index]); + PDstandard2kxz = PDstandardfdOrder22(&kxz[index]); + PDstandard3kxz = PDstandardfdOrder23(&kxz[index]); + PDstandard1kyy = PDstandardfdOrder21(&kyy[index]); + PDstandard2kyy = PDstandardfdOrder22(&kyy[index]); + PDstandard3kyy = PDstandardfdOrder23(&kyy[index]); + PDstandard1kyz = PDstandardfdOrder21(&kyz[index]); + PDstandard2kyz = PDstandardfdOrder22(&kyz[index]); + PDstandard3kyz = PDstandardfdOrder23(&kyz[index]); + PDstandard1kzz = PDstandardfdOrder21(&kzz[index]); + PDstandard2kzz = PDstandardfdOrder22(&kzz[index]); + PDstandard3kzz = PDstandardfdOrder23(&kzz[index]); + break; + + case 4: + PDstandard1gxx = PDstandardfdOrder41(&gxx[index]); + PDstandard2gxx = PDstandardfdOrder42(&gxx[index]); + PDstandard3gxx = PDstandardfdOrder43(&gxx[index]); + PDstandard11gxx = PDstandardfdOrder411(&gxx[index]); + PDstandard22gxx = PDstandardfdOrder422(&gxx[index]); + PDstandard33gxx = PDstandardfdOrder433(&gxx[index]); + PDstandard12gxx = PDstandardfdOrder412(&gxx[index]); + PDstandard13gxx = PDstandardfdOrder413(&gxx[index]); + PDstandard23gxx = PDstandardfdOrder423(&gxx[index]); + PDstandard1gxy = PDstandardfdOrder41(&gxy[index]); + PDstandard2gxy = PDstandardfdOrder42(&gxy[index]); + PDstandard3gxy = PDstandardfdOrder43(&gxy[index]); + PDstandard11gxy = PDstandardfdOrder411(&gxy[index]); + PDstandard22gxy = PDstandardfdOrder422(&gxy[index]); + PDstandard33gxy = PDstandardfdOrder433(&gxy[index]); + PDstandard12gxy = PDstandardfdOrder412(&gxy[index]); + PDstandard13gxy = PDstandardfdOrder413(&gxy[index]); + PDstandard23gxy = PDstandardfdOrder423(&gxy[index]); + PDstandard1gxz = PDstandardfdOrder41(&gxz[index]); + PDstandard2gxz = PDstandardfdOrder42(&gxz[index]); + PDstandard3gxz = PDstandardfdOrder43(&gxz[index]); + PDstandard11gxz = PDstandardfdOrder411(&gxz[index]); + PDstandard22gxz = PDstandardfdOrder422(&gxz[index]); + PDstandard33gxz = PDstandardfdOrder433(&gxz[index]); + PDstandard12gxz = PDstandardfdOrder412(&gxz[index]); + PDstandard13gxz = PDstandardfdOrder413(&gxz[index]); + PDstandard23gxz = PDstandardfdOrder423(&gxz[index]); + PDstandard1gyy = PDstandardfdOrder41(&gyy[index]); + PDstandard2gyy = PDstandardfdOrder42(&gyy[index]); + PDstandard3gyy = PDstandardfdOrder43(&gyy[index]); + PDstandard11gyy = PDstandardfdOrder411(&gyy[index]); + PDstandard22gyy = PDstandardfdOrder422(&gyy[index]); + PDstandard33gyy = PDstandardfdOrder433(&gyy[index]); + PDstandard12gyy = PDstandardfdOrder412(&gyy[index]); + PDstandard13gyy = PDstandardfdOrder413(&gyy[index]); + PDstandard23gyy = PDstandardfdOrder423(&gyy[index]); + PDstandard1gyz = PDstandardfdOrder41(&gyz[index]); + PDstandard2gyz = PDstandardfdOrder42(&gyz[index]); + PDstandard3gyz = PDstandardfdOrder43(&gyz[index]); + PDstandard11gyz = PDstandardfdOrder411(&gyz[index]); + PDstandard22gyz = PDstandardfdOrder422(&gyz[index]); + PDstandard33gyz = PDstandardfdOrder433(&gyz[index]); + PDstandard12gyz = PDstandardfdOrder412(&gyz[index]); + PDstandard13gyz = PDstandardfdOrder413(&gyz[index]); + PDstandard23gyz = PDstandardfdOrder423(&gyz[index]); + PDstandard1gzz = PDstandardfdOrder41(&gzz[index]); + PDstandard2gzz = PDstandardfdOrder42(&gzz[index]); + PDstandard3gzz = PDstandardfdOrder43(&gzz[index]); + PDstandard11gzz = PDstandardfdOrder411(&gzz[index]); + PDstandard22gzz = PDstandardfdOrder422(&gzz[index]); + PDstandard33gzz = PDstandardfdOrder433(&gzz[index]); + PDstandard12gzz = PDstandardfdOrder412(&gzz[index]); + PDstandard13gzz = PDstandardfdOrder413(&gzz[index]); + PDstandard23gzz = PDstandardfdOrder423(&gzz[index]); + PDstandard1kxx = PDstandardfdOrder41(&kxx[index]); + PDstandard2kxx = PDstandardfdOrder42(&kxx[index]); + PDstandard3kxx = PDstandardfdOrder43(&kxx[index]); + PDstandard1kxy = PDstandardfdOrder41(&kxy[index]); + PDstandard2kxy = PDstandardfdOrder42(&kxy[index]); + PDstandard3kxy = PDstandardfdOrder43(&kxy[index]); + PDstandard1kxz = PDstandardfdOrder41(&kxz[index]); + PDstandard2kxz = PDstandardfdOrder42(&kxz[index]); + PDstandard3kxz = PDstandardfdOrder43(&kxz[index]); + PDstandard1kyy = PDstandardfdOrder41(&kyy[index]); + PDstandard2kyy = PDstandardfdOrder42(&kyy[index]); + PDstandard3kyy = PDstandardfdOrder43(&kyy[index]); + PDstandard1kyz = PDstandardfdOrder41(&kyz[index]); + PDstandard2kyz = PDstandardfdOrder42(&kyz[index]); + PDstandard3kyz = PDstandardfdOrder43(&kyz[index]); + PDstandard1kzz = PDstandardfdOrder41(&kzz[index]); + PDstandard2kzz = PDstandardfdOrder42(&kzz[index]); + PDstandard3kzz = PDstandardfdOrder43(&kzz[index]); + break; + + case 6: + PDstandard1gxx = PDstandardfdOrder61(&gxx[index]); + PDstandard2gxx = PDstandardfdOrder62(&gxx[index]); + PDstandard3gxx = PDstandardfdOrder63(&gxx[index]); + PDstandard11gxx = PDstandardfdOrder611(&gxx[index]); + PDstandard22gxx = PDstandardfdOrder622(&gxx[index]); + PDstandard33gxx = PDstandardfdOrder633(&gxx[index]); + PDstandard12gxx = PDstandardfdOrder612(&gxx[index]); + PDstandard13gxx = PDstandardfdOrder613(&gxx[index]); + PDstandard23gxx = PDstandardfdOrder623(&gxx[index]); + PDstandard1gxy = PDstandardfdOrder61(&gxy[index]); + PDstandard2gxy = PDstandardfdOrder62(&gxy[index]); + PDstandard3gxy = PDstandardfdOrder63(&gxy[index]); + PDstandard11gxy = PDstandardfdOrder611(&gxy[index]); + PDstandard22gxy = PDstandardfdOrder622(&gxy[index]); + PDstandard33gxy = PDstandardfdOrder633(&gxy[index]); + PDstandard12gxy = PDstandardfdOrder612(&gxy[index]); + PDstandard13gxy = PDstandardfdOrder613(&gxy[index]); + PDstandard23gxy = PDstandardfdOrder623(&gxy[index]); + PDstandard1gxz = PDstandardfdOrder61(&gxz[index]); + PDstandard2gxz = PDstandardfdOrder62(&gxz[index]); + PDstandard3gxz = PDstandardfdOrder63(&gxz[index]); + PDstandard11gxz = PDstandardfdOrder611(&gxz[index]); + PDstandard22gxz = PDstandardfdOrder622(&gxz[index]); + PDstandard33gxz = PDstandardfdOrder633(&gxz[index]); + PDstandard12gxz = PDstandardfdOrder612(&gxz[index]); + PDstandard13gxz = PDstandardfdOrder613(&gxz[index]); + PDstandard23gxz = PDstandardfdOrder623(&gxz[index]); + PDstandard1gyy = PDstandardfdOrder61(&gyy[index]); + PDstandard2gyy = PDstandardfdOrder62(&gyy[index]); + PDstandard3gyy = PDstandardfdOrder63(&gyy[index]); + PDstandard11gyy = PDstandardfdOrder611(&gyy[index]); + PDstandard22gyy = PDstandardfdOrder622(&gyy[index]); + PDstandard33gyy = PDstandardfdOrder633(&gyy[index]); + PDstandard12gyy = PDstandardfdOrder612(&gyy[index]); + PDstandard13gyy = PDstandardfdOrder613(&gyy[index]); + PDstandard23gyy = PDstandardfdOrder623(&gyy[index]); + PDstandard1gyz = PDstandardfdOrder61(&gyz[index]); + PDstandard2gyz = PDstandardfdOrder62(&gyz[index]); + PDstandard3gyz = PDstandardfdOrder63(&gyz[index]); + PDstandard11gyz = PDstandardfdOrder611(&gyz[index]); + PDstandard22gyz = PDstandardfdOrder622(&gyz[index]); + PDstandard33gyz = PDstandardfdOrder633(&gyz[index]); + PDstandard12gyz = PDstandardfdOrder612(&gyz[index]); + PDstandard13gyz = PDstandardfdOrder613(&gyz[index]); + PDstandard23gyz = PDstandardfdOrder623(&gyz[index]); + PDstandard1gzz = PDstandardfdOrder61(&gzz[index]); + PDstandard2gzz = PDstandardfdOrder62(&gzz[index]); + PDstandard3gzz = PDstandardfdOrder63(&gzz[index]); + PDstandard11gzz = PDstandardfdOrder611(&gzz[index]); + PDstandard22gzz = PDstandardfdOrder622(&gzz[index]); + PDstandard33gzz = PDstandardfdOrder633(&gzz[index]); + PDstandard12gzz = PDstandardfdOrder612(&gzz[index]); + PDstandard13gzz = PDstandardfdOrder613(&gzz[index]); + PDstandard23gzz = PDstandardfdOrder623(&gzz[index]); + PDstandard1kxx = PDstandardfdOrder61(&kxx[index]); + PDstandard2kxx = PDstandardfdOrder62(&kxx[index]); + PDstandard3kxx = PDstandardfdOrder63(&kxx[index]); + PDstandard1kxy = PDstandardfdOrder61(&kxy[index]); + PDstandard2kxy = PDstandardfdOrder62(&kxy[index]); + PDstandard3kxy = PDstandardfdOrder63(&kxy[index]); + PDstandard1kxz = PDstandardfdOrder61(&kxz[index]); + PDstandard2kxz = PDstandardfdOrder62(&kxz[index]); + PDstandard3kxz = PDstandardfdOrder63(&kxz[index]); + PDstandard1kyy = PDstandardfdOrder61(&kyy[index]); + PDstandard2kyy = PDstandardfdOrder62(&kyy[index]); + PDstandard3kyy = PDstandardfdOrder63(&kyy[index]); + PDstandard1kyz = PDstandardfdOrder61(&kyz[index]); + PDstandard2kyz = PDstandardfdOrder62(&kyz[index]); + PDstandard3kyz = PDstandardfdOrder63(&kyz[index]); + PDstandard1kzz = PDstandardfdOrder61(&kzz[index]); + PDstandard2kzz = PDstandardfdOrder62(&kzz[index]); + PDstandard3kzz = PDstandardfdOrder63(&kzz[index]); + break; + + case 8: + PDstandard1gxx = PDstandardfdOrder81(&gxx[index]); + PDstandard2gxx = PDstandardfdOrder82(&gxx[index]); + PDstandard3gxx = PDstandardfdOrder83(&gxx[index]); + PDstandard11gxx = PDstandardfdOrder811(&gxx[index]); + PDstandard22gxx = PDstandardfdOrder822(&gxx[index]); + PDstandard33gxx = PDstandardfdOrder833(&gxx[index]); + PDstandard12gxx = PDstandardfdOrder812(&gxx[index]); + PDstandard13gxx = PDstandardfdOrder813(&gxx[index]); + PDstandard23gxx = PDstandardfdOrder823(&gxx[index]); + PDstandard1gxy = PDstandardfdOrder81(&gxy[index]); + PDstandard2gxy = PDstandardfdOrder82(&gxy[index]); + PDstandard3gxy = PDstandardfdOrder83(&gxy[index]); + PDstandard11gxy = PDstandardfdOrder811(&gxy[index]); + PDstandard22gxy = PDstandardfdOrder822(&gxy[index]); + PDstandard33gxy = PDstandardfdOrder833(&gxy[index]); + PDstandard12gxy = PDstandardfdOrder812(&gxy[index]); + PDstandard13gxy = PDstandardfdOrder813(&gxy[index]); + PDstandard23gxy = PDstandardfdOrder823(&gxy[index]); + PDstandard1gxz = PDstandardfdOrder81(&gxz[index]); + PDstandard2gxz = PDstandardfdOrder82(&gxz[index]); + PDstandard3gxz = PDstandardfdOrder83(&gxz[index]); + PDstandard11gxz = PDstandardfdOrder811(&gxz[index]); + PDstandard22gxz = PDstandardfdOrder822(&gxz[index]); + PDstandard33gxz = PDstandardfdOrder833(&gxz[index]); + PDstandard12gxz = PDstandardfdOrder812(&gxz[index]); + PDstandard13gxz = PDstandardfdOrder813(&gxz[index]); + PDstandard23gxz = PDstandardfdOrder823(&gxz[index]); + PDstandard1gyy = PDstandardfdOrder81(&gyy[index]); + PDstandard2gyy = PDstandardfdOrder82(&gyy[index]); + PDstandard3gyy = PDstandardfdOrder83(&gyy[index]); + PDstandard11gyy = PDstandardfdOrder811(&gyy[index]); + PDstandard22gyy = PDstandardfdOrder822(&gyy[index]); + PDstandard33gyy = PDstandardfdOrder833(&gyy[index]); + PDstandard12gyy = PDstandardfdOrder812(&gyy[index]); + PDstandard13gyy = PDstandardfdOrder813(&gyy[index]); + PDstandard23gyy = PDstandardfdOrder823(&gyy[index]); + PDstandard1gyz = PDstandardfdOrder81(&gyz[index]); + PDstandard2gyz = PDstandardfdOrder82(&gyz[index]); + PDstandard3gyz = PDstandardfdOrder83(&gyz[index]); + PDstandard11gyz = PDstandardfdOrder811(&gyz[index]); + PDstandard22gyz = PDstandardfdOrder822(&gyz[index]); + PDstandard33gyz = PDstandardfdOrder833(&gyz[index]); + PDstandard12gyz = PDstandardfdOrder812(&gyz[index]); + PDstandard13gyz = PDstandardfdOrder813(&gyz[index]); + PDstandard23gyz = PDstandardfdOrder823(&gyz[index]); + PDstandard1gzz = PDstandardfdOrder81(&gzz[index]); + PDstandard2gzz = PDstandardfdOrder82(&gzz[index]); + PDstandard3gzz = PDstandardfdOrder83(&gzz[index]); + PDstandard11gzz = PDstandardfdOrder811(&gzz[index]); + PDstandard22gzz = PDstandardfdOrder822(&gzz[index]); + PDstandard33gzz = PDstandardfdOrder833(&gzz[index]); + PDstandard12gzz = PDstandardfdOrder812(&gzz[index]); + PDstandard13gzz = PDstandardfdOrder813(&gzz[index]); + PDstandard23gzz = PDstandardfdOrder823(&gzz[index]); + PDstandard1kxx = PDstandardfdOrder81(&kxx[index]); + PDstandard2kxx = PDstandardfdOrder82(&kxx[index]); + PDstandard3kxx = PDstandardfdOrder83(&kxx[index]); + PDstandard1kxy = PDstandardfdOrder81(&kxy[index]); + PDstandard2kxy = PDstandardfdOrder82(&kxy[index]); + PDstandard3kxy = PDstandardfdOrder83(&kxy[index]); + PDstandard1kxz = PDstandardfdOrder81(&kxz[index]); + PDstandard2kxz = PDstandardfdOrder82(&kxz[index]); + PDstandard3kxz = PDstandardfdOrder83(&kxz[index]); + PDstandard1kyy = PDstandardfdOrder81(&kyy[index]); + PDstandard2kyy = PDstandardfdOrder82(&kyy[index]); + PDstandard3kyy = PDstandardfdOrder83(&kyy[index]); + PDstandard1kyz = PDstandardfdOrder81(&kyz[index]); + PDstandard2kyz = PDstandardfdOrder82(&kyz[index]); + PDstandard3kyz = PDstandardfdOrder83(&kyz[index]); + PDstandard1kzz = PDstandardfdOrder81(&kzz[index]); + PDstandard2kzz = PDstandardfdOrder82(&kzz[index]); + PDstandard3kzz = PDstandardfdOrder83(&kzz[index]); + break; + } + + /* Calculate temporaries and grid functions */ + CCTK_REAL_VEC JacPDstandard11gxy; + CCTK_REAL_VEC JacPDstandard11gxz; + CCTK_REAL_VEC JacPDstandard11gyy; + CCTK_REAL_VEC JacPDstandard11gyz; + CCTK_REAL_VEC JacPDstandard11gzz; + CCTK_REAL_VEC JacPDstandard12gxx; + CCTK_REAL_VEC JacPDstandard12gxy; + CCTK_REAL_VEC JacPDstandard12gxz; + CCTK_REAL_VEC JacPDstandard12gyy; + CCTK_REAL_VEC JacPDstandard12gyz; + CCTK_REAL_VEC JacPDstandard12gzz; + CCTK_REAL_VEC JacPDstandard13gxx; + CCTK_REAL_VEC JacPDstandard13gxy; + CCTK_REAL_VEC JacPDstandard13gxz; + CCTK_REAL_VEC JacPDstandard13gyy; + CCTK_REAL_VEC JacPDstandard13gyz; + CCTK_REAL_VEC JacPDstandard13gzz; + CCTK_REAL_VEC JacPDstandard1gxx; + CCTK_REAL_VEC JacPDstandard1gxy; + CCTK_REAL_VEC JacPDstandard1gxz; + CCTK_REAL_VEC JacPDstandard1gyy; + CCTK_REAL_VEC JacPDstandard1gyz; + CCTK_REAL_VEC JacPDstandard1gzz; + CCTK_REAL_VEC JacPDstandard1kxy; + CCTK_REAL_VEC JacPDstandard1kxz; + CCTK_REAL_VEC JacPDstandard1kyy; + CCTK_REAL_VEC JacPDstandard1kyz; + CCTK_REAL_VEC JacPDstandard1kzz; + CCTK_REAL_VEC JacPDstandard21gxx; + CCTK_REAL_VEC JacPDstandard21gxy; + CCTK_REAL_VEC JacPDstandard21gxz; + CCTK_REAL_VEC JacPDstandard21gyy; + CCTK_REAL_VEC JacPDstandard21gyz; + CCTK_REAL_VEC JacPDstandard21gzz; + CCTK_REAL_VEC JacPDstandard22gxx; + CCTK_REAL_VEC JacPDstandard22gxy; + CCTK_REAL_VEC JacPDstandard22gxz; + CCTK_REAL_VEC JacPDstandard22gyz; + CCTK_REAL_VEC JacPDstandard22gzz; + CCTK_REAL_VEC JacPDstandard23gxx; + CCTK_REAL_VEC JacPDstandard23gxy; + CCTK_REAL_VEC JacPDstandard23gxz; + CCTK_REAL_VEC JacPDstandard23gyy; + CCTK_REAL_VEC JacPDstandard23gyz; + CCTK_REAL_VEC JacPDstandard23gzz; + CCTK_REAL_VEC JacPDstandard2gxx; + CCTK_REAL_VEC JacPDstandard2gxy; + CCTK_REAL_VEC JacPDstandard2gxz; + CCTK_REAL_VEC JacPDstandard2gyy; + CCTK_REAL_VEC JacPDstandard2gyz; + CCTK_REAL_VEC JacPDstandard2gzz; + CCTK_REAL_VEC JacPDstandard2kxx; + CCTK_REAL_VEC JacPDstandard2kxy; + CCTK_REAL_VEC JacPDstandard2kxz; + CCTK_REAL_VEC JacPDstandard2kyz; + CCTK_REAL_VEC JacPDstandard2kzz; + CCTK_REAL_VEC JacPDstandard31gxx; + CCTK_REAL_VEC JacPDstandard31gxy; + CCTK_REAL_VEC JacPDstandard31gxz; + CCTK_REAL_VEC JacPDstandard31gyy; + CCTK_REAL_VEC JacPDstandard31gyz; + CCTK_REAL_VEC JacPDstandard31gzz; + CCTK_REAL_VEC JacPDstandard32gxx; + CCTK_REAL_VEC JacPDstandard32gxy; + CCTK_REAL_VEC JacPDstandard32gxz; + CCTK_REAL_VEC JacPDstandard32gyy; + CCTK_REAL_VEC JacPDstandard32gyz; + CCTK_REAL_VEC JacPDstandard32gzz; + CCTK_REAL_VEC JacPDstandard33gxx; + CCTK_REAL_VEC JacPDstandard33gxy; + CCTK_REAL_VEC JacPDstandard33gxz; + CCTK_REAL_VEC JacPDstandard33gyy; + CCTK_REAL_VEC JacPDstandard33gyz; + CCTK_REAL_VEC JacPDstandard3gxx; + CCTK_REAL_VEC JacPDstandard3gxy; + CCTK_REAL_VEC JacPDstandard3gxz; + CCTK_REAL_VEC JacPDstandard3gyy; + CCTK_REAL_VEC JacPDstandard3gyz; + CCTK_REAL_VEC JacPDstandard3gzz; + CCTK_REAL_VEC JacPDstandard3kxx; + CCTK_REAL_VEC JacPDstandard3kxy; + CCTK_REAL_VEC JacPDstandard3kxz; + CCTK_REAL_VEC JacPDstandard3kyy; + CCTK_REAL_VEC JacPDstandard3kyz; + + if (use_jacobian) + { + JacPDstandard1gxx = + kmadd(J11L,PDstandard1gxx,kmadd(J21L,PDstandard2gxx,kmul(J31L,PDstandard3gxx))); + + JacPDstandard1gxy = + kmadd(J11L,PDstandard1gxy,kmadd(J21L,PDstandard2gxy,kmul(J31L,PDstandard3gxy))); + + JacPDstandard1gxz = + kmadd(J11L,PDstandard1gxz,kmadd(J21L,PDstandard2gxz,kmul(J31L,PDstandard3gxz))); + + JacPDstandard1gyy = + kmadd(J11L,PDstandard1gyy,kmadd(J21L,PDstandard2gyy,kmul(J31L,PDstandard3gyy))); + + JacPDstandard1gyz = + kmadd(J11L,PDstandard1gyz,kmadd(J21L,PDstandard2gyz,kmul(J31L,PDstandard3gyz))); + + JacPDstandard1gzz = + kmadd(J11L,PDstandard1gzz,kmadd(J21L,PDstandard2gzz,kmul(J31L,PDstandard3gzz))); + + JacPDstandard1kxy = + kmadd(J11L,PDstandard1kxy,kmadd(J21L,PDstandard2kxy,kmul(J31L,PDstandard3kxy))); + + JacPDstandard1kxz = + kmadd(J11L,PDstandard1kxz,kmadd(J21L,PDstandard2kxz,kmul(J31L,PDstandard3kxz))); + + JacPDstandard1kyy = + kmadd(J11L,PDstandard1kyy,kmadd(J21L,PDstandard2kyy,kmul(J31L,PDstandard3kyy))); + + JacPDstandard1kyz = + kmadd(J11L,PDstandard1kyz,kmadd(J21L,PDstandard2kyz,kmul(J31L,PDstandard3kyz))); + + JacPDstandard1kzz = + kmadd(J11L,PDstandard1kzz,kmadd(J21L,PDstandard2kzz,kmul(J31L,PDstandard3kzz))); + + JacPDstandard2gxx = + kmadd(J12L,PDstandard1gxx,kmadd(J22L,PDstandard2gxx,kmul(J32L,PDstandard3gxx))); + + JacPDstandard2gxy = + kmadd(J12L,PDstandard1gxy,kmadd(J22L,PDstandard2gxy,kmul(J32L,PDstandard3gxy))); + + JacPDstandard2gxz = + kmadd(J12L,PDstandard1gxz,kmadd(J22L,PDstandard2gxz,kmul(J32L,PDstandard3gxz))); + + JacPDstandard2gyy = + kmadd(J12L,PDstandard1gyy,kmadd(J22L,PDstandard2gyy,kmul(J32L,PDstandard3gyy))); + + JacPDstandard2gyz = + kmadd(J12L,PDstandard1gyz,kmadd(J22L,PDstandard2gyz,kmul(J32L,PDstandard3gyz))); + + JacPDstandard2gzz = + kmadd(J12L,PDstandard1gzz,kmadd(J22L,PDstandard2gzz,kmul(J32L,PDstandard3gzz))); + + JacPDstandard2kxx = + kmadd(J12L,PDstandard1kxx,kmadd(J22L,PDstandard2kxx,kmul(J32L,PDstandard3kxx))); + + JacPDstandard2kxy = + kmadd(J12L,PDstandard1kxy,kmadd(J22L,PDstandard2kxy,kmul(J32L,PDstandard3kxy))); + + JacPDstandard2kxz = + kmadd(J12L,PDstandard1kxz,kmadd(J22L,PDstandard2kxz,kmul(J32L,PDstandard3kxz))); + + JacPDstandard2kyz = + kmadd(J12L,PDstandard1kyz,kmadd(J22L,PDstandard2kyz,kmul(J32L,PDstandard3kyz))); + + JacPDstandard2kzz = + kmadd(J12L,PDstandard1kzz,kmadd(J22L,PDstandard2kzz,kmul(J32L,PDstandard3kzz))); + + JacPDstandard3gxx = + kmadd(J13L,PDstandard1gxx,kmadd(J23L,PDstandard2gxx,kmul(J33L,PDstandard3gxx))); + + JacPDstandard3gxy = + kmadd(J13L,PDstandard1gxy,kmadd(J23L,PDstandard2gxy,kmul(J33L,PDstandard3gxy))); + + JacPDstandard3gxz = + kmadd(J13L,PDstandard1gxz,kmadd(J23L,PDstandard2gxz,kmul(J33L,PDstandard3gxz))); + + JacPDstandard3gyy = + kmadd(J13L,PDstandard1gyy,kmadd(J23L,PDstandard2gyy,kmul(J33L,PDstandard3gyy))); + + JacPDstandard3gyz = + kmadd(J13L,PDstandard1gyz,kmadd(J23L,PDstandard2gyz,kmul(J33L,PDstandard3gyz))); + + JacPDstandard3gzz = + kmadd(J13L,PDstandard1gzz,kmadd(J23L,PDstandard2gzz,kmul(J33L,PDstandard3gzz))); + + JacPDstandard3kxx = + kmadd(J13L,PDstandard1kxx,kmadd(J23L,PDstandard2kxx,kmul(J33L,PDstandard3kxx))); + + JacPDstandard3kxy = + kmadd(J13L,PDstandard1kxy,kmadd(J23L,PDstandard2kxy,kmul(J33L,PDstandard3kxy))); + + JacPDstandard3kxz = + kmadd(J13L,PDstandard1kxz,kmadd(J23L,PDstandard2kxz,kmul(J33L,PDstandard3kxz))); + + JacPDstandard3kyy = + kmadd(J13L,PDstandard1kyy,kmadd(J23L,PDstandard2kyy,kmul(J33L,PDstandard3kyy))); + + JacPDstandard3kyz = + kmadd(J13L,PDstandard1kyz,kmadd(J23L,PDstandard2kyz,kmul(J33L,PDstandard3kyz))); + + JacPDstandard11gxy = + kmadd(dJ111L,PDstandard1gxy,kmadd(dJ211L,PDstandard2gxy,kmadd(dJ311L,PDstandard3gxy,kmadd(PDstandard11gxy,SQR(J11L),kmadd(PDstandard22gxy,SQR(J21L),kmadd(PDstandard33gxy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gxy,kmul(J31L,PDstandard13gxy)),kmul(J21L,kmul(J31L,PDstandard23gxy))),ToReal(2)))))))); + + JacPDstandard11gxz = + kmadd(dJ111L,PDstandard1gxz,kmadd(dJ211L,PDstandard2gxz,kmadd(dJ311L,PDstandard3gxz,kmadd(PDstandard11gxz,SQR(J11L),kmadd(PDstandard22gxz,SQR(J21L),kmadd(PDstandard33gxz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gxz,kmul(J31L,PDstandard13gxz)),kmul(J21L,kmul(J31L,PDstandard23gxz))),ToReal(2)))))))); + + JacPDstandard11gyy = + kmadd(dJ111L,PDstandard1gyy,kmadd(dJ211L,PDstandard2gyy,kmadd(dJ311L,PDstandard3gyy,kmadd(PDstandard11gyy,SQR(J11L),kmadd(PDstandard22gyy,SQR(J21L),kmadd(PDstandard33gyy,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy)),kmul(J21L,kmul(J31L,PDstandard23gyy))),ToReal(2)))))))); + + JacPDstandard11gyz = + kmadd(dJ111L,PDstandard1gyz,kmadd(dJ211L,PDstandard2gyz,kmadd(dJ311L,PDstandard3gyz,kmadd(PDstandard11gyz,SQR(J11L),kmadd(PDstandard22gyz,SQR(J21L),kmadd(PDstandard33gyz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz)),kmul(J21L,kmul(J31L,PDstandard23gyz))),ToReal(2)))))))); + + JacPDstandard11gzz = + kmadd(dJ111L,PDstandard1gzz,kmadd(dJ211L,PDstandard2gzz,kmadd(dJ311L,PDstandard3gzz,kmadd(PDstandard11gzz,SQR(J11L),kmadd(PDstandard22gzz,SQR(J21L),kmadd(PDstandard33gzz,SQR(J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz)),kmul(J21L,kmul(J31L,PDstandard23gzz))),ToReal(2)))))))); + + JacPDstandard22gxx = + kmadd(dJ122L,PDstandard1gxx,kmadd(dJ222L,PDstandard2gxx,kmadd(dJ322L,PDstandard3gxx,kmadd(PDstandard11gxx,SQR(J12L),kmadd(PDstandard22gxx,SQR(J22L),kmadd(PDstandard33gxx,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx)),kmul(J22L,kmul(J32L,PDstandard23gxx))),ToReal(2)))))))); + + JacPDstandard22gxy = + kmadd(dJ122L,PDstandard1gxy,kmadd(dJ222L,PDstandard2gxy,kmadd(dJ322L,PDstandard3gxy,kmadd(PDstandard11gxy,SQR(J12L),kmadd(PDstandard22gxy,SQR(J22L),kmadd(PDstandard33gxy,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxy,kmul(J32L,PDstandard13gxy)),kmul(J22L,kmul(J32L,PDstandard23gxy))),ToReal(2)))))))); + + JacPDstandard22gxz = + kmadd(dJ122L,PDstandard1gxz,kmadd(dJ222L,PDstandard2gxz,kmadd(dJ322L,PDstandard3gxz,kmadd(PDstandard11gxz,SQR(J12L),kmadd(PDstandard22gxz,SQR(J22L),kmadd(PDstandard33gxz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz)),kmul(J22L,kmul(J32L,PDstandard23gxz))),ToReal(2)))))))); + + JacPDstandard22gyz = + kmadd(dJ122L,PDstandard1gyz,kmadd(dJ222L,PDstandard2gyz,kmadd(dJ322L,PDstandard3gyz,kmadd(PDstandard11gyz,SQR(J12L),kmadd(PDstandard22gyz,SQR(J22L),kmadd(PDstandard33gyz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gyz,kmul(J32L,PDstandard13gyz)),kmul(J22L,kmul(J32L,PDstandard23gyz))),ToReal(2)))))))); + + JacPDstandard22gzz = + kmadd(dJ122L,PDstandard1gzz,kmadd(dJ222L,PDstandard2gzz,kmadd(dJ322L,PDstandard3gzz,kmadd(PDstandard11gzz,SQR(J12L),kmadd(PDstandard22gzz,SQR(J22L),kmadd(PDstandard33gzz,SQR(J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz)),kmul(J22L,kmul(J32L,PDstandard23gzz))),ToReal(2)))))))); + + JacPDstandard33gxx = + kmadd(dJ133L,PDstandard1gxx,kmadd(dJ233L,PDstandard2gxx,kmadd(dJ333L,PDstandard3gxx,kmadd(PDstandard11gxx,SQR(J13L),kmadd(PDstandard22gxx,SQR(J23L),kmadd(PDstandard33gxx,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmul(J23L,kmul(J33L,PDstandard23gxx))),ToReal(2)))))))); + + JacPDstandard33gxy = + kmadd(dJ133L,PDstandard1gxy,kmadd(dJ233L,PDstandard2gxy,kmadd(dJ333L,PDstandard3gxy,kmadd(PDstandard11gxy,SQR(J13L),kmadd(PDstandard22gxy,SQR(J23L),kmadd(PDstandard33gxy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmul(J23L,kmul(J33L,PDstandard23gxy))),ToReal(2)))))))); + + JacPDstandard33gxz = + kmadd(dJ133L,PDstandard1gxz,kmadd(dJ233L,PDstandard2gxz,kmadd(dJ333L,PDstandard3gxz,kmadd(PDstandard11gxz,SQR(J13L),kmadd(PDstandard22gxz,SQR(J23L),kmadd(PDstandard33gxz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxz,kmul(J33L,PDstandard13gxz)),kmul(J23L,kmul(J33L,PDstandard23gxz))),ToReal(2)))))))); + + JacPDstandard33gyy = + kmadd(dJ133L,PDstandard1gyy,kmadd(dJ233L,PDstandard2gyy,kmadd(dJ333L,PDstandard3gyy,kmadd(PDstandard11gyy,SQR(J13L),kmadd(PDstandard22gyy,SQR(J23L),kmadd(PDstandard33gyy,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmul(J23L,kmul(J33L,PDstandard23gyy))),ToReal(2)))))))); + + JacPDstandard33gyz = + kmadd(dJ133L,PDstandard1gyz,kmadd(dJ233L,PDstandard2gyz,kmadd(dJ333L,PDstandard3gyz,kmadd(PDstandard11gyz,SQR(J13L),kmadd(PDstandard22gyz,SQR(J23L),kmadd(PDstandard33gyz,SQR(J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gyz,kmul(J33L,PDstandard13gyz)),kmul(J23L,kmul(J33L,PDstandard23gyz))),ToReal(2)))))))); + + JacPDstandard12gxx = + kmadd(J12L,kmadd(J11L,PDstandard11gxx,kmadd(J21L,PDstandard12gxx,kmul(J31L,PDstandard13gxx))),kmadd(J11L,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx)),kmadd(dJ112L,PDstandard1gxx,kmadd(J22L,kmadd(J21L,PDstandard22gxx,kmul(J31L,PDstandard23gxx)),kmadd(dJ212L,PDstandard2gxx,kmadd(J32L,kmadd(J21L,PDstandard23gxx,kmul(J31L,PDstandard33gxx)),kmul(dJ312L,PDstandard3gxx))))))); + + JacPDstandard12gxy = + kmadd(J12L,kmadd(J11L,PDstandard11gxy,kmadd(J21L,PDstandard12gxy,kmul(J31L,PDstandard13gxy))),kmadd(J11L,kmadd(J22L,PDstandard12gxy,kmul(J32L,PDstandard13gxy)),kmadd(dJ112L,PDstandard1gxy,kmadd(J22L,kmadd(J21L,PDstandard22gxy,kmul(J31L,PDstandard23gxy)),kmadd(dJ212L,PDstandard2gxy,kmadd(J32L,kmadd(J21L,PDstandard23gxy,kmul(J31L,PDstandard33gxy)),kmul(dJ312L,PDstandard3gxy))))))); + + JacPDstandard12gxz = + kmadd(J12L,kmadd(J11L,PDstandard11gxz,kmadd(J21L,PDstandard12gxz,kmul(J31L,PDstandard13gxz))),kmadd(J11L,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz)),kmadd(dJ112L,PDstandard1gxz,kmadd(J22L,kmadd(J21L,PDstandard22gxz,kmul(J31L,PDstandard23gxz)),kmadd(dJ212L,PDstandard2gxz,kmadd(J32L,kmadd(J21L,PDstandard23gxz,kmul(J31L,PDstandard33gxz)),kmul(dJ312L,PDstandard3gxz))))))); + + JacPDstandard12gyy = + kmadd(J12L,kmadd(J11L,PDstandard11gyy,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy))),kmadd(J11L,kmadd(J22L,PDstandard12gyy,kmul(J32L,PDstandard13gyy)),kmadd(dJ112L,PDstandard1gyy,kmadd(J22L,kmadd(J21L,PDstandard22gyy,kmul(J31L,PDstandard23gyy)),kmadd(dJ212L,PDstandard2gyy,kmadd(J32L,kmadd(J21L,PDstandard23gyy,kmul(J31L,PDstandard33gyy)),kmul(dJ312L,PDstandard3gyy))))))); + + JacPDstandard12gyz = + kmadd(J12L,kmadd(J11L,PDstandard11gyz,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz))),kmadd(J11L,kmadd(J22L,PDstandard12gyz,kmul(J32L,PDstandard13gyz)),kmadd(dJ112L,PDstandard1gyz,kmadd(J22L,kmadd(J21L,PDstandard22gyz,kmul(J31L,PDstandard23gyz)),kmadd(dJ212L,PDstandard2gyz,kmadd(J32L,kmadd(J21L,PDstandard23gyz,kmul(J31L,PDstandard33gyz)),kmul(dJ312L,PDstandard3gyz))))))); + + JacPDstandard12gzz = + kmadd(J12L,kmadd(J11L,PDstandard11gzz,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz))),kmadd(J11L,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz)),kmadd(dJ112L,PDstandard1gzz,kmadd(J22L,kmadd(J21L,PDstandard22gzz,kmul(J31L,PDstandard23gzz)),kmadd(dJ212L,PDstandard2gzz,kmadd(J32L,kmadd(J21L,PDstandard23gzz,kmul(J31L,PDstandard33gzz)),kmul(dJ312L,PDstandard3gzz))))))); + + JacPDstandard13gxx = + kmadd(J13L,kmadd(J11L,PDstandard11gxx,kmadd(J21L,PDstandard12gxx,kmul(J31L,PDstandard13gxx))),kmadd(J11L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmadd(dJ113L,PDstandard1gxx,kmadd(J23L,kmadd(J21L,PDstandard22gxx,kmul(J31L,PDstandard23gxx)),kmadd(dJ213L,PDstandard2gxx,kmadd(J33L,kmadd(J21L,PDstandard23gxx,kmul(J31L,PDstandard33gxx)),kmul(dJ313L,PDstandard3gxx))))))); + + JacPDstandard13gxy = + kmadd(J13L,kmadd(J11L,PDstandard11gxy,kmadd(J21L,PDstandard12gxy,kmul(J31L,PDstandard13gxy))),kmadd(J11L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmadd(dJ113L,PDstandard1gxy,kmadd(J23L,kmadd(J21L,PDstandard22gxy,kmul(J31L,PDstandard23gxy)),kmadd(dJ213L,PDstandard2gxy,kmadd(J33L,kmadd(J21L,PDstandard23gxy,kmul(J31L,PDstandard33gxy)),kmul(dJ313L,PDstandard3gxy))))))); + + JacPDstandard13gxz = + kmadd(J13L,kmadd(J11L,PDstandard11gxz,kmadd(J21L,PDstandard12gxz,kmul(J31L,PDstandard13gxz))),kmadd(J11L,kmadd(J23L,PDstandard12gxz,kmul(J33L,PDstandard13gxz)),kmadd(dJ113L,PDstandard1gxz,kmadd(J23L,kmadd(J21L,PDstandard22gxz,kmul(J31L,PDstandard23gxz)),kmadd(dJ213L,PDstandard2gxz,kmadd(J33L,kmadd(J21L,PDstandard23gxz,kmul(J31L,PDstandard33gxz)),kmul(dJ313L,PDstandard3gxz))))))); + + JacPDstandard13gyy = + kmadd(J13L,kmadd(J11L,PDstandard11gyy,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy))),kmadd(J11L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmadd(dJ113L,PDstandard1gyy,kmadd(J23L,kmadd(J21L,PDstandard22gyy,kmul(J31L,PDstandard23gyy)),kmadd(dJ213L,PDstandard2gyy,kmadd(J33L,kmadd(J21L,PDstandard23gyy,kmul(J31L,PDstandard33gyy)),kmul(dJ313L,PDstandard3gyy))))))); + + JacPDstandard13gyz = + kmadd(J13L,kmadd(J11L,PDstandard11gyz,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz))),kmadd(J11L,kmadd(J23L,PDstandard12gyz,kmul(J33L,PDstandard13gyz)),kmadd(dJ113L,PDstandard1gyz,kmadd(J23L,kmadd(J21L,PDstandard22gyz,kmul(J31L,PDstandard23gyz)),kmadd(dJ213L,PDstandard2gyz,kmadd(J33L,kmadd(J21L,PDstandard23gyz,kmul(J31L,PDstandard33gyz)),kmul(dJ313L,PDstandard3gyz))))))); + + JacPDstandard13gzz = + kmadd(J13L,kmadd(J11L,PDstandard11gzz,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz))),kmadd(J11L,kmadd(J23L,PDstandard12gzz,kmul(J33L,PDstandard13gzz)),kmadd(dJ113L,PDstandard1gzz,kmadd(J23L,kmadd(J21L,PDstandard22gzz,kmul(J31L,PDstandard23gzz)),kmadd(dJ213L,PDstandard2gzz,kmadd(J33L,kmadd(J21L,PDstandard23gzz,kmul(J31L,PDstandard33gzz)),kmul(dJ313L,PDstandard3gzz))))))); + + JacPDstandard21gxx = + kmadd(J12L,kmadd(J11L,PDstandard11gxx,kmadd(J21L,PDstandard12gxx,kmul(J31L,PDstandard13gxx))),kmadd(J11L,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx)),kmadd(dJ112L,PDstandard1gxx,kmadd(J22L,kmadd(J21L,PDstandard22gxx,kmul(J31L,PDstandard23gxx)),kmadd(dJ212L,PDstandard2gxx,kmadd(J32L,kmadd(J21L,PDstandard23gxx,kmul(J31L,PDstandard33gxx)),kmul(dJ312L,PDstandard3gxx))))))); + + JacPDstandard21gxy = + kmadd(J12L,kmadd(J11L,PDstandard11gxy,kmadd(J21L,PDstandard12gxy,kmul(J31L,PDstandard13gxy))),kmadd(J11L,kmadd(J22L,PDstandard12gxy,kmul(J32L,PDstandard13gxy)),kmadd(dJ112L,PDstandard1gxy,kmadd(J22L,kmadd(J21L,PDstandard22gxy,kmul(J31L,PDstandard23gxy)),kmadd(dJ212L,PDstandard2gxy,kmadd(J32L,kmadd(J21L,PDstandard23gxy,kmul(J31L,PDstandard33gxy)),kmul(dJ312L,PDstandard3gxy))))))); + + JacPDstandard21gxz = + kmadd(J12L,kmadd(J11L,PDstandard11gxz,kmadd(J21L,PDstandard12gxz,kmul(J31L,PDstandard13gxz))),kmadd(J11L,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz)),kmadd(dJ112L,PDstandard1gxz,kmadd(J22L,kmadd(J21L,PDstandard22gxz,kmul(J31L,PDstandard23gxz)),kmadd(dJ212L,PDstandard2gxz,kmadd(J32L,kmadd(J21L,PDstandard23gxz,kmul(J31L,PDstandard33gxz)),kmul(dJ312L,PDstandard3gxz))))))); + + JacPDstandard21gyy = + kmadd(J12L,kmadd(J11L,PDstandard11gyy,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy))),kmadd(J11L,kmadd(J22L,PDstandard12gyy,kmul(J32L,PDstandard13gyy)),kmadd(dJ112L,PDstandard1gyy,kmadd(J22L,kmadd(J21L,PDstandard22gyy,kmul(J31L,PDstandard23gyy)),kmadd(dJ212L,PDstandard2gyy,kmadd(J32L,kmadd(J21L,PDstandard23gyy,kmul(J31L,PDstandard33gyy)),kmul(dJ312L,PDstandard3gyy))))))); + + JacPDstandard21gyz = + kmadd(J12L,kmadd(J11L,PDstandard11gyz,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz))),kmadd(J11L,kmadd(J22L,PDstandard12gyz,kmul(J32L,PDstandard13gyz)),kmadd(dJ112L,PDstandard1gyz,kmadd(J22L,kmadd(J21L,PDstandard22gyz,kmul(J31L,PDstandard23gyz)),kmadd(dJ212L,PDstandard2gyz,kmadd(J32L,kmadd(J21L,PDstandard23gyz,kmul(J31L,PDstandard33gyz)),kmul(dJ312L,PDstandard3gyz))))))); + + JacPDstandard21gzz = + kmadd(J12L,kmadd(J11L,PDstandard11gzz,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz))),kmadd(J11L,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz)),kmadd(dJ112L,PDstandard1gzz,kmadd(J22L,kmadd(J21L,PDstandard22gzz,kmul(J31L,PDstandard23gzz)),kmadd(dJ212L,PDstandard2gzz,kmadd(J32L,kmadd(J21L,PDstandard23gzz,kmul(J31L,PDstandard33gzz)),kmul(dJ312L,PDstandard3gzz))))))); + + JacPDstandard23gxx = + kmadd(J13L,kmadd(J12L,PDstandard11gxx,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx))),kmadd(J12L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmadd(dJ123L,PDstandard1gxx,kmadd(J23L,kmadd(J22L,PDstandard22gxx,kmul(J32L,PDstandard23gxx)),kmadd(dJ223L,PDstandard2gxx,kmadd(J33L,kmadd(J22L,PDstandard23gxx,kmul(J32L,PDstandard33gxx)),kmul(dJ323L,PDstandard3gxx))))))); + + JacPDstandard23gxy = + kmadd(J13L,kmadd(J12L,PDstandard11gxy,kmadd(J22L,PDstandard12gxy,kmul(J32L,PDstandard13gxy))),kmadd(J12L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmadd(dJ123L,PDstandard1gxy,kmadd(J23L,kmadd(J22L,PDstandard22gxy,kmul(J32L,PDstandard23gxy)),kmadd(dJ223L,PDstandard2gxy,kmadd(J33L,kmadd(J22L,PDstandard23gxy,kmul(J32L,PDstandard33gxy)),kmul(dJ323L,PDstandard3gxy))))))); + + JacPDstandard23gxz = + kmadd(J13L,kmadd(J12L,PDstandard11gxz,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz))),kmadd(J12L,kmadd(J23L,PDstandard12gxz,kmul(J33L,PDstandard13gxz)),kmadd(dJ123L,PDstandard1gxz,kmadd(J23L,kmadd(J22L,PDstandard22gxz,kmul(J32L,PDstandard23gxz)),kmadd(dJ223L,PDstandard2gxz,kmadd(J33L,kmadd(J22L,PDstandard23gxz,kmul(J32L,PDstandard33gxz)),kmul(dJ323L,PDstandard3gxz))))))); + + JacPDstandard23gyy = + kmadd(J13L,kmadd(J12L,PDstandard11gyy,kmadd(J22L,PDstandard12gyy,kmul(J32L,PDstandard13gyy))),kmadd(J12L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmadd(dJ123L,PDstandard1gyy,kmadd(J23L,kmadd(J22L,PDstandard22gyy,kmul(J32L,PDstandard23gyy)),kmadd(dJ223L,PDstandard2gyy,kmadd(J33L,kmadd(J22L,PDstandard23gyy,kmul(J32L,PDstandard33gyy)),kmul(dJ323L,PDstandard3gyy))))))); + + JacPDstandard23gyz = + kmadd(J13L,kmadd(J12L,PDstandard11gyz,kmadd(J22L,PDstandard12gyz,kmul(J32L,PDstandard13gyz))),kmadd(J12L,kmadd(J23L,PDstandard12gyz,kmul(J33L,PDstandard13gyz)),kmadd(dJ123L,PDstandard1gyz,kmadd(J23L,kmadd(J22L,PDstandard22gyz,kmul(J32L,PDstandard23gyz)),kmadd(dJ223L,PDstandard2gyz,kmadd(J33L,kmadd(J22L,PDstandard23gyz,kmul(J32L,PDstandard33gyz)),kmul(dJ323L,PDstandard3gyz))))))); + + JacPDstandard23gzz = + kmadd(J13L,kmadd(J12L,PDstandard11gzz,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz))),kmadd(J12L,kmadd(J23L,PDstandard12gzz,kmul(J33L,PDstandard13gzz)),kmadd(dJ123L,PDstandard1gzz,kmadd(J23L,kmadd(J22L,PDstandard22gzz,kmul(J32L,PDstandard23gzz)),kmadd(dJ223L,PDstandard2gzz,kmadd(J33L,kmadd(J22L,PDstandard23gzz,kmul(J32L,PDstandard33gzz)),kmul(dJ323L,PDstandard3gzz))))))); + + JacPDstandard31gxx = + kmadd(J13L,kmadd(J11L,PDstandard11gxx,kmadd(J21L,PDstandard12gxx,kmul(J31L,PDstandard13gxx))),kmadd(J11L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmadd(dJ113L,PDstandard1gxx,kmadd(J23L,kmadd(J21L,PDstandard22gxx,kmul(J31L,PDstandard23gxx)),kmadd(dJ213L,PDstandard2gxx,kmadd(J33L,kmadd(J21L,PDstandard23gxx,kmul(J31L,PDstandard33gxx)),kmul(dJ313L,PDstandard3gxx))))))); + + JacPDstandard31gxy = + kmadd(J13L,kmadd(J11L,PDstandard11gxy,kmadd(J21L,PDstandard12gxy,kmul(J31L,PDstandard13gxy))),kmadd(J11L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmadd(dJ113L,PDstandard1gxy,kmadd(J23L,kmadd(J21L,PDstandard22gxy,kmul(J31L,PDstandard23gxy)),kmadd(dJ213L,PDstandard2gxy,kmadd(J33L,kmadd(J21L,PDstandard23gxy,kmul(J31L,PDstandard33gxy)),kmul(dJ313L,PDstandard3gxy))))))); + + JacPDstandard31gxz = + kmadd(J13L,kmadd(J11L,PDstandard11gxz,kmadd(J21L,PDstandard12gxz,kmul(J31L,PDstandard13gxz))),kmadd(J11L,kmadd(J23L,PDstandard12gxz,kmul(J33L,PDstandard13gxz)),kmadd(dJ113L,PDstandard1gxz,kmadd(J23L,kmadd(J21L,PDstandard22gxz,kmul(J31L,PDstandard23gxz)),kmadd(dJ213L,PDstandard2gxz,kmadd(J33L,kmadd(J21L,PDstandard23gxz,kmul(J31L,PDstandard33gxz)),kmul(dJ313L,PDstandard3gxz))))))); + + JacPDstandard31gyy = + kmadd(J13L,kmadd(J11L,PDstandard11gyy,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy))),kmadd(J11L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmadd(dJ113L,PDstandard1gyy,kmadd(J23L,kmadd(J21L,PDstandard22gyy,kmul(J31L,PDstandard23gyy)),kmadd(dJ213L,PDstandard2gyy,kmadd(J33L,kmadd(J21L,PDstandard23gyy,kmul(J31L,PDstandard33gyy)),kmul(dJ313L,PDstandard3gyy))))))); + + JacPDstandard31gyz = + kmadd(J13L,kmadd(J11L,PDstandard11gyz,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz))),kmadd(J11L,kmadd(J23L,PDstandard12gyz,kmul(J33L,PDstandard13gyz)),kmadd(dJ113L,PDstandard1gyz,kmadd(J23L,kmadd(J21L,PDstandard22gyz,kmul(J31L,PDstandard23gyz)),kmadd(dJ213L,PDstandard2gyz,kmadd(J33L,kmadd(J21L,PDstandard23gyz,kmul(J31L,PDstandard33gyz)),kmul(dJ313L,PDstandard3gyz))))))); + + JacPDstandard31gzz = + kmadd(J13L,kmadd(J11L,PDstandard11gzz,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz))),kmadd(J11L,kmadd(J23L,PDstandard12gzz,kmul(J33L,PDstandard13gzz)),kmadd(dJ113L,PDstandard1gzz,kmadd(J23L,kmadd(J21L,PDstandard22gzz,kmul(J31L,PDstandard23gzz)),kmadd(dJ213L,PDstandard2gzz,kmadd(J33L,kmadd(J21L,PDstandard23gzz,kmul(J31L,PDstandard33gzz)),kmul(dJ313L,PDstandard3gzz))))))); + + JacPDstandard32gxx = + kmadd(J13L,kmadd(J12L,PDstandard11gxx,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx))),kmadd(J12L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmadd(dJ123L,PDstandard1gxx,kmadd(J23L,kmadd(J22L,PDstandard22gxx,kmul(J32L,PDstandard23gxx)),kmadd(dJ223L,PDstandard2gxx,kmadd(J33L,kmadd(J22L,PDstandard23gxx,kmul(J32L,PDstandard33gxx)),kmul(dJ323L,PDstandard3gxx))))))); + + JacPDstandard32gxy = + kmadd(J13L,kmadd(J12L,PDstandard11gxy,kmadd(J22L,PDstandard12gxy,kmul(J32L,PDstandard13gxy))),kmadd(J12L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmadd(dJ123L,PDstandard1gxy,kmadd(J23L,kmadd(J22L,PDstandard22gxy,kmul(J32L,PDstandard23gxy)),kmadd(dJ223L,PDstandard2gxy,kmadd(J33L,kmadd(J22L,PDstandard23gxy,kmul(J32L,PDstandard33gxy)),kmul(dJ323L,PDstandard3gxy))))))); + + JacPDstandard32gxz = + kmadd(J13L,kmadd(J12L,PDstandard11gxz,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz))),kmadd(J12L,kmadd(J23L,PDstandard12gxz,kmul(J33L,PDstandard13gxz)),kmadd(dJ123L,PDstandard1gxz,kmadd(J23L,kmadd(J22L,PDstandard22gxz,kmul(J32L,PDstandard23gxz)),kmadd(dJ223L,PDstandard2gxz,kmadd(J33L,kmadd(J22L,PDstandard23gxz,kmul(J32L,PDstandard33gxz)),kmul(dJ323L,PDstandard3gxz))))))); + + JacPDstandard32gyy = + kmadd(J13L,kmadd(J12L,PDstandard11gyy,kmadd(J22L,PDstandard12gyy,kmul(J32L,PDstandard13gyy))),kmadd(J12L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmadd(dJ123L,PDstandard1gyy,kmadd(J23L,kmadd(J22L,PDstandard22gyy,kmul(J32L,PDstandard23gyy)),kmadd(dJ223L,PDstandard2gyy,kmadd(J33L,kmadd(J22L,PDstandard23gyy,kmul(J32L,PDstandard33gyy)),kmul(dJ323L,PDstandard3gyy))))))); + + JacPDstandard32gyz = + kmadd(J13L,kmadd(J12L,PDstandard11gyz,kmadd(J22L,PDstandard12gyz,kmul(J32L,PDstandard13gyz))),kmadd(J12L,kmadd(J23L,PDstandard12gyz,kmul(J33L,PDstandard13gyz)),kmadd(dJ123L,PDstandard1gyz,kmadd(J23L,kmadd(J22L,PDstandard22gyz,kmul(J32L,PDstandard23gyz)),kmadd(dJ223L,PDstandard2gyz,kmadd(J33L,kmadd(J22L,PDstandard23gyz,kmul(J32L,PDstandard33gyz)),kmul(dJ323L,PDstandard3gyz))))))); + + JacPDstandard32gzz = + kmadd(J13L,kmadd(J12L,PDstandard11gzz,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz))),kmadd(J12L,kmadd(J23L,PDstandard12gzz,kmul(J33L,PDstandard13gzz)),kmadd(dJ123L,PDstandard1gzz,kmadd(J23L,kmadd(J22L,PDstandard22gzz,kmul(J32L,PDstandard23gzz)),kmadd(dJ223L,PDstandard2gzz,kmadd(J33L,kmadd(J22L,PDstandard23gzz,kmul(J32L,PDstandard33gzz)),kmul(dJ323L,PDstandard3gzz))))))); + } + else + { + JacPDstandard1gxx = PDstandard1gxx; + + JacPDstandard1gxy = PDstandard1gxy; + + JacPDstandard1gxz = PDstandard1gxz; + + JacPDstandard1gyy = PDstandard1gyy; + + JacPDstandard1gyz = PDstandard1gyz; + + JacPDstandard1gzz = PDstandard1gzz; + + JacPDstandard1kxy = PDstandard1kxy; + + JacPDstandard1kxz = PDstandard1kxz; + + JacPDstandard1kyy = PDstandard1kyy; + + JacPDstandard1kyz = PDstandard1kyz; + + JacPDstandard1kzz = PDstandard1kzz; + + JacPDstandard2gxx = PDstandard2gxx; + + JacPDstandard2gxy = PDstandard2gxy; + + JacPDstandard2gxz = PDstandard2gxz; + + JacPDstandard2gyy = PDstandard2gyy; + + JacPDstandard2gyz = PDstandard2gyz; + + JacPDstandard2gzz = PDstandard2gzz; + + JacPDstandard2kxx = PDstandard2kxx; + + JacPDstandard2kxy = PDstandard2kxy; + + JacPDstandard2kxz = PDstandard2kxz; + + JacPDstandard2kyz = PDstandard2kyz; + + JacPDstandard2kzz = PDstandard2kzz; + + JacPDstandard3gxx = PDstandard3gxx; + + JacPDstandard3gxy = PDstandard3gxy; + + JacPDstandard3gxz = PDstandard3gxz; + + JacPDstandard3gyy = PDstandard3gyy; + + JacPDstandard3gyz = PDstandard3gyz; + + JacPDstandard3gzz = PDstandard3gzz; + + JacPDstandard3kxx = PDstandard3kxx; + + JacPDstandard3kxy = PDstandard3kxy; + + JacPDstandard3kxz = PDstandard3kxz; + + JacPDstandard3kyy = PDstandard3kyy; + + JacPDstandard3kyz = PDstandard3kyz; + + JacPDstandard11gxy = PDstandard11gxy; + + JacPDstandard11gxz = PDstandard11gxz; + + JacPDstandard11gyy = PDstandard11gyy; + + JacPDstandard11gyz = PDstandard11gyz; + + JacPDstandard11gzz = PDstandard11gzz; + + JacPDstandard22gxx = PDstandard22gxx; + + JacPDstandard22gxy = PDstandard22gxy; + + JacPDstandard22gxz = PDstandard22gxz; + + JacPDstandard22gyz = PDstandard22gyz; + + JacPDstandard22gzz = PDstandard22gzz; + + JacPDstandard33gxx = PDstandard33gxx; + + JacPDstandard33gxy = PDstandard33gxy; + + JacPDstandard33gxz = PDstandard33gxz; + + JacPDstandard33gyy = PDstandard33gyy; + + JacPDstandard33gyz = PDstandard33gyz; + + JacPDstandard12gxx = PDstandard12gxx; + + JacPDstandard12gxy = PDstandard12gxy; + + JacPDstandard12gxz = PDstandard12gxz; + + JacPDstandard12gyy = PDstandard12gyy; + + JacPDstandard12gyz = PDstandard12gyz; + + JacPDstandard12gzz = PDstandard12gzz; + + JacPDstandard13gxx = PDstandard13gxx; + + JacPDstandard13gxy = PDstandard13gxy; + + JacPDstandard13gxz = PDstandard13gxz; + + JacPDstandard13gyy = PDstandard13gyy; + + JacPDstandard13gyz = PDstandard13gyz; + + JacPDstandard13gzz = PDstandard13gzz; + + JacPDstandard21gxx = PDstandard12gxx; + + JacPDstandard21gxy = PDstandard12gxy; + + JacPDstandard21gxz = PDstandard12gxz; + + JacPDstandard21gyy = PDstandard12gyy; + + JacPDstandard21gyz = PDstandard12gyz; + + JacPDstandard21gzz = PDstandard12gzz; + + JacPDstandard23gxx = PDstandard23gxx; + + JacPDstandard23gxy = PDstandard23gxy; + + JacPDstandard23gxz = PDstandard23gxz; + + JacPDstandard23gyy = PDstandard23gyy; + + JacPDstandard23gyz = PDstandard23gyz; + + JacPDstandard23gzz = PDstandard23gzz; + + JacPDstandard31gxx = PDstandard13gxx; + + JacPDstandard31gxy = PDstandard13gxy; + + JacPDstandard31gxz = PDstandard13gxz; + + JacPDstandard31gyy = PDstandard13gyy; + + JacPDstandard31gyz = PDstandard13gyz; + + JacPDstandard31gzz = PDstandard13gzz; + + JacPDstandard32gxx = PDstandard23gxx; + + JacPDstandard32gxy = PDstandard23gxy; + + JacPDstandard32gxz = PDstandard23gxz; + + JacPDstandard32gyy = PDstandard23gyy; + + JacPDstandard32gyz = PDstandard23gyz; + + JacPDstandard32gzz = PDstandard23gzz; + } + + 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,JacPDstandard1gxx,knmsub(gInv12,JacPDstandard2gxx,kmsub(kmadd(gInv12,JacPDstandard1gxy,kmul(gInv13,JacPDstandard1gxz)),ToReal(2),kmul(gInv13,JacPDstandard3gxx))))); + + CCTK_REAL_VEC gamma211 = + kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard1gxx,knmsub(gInv22,JacPDstandard2gxx,kmsub(kmadd(gInv22,JacPDstandard1gxy,kmul(gInv23,JacPDstandard1gxz)),ToReal(2),kmul(gInv23,JacPDstandard3gxx))))); + + CCTK_REAL_VEC gamma311 = + kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard1gxx,knmsub(gInv32,JacPDstandard2gxx,kmsub(kmadd(gInv32,JacPDstandard1gxy,kmul(gInv33,JacPDstandard1gxz)),ToReal(2),kmul(gInv33,JacPDstandard3gxx))))); + + CCTK_REAL_VEC gamma121 = + kmul(kmadd(gInv12,JacPDstandard1gyy,kmadd(gInv11,JacPDstandard2gxx,kmul(gInv13,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma221 = + kmul(kmadd(gInv22,JacPDstandard1gyy,kmadd(gInv21,JacPDstandard2gxx,kmul(gInv23,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma321 = + kmul(kmadd(gInv32,JacPDstandard1gyy,kmadd(gInv31,JacPDstandard2gxx,kmul(gInv33,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma131 = + kmul(kmadd(gInv13,JacPDstandard1gzz,kmadd(gInv11,JacPDstandard3gxx,kmul(gInv12,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma231 = + kmul(kmadd(gInv23,JacPDstandard1gzz,kmadd(gInv21,JacPDstandard3gxx,kmul(gInv22,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma331 = + kmul(kmadd(gInv33,JacPDstandard1gzz,kmadd(gInv31,JacPDstandard3gxx,kmul(gInv32,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma122 = + kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2gyy,kmadd(gInv11,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv13,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + + CCTK_REAL_VEC gamma222 = + kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2gyy,kmadd(gInv21,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv23,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + + CCTK_REAL_VEC gamma322 = + kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2gyy,kmadd(gInv31,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv33,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + + CCTK_REAL_VEC gamma132 = + kmul(kmadd(gInv13,JacPDstandard2gzz,kmadd(gInv12,JacPDstandard3gyy,kmul(gInv11,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma232 = + kmul(kmadd(gInv23,JacPDstandard2gzz,kmadd(gInv22,JacPDstandard3gyy,kmul(gInv21,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma332 = + kmul(kmadd(gInv33,JacPDstandard2gzz,kmadd(gInv32,JacPDstandard3gyy,kmul(gInv31,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + + CCTK_REAL_VEC gamma133 = + kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard3gzz,kmadd(gInv11,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv12,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + + CCTK_REAL_VEC gamma233 = + kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard3gzz,kmadd(gInv21,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv22,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + + CCTK_REAL_VEC gamma333 = + kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard3gzz,kmadd(gInv31,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv32,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + + 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(JacPDstandard21gxx,JacPDstandard12gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1113 = + kmul(ksub(JacPDstandard31gxx,JacPDstandard13gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1121 = + kmul(ksub(JacPDstandard21gxx,JacPDstandard12gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1122 = ksub(JacPDstandard21gxy,JacPDstandard12gxy); + + CCTK_REAL_VEC R1123 = + kmul(kadd(JacPDstandard21gxz,ksub(JacPDstandard31gxy,kadd(JacPDstandard13gxy,JacPDstandard12gxz))),ToReal(0.5)); + + CCTK_REAL_VEC R1131 = + kmul(ksub(JacPDstandard31gxx,JacPDstandard13gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1132 = + kmul(kadd(JacPDstandard21gxz,ksub(JacPDstandard31gxy,kadd(JacPDstandard13gxy,JacPDstandard12gxz))),ToReal(0.5)); + + CCTK_REAL_VEC R1133 = ksub(JacPDstandard31gxz,JacPDstandard13gxz); + + CCTK_REAL_VEC R1211 = + kmul(ksub(JacPDstandard12gxx,JacPDstandard21gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1212 = + kmul(ToReal(0.5),kadd(JacPDstandard12gxy,kadd(JacPDstandard21gxy,kmadd(kmadd(gamma122,kmul(gamma211,gxyL),kmadd(gamma122,kmul(gamma311,gxzL),kmadd(gamma111,kmadd(gamma122,gxxL,kmadd(gamma222,gxyL,kmul(gamma322,gxzL))),kmadd(gamma211,kmul(gamma222,gyyL),kmadd(gamma222,kmul(gamma311,gyzL),kmadd(gamma211,kmul(gamma322,gyzL),kmul(gamma311,kmul(gamma322,gzzL)))))))),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(gamma121,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmul(gamma221,kmul(gamma321,gyzL))),ToReal(4))))),JacPDstandard22gxx),JacPDstandard11gyy))))); + + CCTK_REAL_VEC R1213 = + kmadd(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard12gxz,kadd(JacPDstandard31gxy,kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma221,gxyL),kmadd(gamma131,kmul(gamma321,gxzL),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma211,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma331,ToReal(2)))))),JacPDstandard23gxx),JacPDstandard11gyz)))))); + + CCTK_REAL_VEC R1221 = + kmadd(gamma111,kmadd(gamma122,gxxL,kmadd(gamma222,gxyL,kmul(gamma322,gxzL))),kmadd(gamma121,kmul(kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard11gyy,kadd(JacPDstandard22gxx,kmadd(gyyL,kmul(SQR(gamma221),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma122,kmul(gamma211,gxyL),kmadd(gamma122,kmul(gamma311,gxzL),kmadd(gamma211,kmul(gamma222,gyyL),kmadd(gamma222,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma322,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma321,ToReal(-4)),kmul(gamma211,kmul(gamma322,ToReal(2)))))),JacPDstandard21gxy),JacPDstandard12gxy))))),kmul(gxxL,SQR(gamma121))))); + + CCTK_REAL_VEC R1222 = + kmul(ksub(JacPDstandard21gyy,JacPDstandard12gyy),ToReal(0.5)); + + CCTK_REAL_VEC R1223 = + kmadd(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard22gxz,kadd(JacPDstandard31gyy,kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma222,gxyL),kmadd(gamma131,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma332,ToReal(-2)),kmul(gamma222,kmul(gamma331,ToReal(2)))))),JacPDstandard23gxy),JacPDstandard12gyz)))))); + + CCTK_REAL_VEC R1231 = + kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard11gyz,kadd(JacPDstandard32gxx,kmadd(kmadd(gamma131,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma331,ToReal(-2)),kmul(gamma211,kmul(gamma332,ToReal(2)))))),JacPDstandard21gxz),JacPDstandard13gxy)))),kmul(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R1232 = + kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard21gyz,kadd(JacPDstandard32gxy,kmadd(kmadd(gamma131,kmadd(gamma222,gxyL,kmul(gamma322,gxzL)),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma222,kmul(gamma331,ToReal(-2)),kmul(gamma221,kmul(gamma332,ToReal(2)))))),JacPDstandard22gxz),JacPDstandard13gyy)))),kmul(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R1233 = + kmul(kadd(JacPDstandard31gyz,ksub(JacPDstandard32gxz,kadd(JacPDstandard23gxz,JacPDstandard13gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R1311 = + kmul(ksub(JacPDstandard13gxx,JacPDstandard31gxx),ToReal(0.5)); + + CCTK_REAL_VEC R1312 = + kmadd(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gxy,kadd(JacPDstandard21gxz,kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma221,gxyL),kmadd(gamma131,kmul(gamma321,gxzL),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma211,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma331,ToReal(2)))))),JacPDstandard32gxx),JacPDstandard11gyz)))))); + + CCTK_REAL_VEC R1313 = + kmul(ToReal(0.5),kadd(JacPDstandard13gxz,kadd(JacPDstandard31gxz,kmadd(kmadd(gamma133,kmul(gamma211,gxyL),kmadd(gamma133,kmul(gamma311,gxzL),kmadd(gamma111,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma211,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma311,gyzL),kmadd(gamma211,kmul(gamma333,gyzL),kmul(gamma311,kmul(gamma333,gzzL)))))))),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(gamma131,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmul(gamma231,kmul(gamma331,gyzL))),ToReal(4))))),JacPDstandard33gxx),JacPDstandard11gzz))))); + + CCTK_REAL_VEC R1321 = + kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard11gyz,kadd(JacPDstandard23gxx,kmadd(kmadd(gamma131,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma331,ToReal(-2)),kmul(gamma211,kmul(gamma332,ToReal(2)))))),JacPDstandard31gxy),JacPDstandard12gxz)))),kmul(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R1322 = + kmul(kadd(JacPDstandard21gyz,ksub(JacPDstandard23gxy,kadd(JacPDstandard32gxy,JacPDstandard12gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R1323 = + kmadd(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard23gxz,kadd(JacPDstandard31gyz,kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmul(gamma231,gxyL),kmadd(gamma132,kmul(gamma331,gxzL),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma333,ToReal(-2)),kmul(gamma231,kmul(gamma332,ToReal(2)))))),JacPDstandard33gxy),JacPDstandard12gzz)))))); + + CCTK_REAL_VEC R1331 = + kmadd(gamma111,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma131,kmul(kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard11gzz,kadd(JacPDstandard33gxx,kmadd(gyyL,kmul(SQR(gamma231),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma133,kmul(gamma211,gxyL),kmadd(gamma133,kmul(gamma311,gxzL),kmadd(gamma211,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma333,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma331,ToReal(-4)),kmul(gamma211,kmul(gamma333,ToReal(2)))))),JacPDstandard31gxz),JacPDstandard13gxz))))),kmul(gxxL,SQR(gamma131))))); + + CCTK_REAL_VEC R1332 = + kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard21gzz,kadd(JacPDstandard33gxy,kmadd(kmadd(gamma132,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma333,ToReal(2)))))),JacPDstandard32gxz),JacPDstandard13gyz)))),kmul(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)))))); + + CCTK_REAL_VEC R1333 = + kmul(ksub(JacPDstandard31gzz,JacPDstandard13gzz),ToReal(0.5)); + + CCTK_REAL_VEC R2111 = + kmul(ksub(JacPDstandard12gxx,JacPDstandard21gxx),ToReal(0.5)); + + CCTK_REAL_VEC R2112 = + kmadd(gamma111,kmadd(gamma122,gxxL,kmadd(gamma222,gxyL,kmul(gamma322,gxzL))),kmadd(gamma121,kmul(kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard11gyy,kadd(JacPDstandard22gxx,kmadd(gyyL,kmul(SQR(gamma221),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma321),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma122,kmul(gamma211,gxyL),kmadd(gamma122,kmul(gamma311,gxzL),kmadd(gamma211,kmul(gamma222,gyyL),kmadd(gamma222,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma322,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma321,ToReal(-4)),kmul(gamma211,kmul(gamma322,ToReal(2)))))),JacPDstandard21gxy),JacPDstandard12gxy))))),kmul(gxxL,SQR(gamma121))))); + + CCTK_REAL_VEC R2113 = + kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard11gyz,kadd(JacPDstandard32gxx,kmadd(kmadd(gamma131,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma331,ToReal(-2)),kmul(gamma211,kmul(gamma332,ToReal(2)))))),JacPDstandard21gxz),JacPDstandard13gxy)))),kmul(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R2121 = + kmul(ToReal(0.5),kadd(JacPDstandard12gxy,kadd(JacPDstandard21gxy,kmadd(kmadd(gamma122,kmul(gamma211,gxyL),kmadd(gamma122,kmul(gamma311,gxzL),kmadd(gamma111,kmadd(gamma122,gxxL,kmadd(gamma222,gxyL,kmul(gamma322,gxzL))),kmadd(gamma211,kmul(gamma222,gyyL),kmadd(gamma222,kmul(gamma311,gyzL),kmadd(gamma211,kmul(gamma322,gyzL),kmul(gamma311,kmul(gamma322,gzzL)))))))),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(gamma121,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmul(gamma221,kmul(gamma321,gyzL))),ToReal(4))))),JacPDstandard22gxx),JacPDstandard11gyy))))); + + CCTK_REAL_VEC R2122 = + kmul(ksub(JacPDstandard21gyy,JacPDstandard12gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2123 = + kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard21gyz,kadd(JacPDstandard32gxy,kmadd(kmadd(gamma131,kmadd(gamma222,gxyL,kmul(gamma322,gxzL)),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma222,kmul(gamma331,ToReal(-2)),kmul(gamma221,kmul(gamma332,ToReal(2)))))),JacPDstandard22gxz),JacPDstandard13gyy)))),kmul(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R2131 = + kmadd(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard12gxz,kadd(JacPDstandard31gxy,kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma221,gxyL),kmadd(gamma131,kmul(gamma321,gxzL),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma211,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma331,ToReal(2)))))),JacPDstandard23gxx),JacPDstandard11gyz)))))); + + CCTK_REAL_VEC R2132 = + kmadd(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard22gxz,kadd(JacPDstandard31gyy,kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma222,gxyL),kmadd(gamma131,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma332,ToReal(-2)),kmul(gamma222,kmul(gamma331,ToReal(2)))))),JacPDstandard23gxy),JacPDstandard12gyz)))))); + + CCTK_REAL_VEC R2133 = + kmul(kadd(JacPDstandard31gyz,ksub(JacPDstandard32gxz,kadd(JacPDstandard23gxz,JacPDstandard13gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R2211 = ksub(JacPDstandard12gxy,JacPDstandard21gxy); + + CCTK_REAL_VEC R2212 = + kmul(ksub(JacPDstandard12gyy,JacPDstandard21gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2213 = + kmul(kadd(JacPDstandard12gyz,ksub(JacPDstandard32gxy,kadd(JacPDstandard23gxy,JacPDstandard21gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R2221 = + kmul(ksub(JacPDstandard12gyy,JacPDstandard21gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2222 = ToReal(0); + + CCTK_REAL_VEC R2223 = + kmul(ksub(JacPDstandard32gyy,JacPDstandard23gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2231 = + kmul(kadd(JacPDstandard12gyz,ksub(JacPDstandard32gxy,kadd(JacPDstandard23gxy,JacPDstandard21gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R2232 = + kmul(ksub(JacPDstandard32gyy,JacPDstandard23gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2233 = ksub(JacPDstandard32gyz,JacPDstandard23gyz); + + CCTK_REAL_VEC R2311 = + kmul(kadd(JacPDstandard12gxz,ksub(JacPDstandard13gxy,kadd(JacPDstandard31gxy,JacPDstandard21gxz))),ToReal(0.5)); + + CCTK_REAL_VEC R2312 = + kmadd(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gyy,kadd(JacPDstandard22gxz,kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma222,gxyL),kmadd(gamma131,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma332,ToReal(-2)),kmul(gamma222,kmul(gamma331,ToReal(2)))))),JacPDstandard32gxy),JacPDstandard21gyz)))))); + + CCTK_REAL_VEC R2313 = + kmadd(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gyz,kadd(JacPDstandard32gxz,kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmul(gamma231,gxyL),kmadd(gamma132,kmul(gamma331,gxzL),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma333,ToReal(-2)),kmul(gamma231,kmul(gamma332,ToReal(2)))))),JacPDstandard33gxy),JacPDstandard21gzz)))))); + + CCTK_REAL_VEC R2321 = + kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard12gyz,kadd(JacPDstandard23gxy,kmadd(kmadd(gamma131,kmadd(gamma222,gxyL,kmul(gamma322,gxzL)),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma222,kmul(gamma331,ToReal(-2)),kmul(gamma221,kmul(gamma332,ToReal(2)))))),JacPDstandard31gyy),JacPDstandard22gxz)))),kmul(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R2322 = + kmul(ksub(JacPDstandard23gyy,JacPDstandard32gyy),ToReal(0.5)); + + CCTK_REAL_VEC R2323 = + kmul(ToReal(0.5),kadd(JacPDstandard23gyz,kadd(JacPDstandard32gyz,kmadd(kmadd(gamma133,kmul(gamma222,gxyL),kmadd(gamma133,kmul(gamma322,gxzL),kmadd(gamma122,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma222,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma322,gyzL),kmadd(gamma222,kmul(gamma333,gyzL),kmul(gamma322,kmul(gamma333,gzzL)))))))),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(gamma132,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)),kmul(gamma232,kmul(gamma332,gyzL))),ToReal(4))))),JacPDstandard33gyy),JacPDstandard22gzz))))); + + CCTK_REAL_VEC R2331 = + kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard12gzz,kadd(JacPDstandard33gxy,kmadd(kmadd(gamma132,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma333,ToReal(2)))))),JacPDstandard31gyz),JacPDstandard23gxz)))),kmul(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)))))); + + CCTK_REAL_VEC R2332 = + kmadd(gamma122,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma132,kmul(kmadd(gamma232,gxyL,kmul(gamma332,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard22gzz,kadd(JacPDstandard33gyy,kmadd(gyyL,kmul(SQR(gamma232),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma133,kmul(gamma222,gxyL),kmadd(gamma133,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma333,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma232,kmul(gamma332,ToReal(-4)),kmul(gamma222,kmul(gamma333,ToReal(2)))))),JacPDstandard32gyz),JacPDstandard23gyz))))),kmul(gxxL,SQR(gamma132))))); + + CCTK_REAL_VEC R2333 = + kmul(ksub(JacPDstandard32gzz,JacPDstandard23gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3111 = + kmul(ksub(JacPDstandard13gxx,JacPDstandard31gxx),ToReal(0.5)); + + CCTK_REAL_VEC R3112 = + kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard11gyz,kadd(JacPDstandard23gxx,kmadd(kmadd(gamma131,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma331,ToReal(-2)),kmul(gamma211,kmul(gamma332,ToReal(2)))))),JacPDstandard31gxy),JacPDstandard12gxz)))),kmul(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R3113 = + kmadd(gamma111,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma131,kmul(kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard11gzz,kadd(JacPDstandard33gxx,kmadd(gyyL,kmul(SQR(gamma231),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma331),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma133,kmul(gamma211,gxyL),kmadd(gamma133,kmul(gamma311,gxzL),kmadd(gamma211,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma333,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma331,ToReal(-4)),kmul(gamma211,kmul(gamma333,ToReal(2)))))),JacPDstandard31gxz),JacPDstandard13gxz))))),kmul(gxxL,SQR(gamma131))))); + + CCTK_REAL_VEC R3121 = + kmadd(gamma121,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gxy,kadd(JacPDstandard21gxz,kmadd(kmadd(gamma132,kmadd(gamma211,gxyL,kmul(gamma311,gxzL)),kmadd(gamma111,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma211,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma311,gyzL),kmul(gamma311,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma221,gxyL),kmadd(gamma131,kmul(gamma321,gxzL),kmadd(gamma221,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma211,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma331,ToReal(2)))))),JacPDstandard32gxx),JacPDstandard11gyz)))))); + + CCTK_REAL_VEC R3122 = + kmul(kadd(JacPDstandard21gyz,ksub(JacPDstandard23gxy,kadd(JacPDstandard32gxy,JacPDstandard12gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R3123 = + kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard21gzz,kadd(JacPDstandard33gxy,kmadd(kmadd(gamma132,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma333,ToReal(2)))))),JacPDstandard32gxz),JacPDstandard13gyz)))),kmul(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)))))); + + CCTK_REAL_VEC R3131 = + kmul(ToReal(0.5),kadd(JacPDstandard13gxz,kadd(JacPDstandard31gxz,kmadd(kmadd(gamma133,kmul(gamma211,gxyL),kmadd(gamma133,kmul(gamma311,gxzL),kmadd(gamma111,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma211,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma311,gyzL),kmadd(gamma211,kmul(gamma333,gyzL),kmul(gamma311,kmul(gamma333,gzzL)))))))),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(gamma131,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmul(gamma231,kmul(gamma331,gyzL))),ToReal(4))))),JacPDstandard33gxx),JacPDstandard11gzz))))); + + CCTK_REAL_VEC R3132 = + kmadd(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard23gxz,kadd(JacPDstandard31gyz,kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmul(gamma231,gxyL),kmadd(gamma132,kmul(gamma331,gxzL),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma333,ToReal(-2)),kmul(gamma231,kmul(gamma332,ToReal(2)))))),JacPDstandard33gxy),JacPDstandard12gzz)))))); + + CCTK_REAL_VEC R3133 = + kmul(ksub(JacPDstandard31gzz,JacPDstandard13gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3211 = + kmul(kadd(JacPDstandard12gxz,ksub(JacPDstandard13gxy,kadd(JacPDstandard31gxy,JacPDstandard21gxz))),ToReal(0.5)); + + CCTK_REAL_VEC R3212 = + kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard12gyz,kadd(JacPDstandard23gxy,kmadd(kmadd(gamma131,kmadd(gamma222,gxyL,kmul(gamma322,gxzL)),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma222,kmul(gamma331,ToReal(-2)),kmul(gamma221,kmul(gamma332,ToReal(2)))))),JacPDstandard31gyy),JacPDstandard22gxz)))),kmul(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)))))); + + CCTK_REAL_VEC R3213 = + kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmsub(ToReal(0.5),kadd(JacPDstandard12gzz,kadd(JacPDstandard33gxy,kmadd(kmadd(gamma132,kmadd(gamma231,gxyL,kmul(gamma331,gxzL)),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL))))),ToReal(2),kmul(gyzL,kmadd(gamma231,kmul(gamma332,ToReal(-2)),kmul(gamma221,kmul(gamma333,ToReal(2)))))),JacPDstandard31gyz),JacPDstandard23gxz)))),kmul(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)))))); + + CCTK_REAL_VEC R3221 = + kmadd(gamma122,kmadd(gamma131,gxxL,kmadd(gamma231,gxyL,kmul(gamma331,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gyy,kadd(JacPDstandard22gxz,kmadd(kmadd(gamma132,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmadd(gamma221,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma332,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma131,kmul(gamma222,gxyL),kmadd(gamma131,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma231,gyyL),kmadd(gamma231,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma331,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma332,ToReal(-2)),kmul(gamma222,kmul(gamma331,ToReal(2)))))),JacPDstandard32gxy),JacPDstandard21gyz)))))); + + CCTK_REAL_VEC R3222 = + kmul(ksub(JacPDstandard23gyy,JacPDstandard32gyy),ToReal(0.5)); + + CCTK_REAL_VEC R3223 = + kmadd(gamma122,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma132,kmul(kmadd(gamma232,gxyL,kmul(gamma332,gxzL)),ToReal(-2)),kmsub(ToReal(0.5),kadd(JacPDstandard22gzz,kadd(JacPDstandard33gyy,kmadd(gyyL,kmul(SQR(gamma232),ToReal(-2)),kmadd(gzzL,kmul(SQR(gamma332),ToReal(-2)),ksub(ksub(kmadd(kmadd(gamma133,kmul(gamma222,gxyL),kmadd(gamma133,kmul(gamma322,gxzL),kmadd(gamma222,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma322,gyzL),kmul(gamma322,kmul(gamma333,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma232,kmul(gamma332,ToReal(-4)),kmul(gamma222,kmul(gamma333,ToReal(2)))))),JacPDstandard32gyz),JacPDstandard23gyz))))),kmul(gxxL,SQR(gamma132))))); + + CCTK_REAL_VEC R3231 = + kmadd(gamma131,kmadd(gamma132,gxxL,kmadd(gamma232,gxyL,kmul(gamma332,gxzL))),kmul(ToReal(0.5),kadd(JacPDstandard13gyz,kadd(JacPDstandard32gxz,kmadd(kmadd(gamma133,kmadd(gamma221,gxyL,kmul(gamma321,gxzL)),kmadd(gamma121,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma221,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma321,gyzL),kmul(gamma321,kmul(gamma333,gzzL)))))),ToReal(-2),ksub(ksub(kmadd(kmadd(gamma132,kmul(gamma231,gxyL),kmadd(gamma132,kmul(gamma331,gxzL),kmadd(gamma231,kmul(gamma232,gyyL),kmadd(gamma232,kmul(gamma331,gyzL),kmul(gamma331,kmul(gamma332,gzzL)))))),ToReal(2),kmul(gyzL,kmadd(gamma221,kmul(gamma333,ToReal(-2)),kmul(gamma231,kmul(gamma332,ToReal(2)))))),JacPDstandard33gxy),JacPDstandard21gzz)))))); + + CCTK_REAL_VEC R3232 = + kmul(ToReal(0.5),kadd(JacPDstandard23gyz,kadd(JacPDstandard32gyz,kmadd(kmadd(gamma133,kmul(gamma222,gxyL),kmadd(gamma133,kmul(gamma322,gxzL),kmadd(gamma122,kmadd(gamma133,gxxL,kmadd(gamma233,gxyL,kmul(gamma333,gxzL))),kmadd(gamma222,kmul(gamma233,gyyL),kmadd(gamma233,kmul(gamma322,gyzL),kmadd(gamma222,kmul(gamma333,gyzL),kmul(gamma322,kmul(gamma333,gzzL)))))))),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(gamma132,kmadd(gamma232,gxyL,kmul(gamma332,gxzL)),kmul(gamma232,kmul(gamma332,gyzL))),ToReal(4))))),JacPDstandard33gyy),JacPDstandard22gzz))))); + + CCTK_REAL_VEC R3233 = + kmul(ksub(JacPDstandard32gzz,JacPDstandard23gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3311 = ksub(JacPDstandard13gxz,JacPDstandard31gxz); + + CCTK_REAL_VEC R3312 = + kmul(kadd(JacPDstandard13gyz,ksub(JacPDstandard23gxz,kadd(JacPDstandard32gxz,JacPDstandard31gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R3313 = + kmul(ksub(JacPDstandard13gzz,JacPDstandard31gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3321 = + kmul(kadd(JacPDstandard13gyz,ksub(JacPDstandard23gxz,kadd(JacPDstandard32gxz,JacPDstandard31gyz))),ToReal(0.5)); + + CCTK_REAL_VEC R3322 = ksub(JacPDstandard23gyz,JacPDstandard32gyz); + + CCTK_REAL_VEC R3323 = + kmul(ksub(JacPDstandard23gzz,JacPDstandard32gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3331 = + kmul(ksub(JacPDstandard13gzz,JacPDstandard31gzz),ToReal(0.5)); + + CCTK_REAL_VEC R3332 = + kmul(ksub(JacPDstandard23gzz,JacPDstandard32gzz),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 = + kadd(JacPDstandard1kxy,kmadd(gamma121,kxxL,kmadd(gamma321,kxzL,ksub(knmsub(gamma211,kyyL,kmsub(kxyL,ksub(gamma221,gamma111),kmul(gamma311,kyzL))),JacPDstandard2kxx)))); + + CCTK_REAL_VEC Ro113 = + kadd(JacPDstandard1kxz,kmadd(gamma131,kxxL,kmadd(gamma231,kxyL,ksub(knmsub(gamma211,kyzL,kmsub(kxzL,ksub(gamma331,gamma111),kmul(gamma311,kzzL))),JacPDstandard3kxx)))); + + CCTK_REAL_VEC Ro121 = + kadd(JacPDstandard2kxx,kmadd(gamma211,kyyL,kmadd(gamma311,kyzL,ksub(knmsub(gamma121,kxxL,kmsub(kxyL,ksub(gamma111,gamma221),kmul(gamma321,kxzL))),JacPDstandard1kxy)))); + + CCTK_REAL_VEC Ro122 = ToReal(0); + + CCTK_REAL_VEC Ro123 = + kadd(JacPDstandard2kxz,kmadd(gamma131,kxyL,kmadd(gamma231,kyyL,ksub(knmsub(gamma121,kxzL,kmsub(kyzL,ksub(gamma331,gamma221),kmul(gamma321,kzzL))),JacPDstandard3kxy)))); + + CCTK_REAL_VEC Ro131 = + kadd(JacPDstandard3kxx,kmadd(gamma211,kyzL,kmadd(gamma311,kzzL,ksub(knmsub(gamma131,kxxL,kmsub(kxzL,ksub(gamma111,gamma331),kmul(gamma231,kxyL))),JacPDstandard1kxz)))); + + CCTK_REAL_VEC Ro132 = + kadd(JacPDstandard3kxy,kmadd(gamma121,kxzL,kmadd(gamma321,kzzL,ksub(knmsub(gamma131,kxyL,kmsub(kyzL,ksub(gamma221,gamma331),kmul(gamma231,kyyL))),JacPDstandard2kxz)))); + + CCTK_REAL_VEC Ro133 = ToReal(0); + + CCTK_REAL_VEC Ro211 = ToReal(0); + + CCTK_REAL_VEC Ro212 = + kadd(JacPDstandard1kyy,kmadd(gamma122,kxxL,kmadd(gamma322,kxzL,ksub(knmsub(gamma221,kyyL,kmsub(kxyL,ksub(gamma222,gamma121),kmul(gamma321,kyzL))),JacPDstandard2kxy)))); + + CCTK_REAL_VEC Ro213 = + kadd(JacPDstandard1kyz,kmadd(gamma132,kxxL,kmadd(gamma232,kxyL,ksub(knmsub(gamma221,kyzL,kmsub(kxzL,ksub(gamma332,gamma121),kmul(gamma321,kzzL))),JacPDstandard3kxy)))); + + CCTK_REAL_VEC Ro221 = + kadd(JacPDstandard2kxy,kmadd(gamma221,kyyL,kmadd(gamma321,kyzL,ksub(knmsub(gamma122,kxxL,kmsub(kxyL,ksub(gamma121,gamma222),kmul(gamma322,kxzL))),JacPDstandard1kyy)))); + + CCTK_REAL_VEC Ro222 = ToReal(0); + + CCTK_REAL_VEC Ro223 = + kadd(JacPDstandard2kyz,kmadd(gamma132,kxyL,kmadd(gamma232,kyyL,ksub(knmsub(gamma122,kxzL,kmsub(kyzL,ksub(gamma332,gamma222),kmul(gamma322,kzzL))),JacPDstandard3kyy)))); + + CCTK_REAL_VEC Ro231 = + kadd(JacPDstandard3kxy,kmadd(gamma221,kyzL,kmadd(gamma321,kzzL,ksub(knmsub(gamma132,kxxL,kmsub(kxzL,ksub(gamma121,gamma332),kmul(gamma232,kxyL))),JacPDstandard1kyz)))); + + CCTK_REAL_VEC Ro232 = + kadd(JacPDstandard3kyy,kmadd(gamma122,kxzL,kmadd(gamma322,kzzL,ksub(knmsub(gamma132,kxyL,kmsub(kyzL,ksub(gamma222,gamma332),kmul(gamma232,kyyL))),JacPDstandard2kyz)))); + + CCTK_REAL_VEC Ro233 = ToReal(0); + + CCTK_REAL_VEC Ro311 = ToReal(0); + + CCTK_REAL_VEC Ro312 = + kadd(JacPDstandard1kyz,kmadd(gamma132,kxxL,kmadd(gamma332,kxzL,ksub(knmsub(gamma231,kyyL,kmsub(kxyL,ksub(gamma232,gamma131),kmul(gamma331,kyzL))),JacPDstandard2kxz)))); + + CCTK_REAL_VEC Ro313 = + kadd(JacPDstandard1kzz,kmadd(gamma133,kxxL,kmadd(gamma233,kxyL,ksub(knmsub(gamma231,kyzL,kmsub(kxzL,ksub(gamma333,gamma131),kmul(gamma331,kzzL))),JacPDstandard3kxz)))); + + CCTK_REAL_VEC Ro321 = + kadd(JacPDstandard2kxz,kmadd(gamma231,kyyL,kmadd(gamma331,kyzL,ksub(knmsub(gamma132,kxxL,kmsub(kxyL,ksub(gamma131,gamma232),kmul(gamma332,kxzL))),JacPDstandard1kyz)))); + + CCTK_REAL_VEC Ro322 = ToReal(0); + + CCTK_REAL_VEC Ro323 = + kadd(JacPDstandard2kzz,kmadd(gamma133,kxyL,kmadd(gamma233,kyyL,ksub(knmsub(gamma132,kxzL,kmsub(kyzL,ksub(gamma333,gamma232),kmul(gamma332,kzzL))),JacPDstandard3kyz)))); + + CCTK_REAL_VEC Ro331 = + kadd(JacPDstandard3kxz,kmadd(gamma231,kyzL,kmadd(gamma331,kzzL,ksub(knmsub(gamma133,kxxL,kmsub(kxzL,ksub(gamma131,gamma333),kmul(gamma233,kxyL))),JacPDstandard1kzz)))); + + CCTK_REAL_VEC Ro332 = + kadd(JacPDstandard3kyz,kmadd(gamma132,kxzL,kmadd(gamma332,kzzL,ksub(knmsub(gamma133,kxyL,kmsub(kyzL,ksub(gamma232,gamma333),kmul(gamma233,kyyL))),JacPDstandard2kzz)))); + + 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,kmul(SQR(n1),SQR(rmbar1)),kmadd(R4p2121,kmul(SQR(n2),SQR(rmbar1)),kmadd(R4p3131,kmul(SQR(n3),SQR(rmbar1)),kmadd(Rojo11,kmul(SQR(nn),SQR(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,kmul(SQR(n1),SQR(rmbar2)),kmadd(R4p2222,kmul(SQR(n2),SQR(rmbar2)),kmadd(R4p3232,kmul(SQR(n3),SQR(rmbar2)),kmadd(Rojo22,kmul(SQR(nn),SQR(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(R4p1313,kmul(SQR(n1),SQR(rmbar3)),kmadd(R4p2323,kmul(SQR(n2),SQR(rmbar3)),kmadd(R4p3333,kmul(SQR(n3),SQR(rmbar3)),kmadd(Rojo33,kmul(SQR(nn),SQR(rmbar3)),knmsub(n3,kmul(SQR(imbar3),kmadd(n1,R4p3313,kmul(n2,R4p2333))),kmadd(SQR(imbar3),kmsub(n2,kmul(nn,kmul(Ro323,ToReal(-2))),kmul(R4p2323,SQR(n2))),kmadd(SQR(imbar3),kmsub(n3,kmul(nn,kmul(Ro333,ToReal(-2))),kmul(R4p3333,SQR(n3))),kmadd(SQR(imbar3),kmsub(n1,kmul(nn,kmul(Ro313,ToReal(-2))),kmadd(n2,kmul(n3,R4p3323),kmadd(Rojo33,SQR(nn),kmul(R4p1313,SQR(n1))))),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)))))))))))))))))))))))))); + + /* If necessary, store only partial vectors after the first iteration */ + + if (CCTK_REAL_VEC_SIZE > 2 && CCTK_BUILTIN_EXPECT(i < lc_imin && i+CCTK_REAL_VEC_SIZE > lc_imax, 0)) + { + ptrdiff_t const elt_count_lo = lc_imin-i; + ptrdiff_t const elt_count_hi = lc_imax-i; + vec_store_nta_partial_mid(Psi4i[index],Psi4iL,elt_count_lo,elt_count_hi); + vec_store_nta_partial_mid(Psi4r[index],Psi4rL,elt_count_lo,elt_count_hi); + break; + } + + /* If necessary, store only partial vectors after the first iteration */ + + if (CCTK_REAL_VEC_SIZE > 1 && CCTK_BUILTIN_EXPECT(i < lc_imin, 0)) + { + ptrdiff_t const elt_count = lc_imin-i; + vec_store_nta_partial_hi(Psi4i[index],Psi4iL,elt_count); + vec_store_nta_partial_hi(Psi4r[index],Psi4rL,elt_count); + continue; + } + + /* If necessary, store only partial vectors after the last iteration */ + + if (CCTK_REAL_VEC_SIZE > 1 && CCTK_BUILTIN_EXPECT(i+CCTK_REAL_VEC_SIZE > lc_imax, 0)) + { + ptrdiff_t const elt_count = lc_imax-i; + vec_store_nta_partial_lo(Psi4i[index],Psi4iL,elt_count); + vec_store_nta_partial_lo(Psi4r[index],Psi4rL,elt_count); + break; + } + vec_store_nta(Psi4i[index],Psi4iL); + vec_store_nta(Psi4r[index],Psi4rL); + } + LC_ENDLOOP3VEC (psi4_calc_Nth); +} + +extern "C" void psi4_calc_Nth(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + + if (verbose > 1) + { + CCTK_VInfo(CCTK_THORNSTRING,"Entering psi4_calc_Nth_Body"); + } + + if (cctk_iteration % psi4_calc_Nth_calc_every != psi4_calc_Nth_calc_offset) + { + return; + } + + const char *groups[] = {"admbase::curv","admbase::metric","grid::coordinates","WeylScal4::Psi4i_group","WeylScal4::Psi4r_group"}; + GenericFD_AssertGroupStorage(cctkGH, "psi4_calc_Nth", 5, groups); + + switch(fdOrder) + { + case 2: + GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_Nth", 1, 1, 1); + break; + + case 4: + GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_Nth", 2, 2, 2); + break; + + case 6: + GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_Nth", 3, 3, 3); + break; + + case 8: + GenericFD_EnsureStencilFits(cctkGH, "psi4_calc_Nth", 4, 4, 4); + break; + } + + GenericFD_LoopOverInterior(cctkGH, &psi4_calc_Nth_Body); + + if (verbose > 1) + { + CCTK_VInfo(CCTK_THORNSTRING,"Leaving psi4_calc_Nth_Body"); + } +} |