diff options
Diffstat (limited to 'src/WeylScal4_psis_calc_Nth.cc')
-rw-r--r-- | src/WeylScal4_psis_calc_Nth.cc | 374 |
1 files changed, 193 insertions, 181 deletions
diff --git a/src/WeylScal4_psis_calc_Nth.cc b/src/WeylScal4_psis_calc_Nth.cc index ed4af44..a136d5e 100644 --- a/src/WeylScal4_psis_calc_Nth.cc +++ b/src/WeylScal4_psis_calc_Nth.cc @@ -75,69 +75,69 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int /* Include user-supplied include files */ /* Initialise finite differencing variables */ - ptrdiff_t const di CCTK_ATTRIBUTE_UNUSED = 1; - ptrdiff_t const dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); - ptrdiff_t const dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); - ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di; - ptrdiff_t const cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj; - ptrdiff_t const cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk; - CCTK_REAL_VEC const dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0)); - CCTK_REAL_VEC const dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1)); - CCTK_REAL_VEC const dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2)); - CCTK_REAL_VEC const dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME); - CCTK_REAL_VEC const t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time); - CCTK_REAL_VEC const dxi CCTK_ATTRIBUTE_UNUSED = INV(dx); - CCTK_REAL_VEC const dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); - CCTK_REAL_VEC const dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); - CCTK_REAL_VEC const khalf CCTK_ATTRIBUTE_UNUSED = ToReal(0.5); - CCTK_REAL_VEC const kthird CCTK_ATTRIBUTE_UNUSED = ToReal(1.0/3.0); - CCTK_REAL_VEC const ktwothird CCTK_ATTRIBUTE_UNUSED = ToReal(2.0/3.0); - CCTK_REAL_VEC const kfourthird CCTK_ATTRIBUTE_UNUSED = ToReal(4.0/3.0); - CCTK_REAL_VEC const keightthird CCTK_ATTRIBUTE_UNUSED = ToReal(8.0/3.0); - CCTK_REAL_VEC const hdxi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dxi); - CCTK_REAL_VEC const hdyi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dyi); - CCTK_REAL_VEC const hdzi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dzi); + ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1; + ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); + ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); + ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di; + ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj; + ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk; + CCTK_REAL_VEC /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0)); + CCTK_REAL_VEC /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1)); + CCTK_REAL_VEC /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2)); + CCTK_REAL_VEC /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME); + CCTK_REAL_VEC /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time); + CCTK_REAL_VEC /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL_VEC /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL_VEC /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL_VEC /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = ToReal(0.5); + CCTK_REAL_VEC /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = ToReal(1.0/3.0); + CCTK_REAL_VEC /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = ToReal(2.0/3.0); + CCTK_REAL_VEC /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = ToReal(4.0/3.0); + CCTK_REAL_VEC /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = ToReal(8.0/3.0); + CCTK_REAL_VEC /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dxi); + CCTK_REAL_VEC /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dyi); + CCTK_REAL_VEC /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dzi); /* Initialize predefined quantities */ - CCTK_REAL_VEC const p1o12dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dx); - CCTK_REAL_VEC const p1o12dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dy); - CCTK_REAL_VEC const p1o12dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dz); - CCTK_REAL_VEC const p1o144dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dy,dx)); - CCTK_REAL_VEC const p1o144dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dx)); - CCTK_REAL_VEC const p1o144dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dy)); - CCTK_REAL_VEC const p1o180dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dx,dx)); - CCTK_REAL_VEC const p1o180dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dy,dy)); - CCTK_REAL_VEC const p1o180dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dz,dz)); - CCTK_REAL_VEC const p1o2dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dx); - CCTK_REAL_VEC const p1o2dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dy); - CCTK_REAL_VEC const p1o2dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dz); - CCTK_REAL_VEC const p1o3600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dy,dx)); - CCTK_REAL_VEC const p1o3600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dx)); - CCTK_REAL_VEC const p1o3600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dy)); - CCTK_REAL_VEC const p1o4dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dy,dx)); - CCTK_REAL_VEC const p1o4dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dx)); - CCTK_REAL_VEC const p1o4dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dy)); - CCTK_REAL_VEC const p1o5040dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dx,dx)); - CCTK_REAL_VEC const p1o5040dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dy,dy)); - CCTK_REAL_VEC const p1o5040dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dz,dz)); - CCTK_REAL_VEC const p1o60dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dx); - CCTK_REAL_VEC const p1o60dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dy); - CCTK_REAL_VEC const p1o60dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dz); - CCTK_REAL_VEC const p1o705600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dy,dx)); - CCTK_REAL_VEC const p1o705600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dx)); - CCTK_REAL_VEC const p1o705600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dy)); - CCTK_REAL_VEC const p1o840dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dx); - CCTK_REAL_VEC const p1o840dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dy); - CCTK_REAL_VEC const p1o840dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dz); - CCTK_REAL_VEC const p1odx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dx,dx)); - CCTK_REAL_VEC const p1ody2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dy,dy)); - CCTK_REAL_VEC const p1odz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dz,dz)); - CCTK_REAL_VEC const pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx)); - CCTK_REAL_VEC const pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy)); - CCTK_REAL_VEC const pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz)); + CCTK_REAL_VEC /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dx); + CCTK_REAL_VEC /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dy); + CCTK_REAL_VEC /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dz); + CCTK_REAL_VEC /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dy,dx)); + CCTK_REAL_VEC /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dx)); + CCTK_REAL_VEC /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dy)); + CCTK_REAL_VEC /*const*/ p1o180dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dx,dx)); + CCTK_REAL_VEC /*const*/ p1o180dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dy,dy)); + CCTK_REAL_VEC /*const*/ p1o180dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dz,dz)); + CCTK_REAL_VEC /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dx); + CCTK_REAL_VEC /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dy); + CCTK_REAL_VEC /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dz); + CCTK_REAL_VEC /*const*/ p1o3600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dy,dx)); + CCTK_REAL_VEC /*const*/ p1o3600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dx)); + CCTK_REAL_VEC /*const*/ p1o3600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dy)); + CCTK_REAL_VEC /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dy,dx)); + CCTK_REAL_VEC /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dx)); + CCTK_REAL_VEC /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dy)); + CCTK_REAL_VEC /*const*/ p1o5040dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dx,dx)); + CCTK_REAL_VEC /*const*/ p1o5040dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dy,dy)); + CCTK_REAL_VEC /*const*/ p1o5040dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dz,dz)); + CCTK_REAL_VEC /*const*/ p1o60dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dx); + CCTK_REAL_VEC /*const*/ p1o60dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dy); + CCTK_REAL_VEC /*const*/ p1o60dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dz); + CCTK_REAL_VEC /*const*/ p1o705600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dy,dx)); + CCTK_REAL_VEC /*const*/ p1o705600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dx)); + CCTK_REAL_VEC /*const*/ p1o705600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dy)); + CCTK_REAL_VEC /*const*/ p1o840dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dx); + CCTK_REAL_VEC /*const*/ p1o840dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dy); + CCTK_REAL_VEC /*const*/ p1o840dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dz); + CCTK_REAL_VEC /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dx,dx)); + CCTK_REAL_VEC /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dy,dy)); + CCTK_REAL_VEC /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dz,dz)); + CCTK_REAL_VEC /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx)); + CCTK_REAL_VEC /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy)); + CCTK_REAL_VEC /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz)); /* Jacobian variable pointers */ - bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) + bool const use_jacobian CCTK_ATTRIBUTE_UNUSED = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) && strlen(jacobian_group) > 0; if (use_jacobian && strlen(jacobian_derivative_group) == 0) { @@ -191,12 +191,12 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int /* Loop over the grid points */ #pragma omp parallel - LC_LOOP3VEC(WeylScal4_psis_calc_Nth, + CCTK_LOOP3STR(WeylScal4_psis_calc_Nth, i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], cctk_ash[0],cctk_ash[1],cctk_ash[2], - CCTK_REAL_VEC_SIZE) + imin,imax, CCTK_REAL_VEC_SIZE) { - ptrdiff_t const index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; + ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; /* Assign local copies of grid functions */ @@ -790,31 +790,31 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int kmadd(J13L,PDstandard1kyz,kmadd(J23L,PDstandard2kyz,kmul(J33L,PDstandard3kyz))); JacPDstandard11gyy = - kmadd(dJ111L,PDstandard1gyy,kmadd(dJ211L,PDstandard2gyy,kmadd(dJ311L,PDstandard3gyy,kmadd(PDstandard11gyy,kmul(J11L,J11L),kmadd(PDstandard22gyy,kmul(J21L,J21L),kmadd(PDstandard33gyy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gyy,kmul(J31L,PDstandard13gyy)),kmul(J21L,kmul(J31L,PDstandard23gyy))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard1gyy,kmadd(dJ211L,PDstandard2gyy,kmadd(dJ311L,PDstandard3gyy,kmadd(PDstandard11gyy,kmul(J11L,J11L),kmadd(PDstandard22gyy,kmul(J21L,J21L),kmadd(PDstandard33gyy,kmul(J31L,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,kmul(J11L,J11L),kmadd(PDstandard22gyz,kmul(J21L,J21L),kmadd(PDstandard33gyz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gyz,kmul(J31L,PDstandard13gyz)),kmul(J21L,kmul(J31L,PDstandard23gyz))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard1gyz,kmadd(dJ211L,PDstandard2gyz,kmadd(dJ311L,PDstandard3gyz,kmadd(PDstandard11gyz,kmul(J11L,J11L),kmadd(PDstandard22gyz,kmul(J21L,J21L),kmadd(PDstandard33gyz,kmul(J31L,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,kmul(J11L,J11L),kmadd(PDstandard22gzz,kmul(J21L,J21L),kmadd(PDstandard33gzz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard12gzz,kmul(J31L,PDstandard13gzz)),kmul(J21L,kmul(J31L,PDstandard23gzz))),ToReal(2)))))))); + kmadd(dJ111L,PDstandard1gzz,kmadd(dJ211L,PDstandard2gzz,kmadd(dJ311L,PDstandard3gzz,kmadd(PDstandard11gzz,kmul(J11L,J11L),kmadd(PDstandard22gzz,kmul(J21L,J21L),kmadd(PDstandard33gzz,kmul(J31L,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,kmul(J12L,J12L),kmadd(PDstandard22gxx,kmul(J22L,J22L),kmadd(PDstandard33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx)),kmul(J22L,kmul(J32L,PDstandard23gxx))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard1gxx,kmadd(dJ222L,PDstandard2gxx,kmadd(dJ322L,PDstandard3gxx,kmadd(PDstandard11gxx,kmul(J12L,J12L),kmadd(PDstandard22gxx,kmul(J22L,J22L),kmadd(PDstandard33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxx,kmul(J32L,PDstandard13gxx)),kmul(J22L,kmul(J32L,PDstandard23gxx))),ToReal(2.)))))))); JacPDstandard22gxz = - kmadd(dJ122L,PDstandard1gxz,kmadd(dJ222L,PDstandard2gxz,kmadd(dJ322L,PDstandard3gxz,kmadd(PDstandard11gxz,kmul(J12L,J12L),kmadd(PDstandard22gxz,kmul(J22L,J22L),kmadd(PDstandard33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz)),kmul(J22L,kmul(J32L,PDstandard23gxz))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard1gxz,kmadd(dJ222L,PDstandard2gxz,kmadd(dJ322L,PDstandard3gxz,kmadd(PDstandard11gxz,kmul(J12L,J12L),kmadd(PDstandard22gxz,kmul(J22L,J22L),kmadd(PDstandard33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gxz,kmul(J32L,PDstandard13gxz)),kmul(J22L,kmul(J32L,PDstandard23gxz))),ToReal(2.)))))))); JacPDstandard22gzz = - kmadd(dJ122L,PDstandard1gzz,kmadd(dJ222L,PDstandard2gzz,kmadd(dJ322L,PDstandard3gzz,kmadd(PDstandard11gzz,kmul(J12L,J12L),kmadd(PDstandard22gzz,kmul(J22L,J22L),kmadd(PDstandard33gzz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard12gzz,kmul(J32L,PDstandard13gzz)),kmul(J22L,kmul(J32L,PDstandard23gzz))),ToReal(2)))))))); + kmadd(dJ122L,PDstandard1gzz,kmadd(dJ222L,PDstandard2gzz,kmadd(dJ322L,PDstandard3gzz,kmadd(PDstandard11gzz,kmul(J12L,J12L),kmadd(PDstandard22gzz,kmul(J22L,J22L),kmadd(PDstandard33gzz,kmul(J32L,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,kmul(J13L,J13L),kmadd(PDstandard22gxx,kmul(J23L,J23L),kmadd(PDstandard33gxx,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxx,kmul(J33L,PDstandard13gxx)),kmul(J23L,kmul(J33L,PDstandard23gxx))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard1gxx,kmadd(dJ233L,PDstandard2gxx,kmadd(dJ333L,PDstandard3gxx,kmadd(PDstandard11gxx,kmul(J13L,J13L),kmadd(PDstandard22gxx,kmul(J23L,J23L),kmadd(PDstandard33gxx,kmul(J33L,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,kmul(J13L,J13L),kmadd(PDstandard22gxy,kmul(J23L,J23L),kmadd(PDstandard33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmul(J23L,kmul(J33L,PDstandard23gxy))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard1gxy,kmadd(dJ233L,PDstandard2gxy,kmadd(dJ333L,PDstandard3gxy,kmadd(PDstandard11gxy,kmul(J13L,J13L),kmadd(PDstandard22gxy,kmul(J23L,J23L),kmadd(PDstandard33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gxy,kmul(J33L,PDstandard13gxy)),kmul(J23L,kmul(J33L,PDstandard23gxy))),ToReal(2.)))))))); JacPDstandard33gyy = - kmadd(dJ133L,PDstandard1gyy,kmadd(dJ233L,PDstandard2gyy,kmadd(dJ333L,PDstandard3gyy,kmadd(PDstandard11gyy,kmul(J13L,J13L),kmadd(PDstandard22gyy,kmul(J23L,J23L),kmadd(PDstandard33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmul(J23L,kmul(J33L,PDstandard23gyy))),ToReal(2)))))))); + kmadd(dJ133L,PDstandard1gyy,kmadd(dJ233L,PDstandard2gyy,kmadd(dJ333L,PDstandard3gyy,kmadd(PDstandard11gyy,kmul(J13L,J13L),kmadd(PDstandard22gyy,kmul(J23L,J23L),kmadd(PDstandard33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard12gyy,kmul(J33L,PDstandard13gyy)),kmul(J23L,kmul(J33L,PDstandard23gyy))),ToReal(2.)))))))); 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))))))); @@ -979,102 +979,105 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int } CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED detg = - knmsub(gyyL,kmul(gxzL,gxzL),knmsub(gxxL,kmul(gyzL,gyzL),kmadd(gzzL,kmsub(gxxL,gyyL,kmul(gxyL,gxyL)),kmul(gxyL,kmul(gxzL,kmul(gyzL,ToReal(2))))))); + kmadd(kmadd(gyyL,kmul(gxzL,gxzL),kmul(gxxL,kmul(gyzL,gyzL))),ToReal(-1.),kmadd(gzzL,kmadd(gxxL,gyyL,kmul(kmul(gxyL,gxyL),ToReal(-1.))),kmul(gxyL,kmul(gxzL,kmul(gyzL,ToReal(2.)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED invdetg = kdiv(ToReal(1),detg); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv11 = - kmul(invdetg,kmsub(gyyL,gzzL,kmul(gyzL,gyzL))); + kmul(invdetg,kmadd(gyyL,gzzL,kmul(kmul(gyzL,gyzL),ToReal(-1.)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv12 = - kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); + kmul(invdetg,kmadd(gxzL,gyzL,kmul(gxyL,kmul(gzzL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv13 = - kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); + kmul(invdetg,kmadd(gxyL,gyzL,kmul(gxzL,kmul(gyyL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv21 = - kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); + kmul(invdetg,kmadd(gxzL,gyzL,kmul(gxyL,kmul(gzzL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv22 = - kmul(invdetg,kmsub(gxxL,gzzL,kmul(gxzL,gxzL))); + kmul(invdetg,kmadd(gxxL,gzzL,kmul(kmul(gxzL,gxzL),ToReal(-1.)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv23 = - kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); + kmul(invdetg,kmadd(gxyL,gxzL,kmul(gxxL,kmul(gyzL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv31 = - kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); + kmul(invdetg,kmadd(gxyL,gyzL,kmul(gxzL,kmul(gyyL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv32 = - kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); + kmul(invdetg,kmadd(gxyL,gxzL,kmul(gxxL,kmul(gyzL,ToReal(-1.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv33 = - kmul(invdetg,kmsub(gxxL,gyyL,kmul(gxyL,gxyL))); + kmul(invdetg,kmadd(gxxL,gyyL,kmul(kmul(gxyL,gxyL),ToReal(-1.)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma111 = - kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard1gxx,knmsub(gInv12,JacPDstandard2gxx,kmsub(kmadd(gInv12,JacPDstandard1gxy,kmul(gInv13,JacPDstandard1gxz)),ToReal(2),kmul(gInv13,JacPDstandard3gxx))))); + kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard1gxx,kmadd(kmadd(gInv12,JacPDstandard2gxx,kmul(gInv13,JacPDstandard3gxx)),ToReal(-1.),kmul(kmadd(gInv12,JacPDstandard1gxy,kmul(gInv13,JacPDstandard1gxz)),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma211 = - kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard1gxx,knmsub(gInv22,JacPDstandard2gxx,kmsub(kmadd(gInv22,JacPDstandard1gxy,kmul(gInv23,JacPDstandard1gxz)),ToReal(2),kmul(gInv23,JacPDstandard3gxx))))); + kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard1gxx,kmadd(kmadd(gInv22,JacPDstandard2gxx,kmul(gInv23,JacPDstandard3gxx)),ToReal(-1.),kmul(kmadd(gInv22,JacPDstandard1gxy,kmul(gInv23,JacPDstandard1gxz)),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma311 = - kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard1gxx,knmsub(gInv32,JacPDstandard2gxx,kmsub(kmadd(gInv32,JacPDstandard1gxy,kmul(gInv33,JacPDstandard1gxz)),ToReal(2),kmul(gInv33,JacPDstandard3gxx))))); + kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard1gxx,kmadd(kmadd(gInv32,JacPDstandard2gxx,kmul(gInv33,JacPDstandard3gxx)),ToReal(-1.),kmul(kmadd(gInv32,JacPDstandard1gxy,kmul(gInv33,JacPDstandard1gxz)),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma121 = - kmul(kmadd(gInv12,JacPDstandard1gyy,kmadd(gInv11,JacPDstandard2gxx,kmul(gInv13,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv12,JacPDstandard1gyy,kmadd(gInv11,JacPDstandard2gxx,kmul(gInv13,kadd(JacPDstandard1gyz,kmadd(JacPDstandard3gxy,ToReal(-1.),JacPDstandard2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma221 = - kmul(kmadd(gInv22,JacPDstandard1gyy,kmadd(gInv21,JacPDstandard2gxx,kmul(gInv23,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv22,JacPDstandard1gyy,kmadd(gInv21,JacPDstandard2gxx,kmul(gInv23,kadd(JacPDstandard1gyz,kmadd(JacPDstandard3gxy,ToReal(-1.),JacPDstandard2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma321 = - kmul(kmadd(gInv32,JacPDstandard1gyy,kmadd(gInv31,JacPDstandard2gxx,kmul(gInv33,kadd(JacPDstandard1gyz,ksub(JacPDstandard2gxz,JacPDstandard3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv32,JacPDstandard1gyy,kmadd(gInv31,JacPDstandard2gxx,kmul(gInv33,kadd(JacPDstandard1gyz,kmadd(JacPDstandard3gxy,ToReal(-1.),JacPDstandard2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma131 = - kmul(kmadd(gInv13,JacPDstandard1gzz,kmadd(gInv11,JacPDstandard3gxx,kmul(gInv12,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv13,JacPDstandard1gzz,kmadd(gInv11,JacPDstandard3gxx,kmul(gInv12,kadd(JacPDstandard1gyz,kmadd(JacPDstandard2gxz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma231 = - kmul(kmadd(gInv23,JacPDstandard1gzz,kmadd(gInv21,JacPDstandard3gxx,kmul(gInv22,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv23,JacPDstandard1gzz,kmadd(gInv21,JacPDstandard3gxx,kmul(gInv22,kadd(JacPDstandard1gyz,kmadd(JacPDstandard2gxz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma331 = - kmul(kmadd(gInv33,JacPDstandard1gzz,kmadd(gInv31,JacPDstandard3gxx,kmul(gInv32,kadd(JacPDstandard1gyz,ksub(JacPDstandard3gxy,JacPDstandard2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv33,JacPDstandard1gzz,kmadd(gInv31,JacPDstandard3gxx,kmul(gInv32,kadd(JacPDstandard1gyz,kmadd(JacPDstandard2gxz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma122 = - kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2gyy,kmadd(gInv11,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv13,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2gyy,kmadd(gInv11,kmadd(JacPDstandard1gyy,ToReal(-1.),kmul(JacPDstandard2gxy,ToReal(2.))),kmul(gInv13,kmadd(JacPDstandard3gyy,ToReal(-1.),kmul(JacPDstandard2gyz,ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma222 = - kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2gyy,kmadd(gInv21,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv23,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2gyy,kmadd(gInv21,kmadd(JacPDstandard1gyy,ToReal(-1.),kmul(JacPDstandard2gxy,ToReal(2.))),kmul(gInv23,kmadd(JacPDstandard3gyy,ToReal(-1.),kmul(JacPDstandard2gyz,ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma322 = - kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2gyy,kmadd(gInv31,kmsub(JacPDstandard2gxy,ToReal(2),JacPDstandard1gyy),kmul(gInv33,kmsub(JacPDstandard2gyz,ToReal(2),JacPDstandard3gyy))))); + kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2gyy,kmadd(gInv31,kmadd(JacPDstandard1gyy,ToReal(-1.),kmul(JacPDstandard2gxy,ToReal(2.))),kmul(gInv33,kmadd(JacPDstandard3gyy,ToReal(-1.),kmul(JacPDstandard2gyz,ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma132 = - kmul(kmadd(gInv13,JacPDstandard2gzz,kmadd(gInv12,JacPDstandard3gyy,kmul(gInv11,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + kmul(kmadd(gInv13,JacPDstandard2gzz,kmadd(gInv12,JacPDstandard3gyy,kmul(gInv11,kadd(JacPDstandard2gxz,kmadd(JacPDstandard1gyz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma232 = - kmul(kmadd(gInv23,JacPDstandard2gzz,kmadd(gInv22,JacPDstandard3gyy,kmul(gInv21,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + kmul(kmadd(gInv23,JacPDstandard2gzz,kmadd(gInv22,JacPDstandard3gyy,kmul(gInv21,kadd(JacPDstandard2gxz,kmadd(JacPDstandard1gyz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma332 = - kmul(kmadd(gInv33,JacPDstandard2gzz,kmadd(gInv32,JacPDstandard3gyy,kmul(gInv31,kadd(JacPDstandard2gxz,ksub(JacPDstandard3gxy,JacPDstandard1gyz))))),ToReal(0.5)); + kmul(kmadd(gInv33,JacPDstandard2gzz,kmadd(gInv32,JacPDstandard3gyy,kmul(gInv31,kadd(JacPDstandard2gxz,kmadd(JacPDstandard1gyz,ToReal(-1.),JacPDstandard3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma133 = - kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard3gzz,kmadd(gInv11,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv12,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard3gzz,kmadd(kmadd(gInv11,JacPDstandard1gzz,kmul(gInv12,JacPDstandard2gzz)),ToReal(-1.),kmul(kmadd(gInv11,JacPDstandard3gxz,kmul(gInv12,JacPDstandard3gyz)),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma233 = - kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard3gzz,kmadd(gInv21,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv22,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard3gzz,kmadd(kmadd(gInv21,JacPDstandard1gzz,kmul(gInv22,JacPDstandard2gzz)),ToReal(-1.),kmul(kmadd(gInv21,JacPDstandard3gxz,kmul(gInv22,JacPDstandard3gyz)),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma333 = - kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard3gzz,kmadd(gInv31,kmsub(JacPDstandard3gxz,ToReal(2),JacPDstandard1gzz),kmul(gInv32,kmsub(JacPDstandard3gyz,ToReal(2),JacPDstandard2gzz))))); + kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard3gzz,kmadd(kmadd(gInv31,JacPDstandard1gzz,kmul(gInv32,JacPDstandard2gzz)),ToReal(-1.),kmul(kmadd(gInv31,JacPDstandard3gxz,kmul(gInv32,JacPDstandard3gyz)),ToReal(2.))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED xmoved = ksub(xL,ToReal(xorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED xmoved = + kadd(xL,ToReal(-1.*xorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ymoved = ksub(yL,ToReal(yorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ymoved = + kadd(yL,ToReal(-1.*yorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED zmoved = ksub(zL,ToReal(zorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED zmoved = + kadd(zL,ToReal(-1.*zorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va1 = kneg(ymoved); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va1 = kmul(ymoved,ToReal(-1.)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va2 = kadd(xmoved,ToReal(offset)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va3 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va3 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED vb1 = kadd(xmoved,ToReal(offset)); @@ -1083,13 +1086,13 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED vb3 = zmoved; CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED 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)))))); + kmul(ksqrt(detg),kmadd(vb3,kmadd(gInv11,va2,kmul(gInv12,kmul(va1,ToReal(-1.)))),kmadd(vb1,kmadd(gInv12,va3,kmul(gInv13,kmul(va2,ToReal(-1.)))),kmul(vb2,kmadd(gInv13,va1,kmul(gInv11,kmul(va3,ToReal(-1.)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED 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)))))); + kmul(ksqrt(detg),kmadd(vb3,kmadd(gInv21,va2,kmul(gInv22,kmul(va1,ToReal(-1.)))),kmadd(vb1,kmadd(gInv22,va3,kmul(gInv23,kmul(va2,ToReal(-1.)))),kmul(vb2,kmadd(gInv23,va1,kmul(gInv21,kmul(va3,ToReal(-1.)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED 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)))))); + kmul(ksqrt(detg),kmadd(vb3,kmadd(gInv31,va2,kmul(gInv32,kmul(va1,ToReal(-1.)))),kmadd(vb1,kmadd(gInv32,va3,kmul(gInv33,kmul(va2,ToReal(-1.)))),kmul(vb2,kmadd(gInv33,va1,kmul(gInv31,kmul(va3,ToReal(-1.)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wa1 = va1; @@ -1098,7 +1101,7 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wa3 = va3; CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega11 = - kmadd(gxxL,kmul(wa1,wa1),kmadd(gyyL,kmul(wa2,wa2),kmadd(gzzL,kmul(wa3,wa3),kmul(kmadd(gyzL,kmul(wa2,wa3),kmul(wa1,kmadd(gxyL,wa2,kmul(gxzL,wa3)))),ToReal(2))))); + kmadd(gxxL,kmul(wa1,wa1),kmadd(gyyL,kmul(wa2,wa2),kmadd(gzzL,kmul(wa3,wa3),kmul(kmadd(gyzL,kmul(wa2,wa3),kmul(wa1,kmadd(gxyL,wa2,kmul(gxzL,wa3)))),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ea1 = kdiv(wa1,ksqrt(omega11)); @@ -1109,14 +1112,17 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED 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 CCTK_ATTRIBUTE_UNUSED wb1 = knmsub(ea1,omega12,vb1); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb1 = + kmadd(ea1,kmul(omega12,ToReal(-1.)),vb1); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb2 = knmsub(ea2,omega12,vb2); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb2 = + kmadd(ea2,kmul(omega12,ToReal(-1.)),vb2); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb3 = knmsub(ea3,omega12,vb3); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb3 = + kmadd(ea3,kmul(omega12,ToReal(-1.)),vb3); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega22 = - kmadd(gxxL,kmul(wb1,wb1),kmadd(gyyL,kmul(wb2,wb2),kmadd(gzzL,kmul(wb3,wb3),kmul(kmadd(gyzL,kmul(wb2,wb3),kmul(wb1,kmadd(gxyL,wb2,kmul(gxzL,wb3)))),ToReal(2))))); + kmadd(gxxL,kmul(wb1,wb1),kmadd(gyyL,kmul(wb2,wb2),kmadd(gzzL,kmul(wb3,wb3),kmul(kmadd(gyzL,kmul(wb2,wb3),kmul(wb1,kmadd(gxyL,wb2,kmul(gxzL,wb3)))),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED eb1 = kdiv(wb1,ksqrt(omega22)); @@ -1131,16 +1137,16 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int 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 CCTK_ATTRIBUTE_UNUSED wc1 = - ksub(vc1,kmadd(eb1,omega23,kmul(ea1,omega13))); + kmadd(kmadd(ea1,omega13,kmul(eb1,omega23)),ToReal(-1.),vc1); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wc2 = - ksub(vc2,kmadd(eb2,omega23,kmul(ea2,omega13))); + kmadd(kmadd(ea2,omega13,kmul(eb2,omega23)),ToReal(-1.),vc2); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wc3 = - ksub(vc3,kmadd(eb3,omega23,kmul(ea3,omega13))); + kmadd(kmadd(ea3,omega13,kmul(eb3,omega23)),ToReal(-1.),vc3); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega33 = - kmadd(gxxL,kmul(wc1,wc1),kmadd(gyyL,kmul(wc2,wc2),kmadd(gzzL,kmul(wc3,wc3),kmul(kmadd(gyzL,kmul(wc2,wc3),kmul(wc1,kmadd(gxyL,wc2,kmul(gxzL,wc3)))),ToReal(2))))); + kmadd(gxxL,kmul(wc1,wc1),kmadd(gyyL,kmul(wc2,wc2),kmadd(gzzL,kmul(wc3,wc3),kmul(kmadd(gyzL,kmul(wc2,wc3),kmul(wc1,kmadd(gxyL,wc2,kmul(gxzL,wc3)))),ToReal(2.))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ec1 = kdiv(wc1,ksqrt(omega33)); @@ -1157,11 +1163,14 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ltet3 = kmul(eb3,isqrt2); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n1 = kneg(kmul(eb1,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n1 = + kmul(eb1,kmul(isqrt2,ToReal(-1.))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n2 = kneg(kmul(eb2,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n2 = + kmul(eb2,kmul(isqrt2,ToReal(-1.))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n3 = kneg(kmul(eb3,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n3 = + kmul(eb3,kmul(isqrt2,ToReal(-1.))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rm1 = kmul(ec1,isqrt2); @@ -1181,178 +1190,181 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rmbar3 = kmul(ec3,isqrt2); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar1 = kneg(kmul(ea1,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar1 = + kmul(ea1,kmul(isqrt2,ToReal(-1.))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar2 = kneg(kmul(ea2,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar2 = + kmul(ea2,kmul(isqrt2,ToReal(-1.))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar3 = kneg(kmul(ea3,isqrt2)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar3 = + kmul(ea3,kmul(isqrt2,ToReal(-1.))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED nn = isqrt2; CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1212 = - kmul(ToReal(0.5),kadd(JacPDstandard12gxy,kadd(JacPDstandard21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma121,gamma121),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4))))),JacPDstandard22gxx),JacPDstandard11gyy))))); + kmul(ToReal(0.5),kadd(JacPDstandard12gxy,kadd(JacPDstandard21gxy,kmadd(kmadd(gamma122,kmadd(gxxL,gamma111,kmadd(gxyL,gamma211,kmul(gxzL,gamma311))),kmadd(gamma222,kmadd(gxyL,gamma111,kmadd(gyyL,gamma211,kmul(gyzL,gamma311))),kmul(kmadd(gxzL,gamma111,kmadd(gyzL,gamma211,kmul(gzzL,gamma311))),gamma322))),ToReal(-2.),kmadd(kadd(JacPDstandard11gyy,JacPDstandard22gxx),ToReal(-1.),kmul(kmadd(gamma121,kmadd(gxxL,gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321))),kmadd(gamma221,kmadd(gxyL,gamma121,kmadd(gyyL,gamma221,kmul(gyzL,gamma321))),kmul(gamma321,kmadd(gxzL,gamma121,kmadd(gyzL,gamma221,kmul(gzzL,gamma321)))))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1213 = - kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard12gxz,kadd(JacPDstandard31gxy,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2),JacPDstandard23gxx),JacPDstandard11gyz)))))); + kmul(ToReal(0.5),kadd(JacPDstandard12gxz,kadd(JacPDstandard31gxy,kmadd(kmadd(gamma132,kmadd(gxxL,gamma111,kmadd(gxyL,gamma211,kmul(gxzL,gamma311))),kmadd(gamma232,kmadd(gxyL,gamma111,kmadd(gyyL,gamma211,kmul(gyzL,gamma311))),kmul(kmadd(gxzL,gamma111,kmadd(gyzL,gamma211,kmul(gzzL,gamma311))),gamma332))),ToReal(-2.),kmadd(kadd(JacPDstandard11gyz,JacPDstandard23gxx),ToReal(-1.),kmul(kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmadd(gamma221,kmadd(gxyL,gamma131,kmadd(gyyL,gamma231,kmul(gyzL,gamma331))),kmul(gamma321,kmadd(gxzL,gamma131,kmadd(gyzL,gamma231,kmul(gzzL,gamma331)))))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1223 = - kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard22gxz,kadd(JacPDstandard31gyy,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2),JacPDstandard23gxy),JacPDstandard12gyz)))))); + kmul(ToReal(0.5),kadd(JacPDstandard22gxz,kadd(JacPDstandard31gyy,kmadd(kmadd(gamma132,kmadd(gxxL,gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321))),kmadd(gamma232,kmadd(gxyL,gamma121,kmadd(gyyL,gamma221,kmul(gyzL,gamma321))),kmul(kmadd(gxzL,gamma121,kmadd(gyzL,gamma221,kmul(gzzL,gamma321))),gamma332))),ToReal(-2.),kmadd(kadd(JacPDstandard12gyz,JacPDstandard23gxy),ToReal(-1.),kmul(kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmadd(gamma222,kmadd(gxyL,gamma131,kmadd(gyyL,gamma231,kmul(gyzL,gamma331))),kmul(gamma322,kmadd(gxzL,gamma131,kmadd(gyzL,gamma231,kmul(gzzL,gamma331)))))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1313 = - kmul(ToReal(0.5),kadd(JacPDstandard13gxz,kadd(JacPDstandard31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma131,gamma131),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4))))),JacPDstandard33gxx),JacPDstandard11gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard13gxz,kadd(JacPDstandard31gxz,kmadd(kmadd(gamma133,kmadd(gxxL,gamma111,kmadd(gxyL,gamma211,kmul(gxzL,gamma311))),kmadd(gamma233,kmadd(gxyL,gamma111,kmadd(gyyL,gamma211,kmul(gyzL,gamma311))),kmul(kmadd(gxzL,gamma111,kmadd(gyzL,gamma211,kmul(gzzL,gamma311))),gamma333))),ToReal(-2.),kmadd(kadd(JacPDstandard11gzz,JacPDstandard33gxx),ToReal(-1.),kmul(kmadd(gamma131,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmadd(gamma231,kmadd(gxyL,gamma131,kmadd(gyyL,gamma231,kmul(gyzL,gamma331))),kmul(gamma331,kmadd(gxzL,gamma131,kmadd(gyzL,gamma231,kmul(gzzL,gamma331)))))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1323 = - kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard23gxz,kadd(JacPDstandard31gyz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2),JacPDstandard33gxy),JacPDstandard12gzz)))))); + kmul(ToReal(0.5),kadd(JacPDstandard23gxz,kadd(JacPDstandard31gyz,kmadd(kmadd(gamma133,kmadd(gxxL,gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321))),kmadd(gamma233,kmadd(gxyL,gamma121,kmadd(gyyL,gamma221,kmul(gyzL,gamma321))),kmul(kmadd(gxzL,gamma121,kmadd(gyzL,gamma221,kmul(gzzL,gamma321))),gamma333))),ToReal(-2.),kmadd(kadd(JacPDstandard12gzz,JacPDstandard33gxy),ToReal(-1.),kmul(kmadd(gamma132,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmadd(gamma232,kmadd(gxyL,gamma131,kmadd(gyyL,gamma231,kmul(gyzL,gamma331))),kmul(kmadd(gxzL,gamma131,kmadd(gyzL,gamma231,kmul(gzzL,gamma331))),gamma332))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R2323 = - kmul(ToReal(0.5),kadd(JacPDstandard23gyz,kadd(JacPDstandard32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma132,gamma132),ToReal(2)),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(2)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(2)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4))))),JacPDstandard33gyy),JacPDstandard22gzz))))); + kmul(ToReal(0.5),kadd(JacPDstandard23gyz,kadd(JacPDstandard32gyz,kmadd(kmadd(gamma133,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma233,kmadd(gxyL,gamma122,kmadd(gyyL,gamma222,kmul(gyzL,gamma322))),kmul(kmadd(gxzL,gamma122,kmadd(gyzL,gamma222,kmul(gzzL,gamma322))),gamma333))),ToReal(-2.),kmadd(kadd(JacPDstandard22gzz,JacPDstandard33gyy),ToReal(-1.),kmul(kmadd(gamma132,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmadd(gamma232,kmadd(gxyL,gamma132,kmadd(gyyL,gamma232,kmul(gyzL,gamma332))),kmul(gamma332,kmadd(gxzL,gamma132,kmadd(gyzL,gamma232,kmul(gzzL,gamma332)))))),ToReal(2.))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1212 = - kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)); + kmadd(kxxL,kyyL,kmadd(kmul(kxyL,kxyL),ToReal(-1.),R1212)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1213 = - kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)); + kmadd(kxxL,kyzL,kmadd(kxyL,kmul(kxzL,ToReal(-1.)),R1213)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1223 = - kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1223)); + kmadd(kxyL,kyzL,kmadd(kxzL,kmul(kyyL,ToReal(-1.)),R1223)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1313 = - kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)); + kmadd(kxxL,kzzL,kmadd(kmul(kxzL,kxzL),ToReal(-1.),R1313)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1323 = - kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)); + kmadd(kxyL,kzzL,kmadd(kxzL,kmul(kyzL,ToReal(-1.)),R1323)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p2323 = - kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323)); + kmadd(kyyL,kzzL,kmadd(kmul(kyzL,kyzL),ToReal(-1.),R2323)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro111 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro111 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro112 = - kmadd(kxxL,gamma121,kmadd(kxzL,gamma321,kadd(JacPDstandard1kxy,knmsub(kyyL,gamma211,knmsub(kyzL,gamma311,kmsub(kxyL,ksub(gamma221,gamma111),JacPDstandard2kxx)))))); + kmadd(kxxL,gamma121,kmadd(kxyL,gamma221,kmadd(kxzL,gamma321,kmadd(kmadd(kxyL,gamma111,kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,JacPDstandard2kxx))),ToReal(-1.),JacPDstandard1kxy)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro113 = - kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,kadd(JacPDstandard1kxz,knmsub(kyzL,gamma211,knmsub(kzzL,gamma311,kmsub(kxzL,ksub(gamma331,gamma111),JacPDstandard3kxx)))))); + kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,kmadd(kxzL,gamma331,kmadd(kmadd(kxzL,gamma111,kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,JacPDstandard3kxx))),ToReal(-1.),JacPDstandard1kxz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro121 = - kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,kadd(JacPDstandard2kxx,knmsub(kxxL,gamma121,knmsub(kxzL,gamma321,kmsub(kxyL,ksub(gamma111,gamma221),JacPDstandard1kxy)))))); + kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,kadd(JacPDstandard2kxx,kmadd(kmadd(kxxL,gamma121,kmadd(kxzL,gamma321,JacPDstandard1kxy)),ToReal(-1.),kmul(kxyL,kmadd(gamma221,ToReal(-1.),gamma111)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro122 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro122 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro123 = - kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kadd(JacPDstandard2kxz,knmsub(kxzL,gamma121,knmsub(kzzL,gamma321,kmsub(kyzL,ksub(gamma331,gamma221),JacPDstandard3kxy)))))); + kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kmadd(kmadd(kxzL,gamma121,kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,JacPDstandard3kxy))),ToReal(-1.),JacPDstandard2kxz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro131 = - kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,kadd(JacPDstandard3kxx,knmsub(kxxL,gamma131,knmsub(kxyL,gamma231,kmsub(kxzL,ksub(gamma111,gamma331),JacPDstandard1kxz)))))); + kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,kadd(JacPDstandard3kxx,kmadd(kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,JacPDstandard1kxz)),ToReal(-1.),kmul(kxzL,kmadd(gamma331,ToReal(-1.),gamma111)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro132 = - kmadd(kxzL,gamma121,kmadd(kzzL,gamma321,kadd(JacPDstandard3kxy,knmsub(kxyL,gamma131,knmsub(kyyL,gamma231,kmsub(kyzL,ksub(gamma221,gamma331),JacPDstandard2kxz)))))); + kmadd(kxzL,gamma121,kmadd(kzzL,gamma321,kadd(JacPDstandard3kxy,kmadd(kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,JacPDstandard2kxz)),ToReal(-1.),kmul(kyzL,kmadd(gamma331,ToReal(-1.),gamma221)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro133 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro133 = ToReal(0.); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro211 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro211 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro212 = - kmadd(kxxL,gamma122,kmadd(kxzL,gamma322,kadd(JacPDstandard1kyy,knmsub(kyyL,gamma221,knmsub(kyzL,gamma321,kmsub(kxyL,ksub(gamma222,gamma121),JacPDstandard2kxy)))))); + kmadd(kxxL,gamma122,kmadd(kxyL,gamma222,kmadd(kxzL,gamma322,kmadd(kmadd(kxyL,gamma121,kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,JacPDstandard2kxy))),ToReal(-1.),JacPDstandard1kyy)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro213 = - kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kadd(JacPDstandard1kyz,knmsub(kyzL,gamma221,knmsub(kzzL,gamma321,kmsub(kxzL,ksub(gamma332,gamma121),JacPDstandard3kxy)))))); + kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kmadd(kxzL,gamma332,kmadd(kmadd(kxzL,gamma121,kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,JacPDstandard3kxy))),ToReal(-1.),JacPDstandard1kyz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro221 = - kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,kadd(JacPDstandard2kxy,knmsub(kxxL,gamma122,knmsub(kxzL,gamma322,kmsub(kxyL,ksub(gamma121,gamma222),JacPDstandard1kyy)))))); + kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,kadd(JacPDstandard2kxy,kmadd(kmadd(kxxL,gamma122,kmadd(kxzL,gamma322,JacPDstandard1kyy)),ToReal(-1.),kmul(kxyL,kmadd(gamma222,ToReal(-1.),gamma121)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro222 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro222 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro223 = - kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,kadd(JacPDstandard2kyz,knmsub(kxzL,gamma122,knmsub(kzzL,gamma322,kmsub(kyzL,ksub(gamma332,gamma222),JacPDstandard3kyy)))))); + kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,kmadd(kyzL,gamma332,kmadd(kmadd(kxzL,gamma122,kmadd(kyzL,gamma222,kmadd(kzzL,gamma322,JacPDstandard3kyy))),ToReal(-1.),JacPDstandard2kyz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro231 = - kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,kadd(JacPDstandard3kxy,knmsub(kxxL,gamma132,knmsub(kxyL,gamma232,kmsub(kxzL,ksub(gamma121,gamma332),JacPDstandard1kyz)))))); + kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,kadd(JacPDstandard3kxy,kmadd(kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,JacPDstandard1kyz)),ToReal(-1.),kmul(kxzL,kmadd(gamma332,ToReal(-1.),gamma121)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro232 = - kmadd(kxzL,gamma122,kmadd(kzzL,gamma322,kadd(JacPDstandard3kyy,knmsub(kxyL,gamma132,knmsub(kyyL,gamma232,kmsub(kyzL,ksub(gamma222,gamma332),JacPDstandard2kyz)))))); + kmadd(kxzL,gamma122,kmadd(kzzL,gamma322,kadd(JacPDstandard3kyy,kmadd(kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,JacPDstandard2kyz)),ToReal(-1.),kmul(kyzL,kmadd(gamma332,ToReal(-1.),gamma222)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro233 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro233 = ToReal(0.); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro311 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro311 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro312 = - kmadd(kxxL,gamma132,kmadd(kxzL,gamma332,kadd(JacPDstandard1kyz,knmsub(kyyL,gamma231,knmsub(kyzL,gamma331,kmsub(kxyL,ksub(gamma232,gamma131),JacPDstandard2kxz)))))); + kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kmadd(kxzL,gamma332,kmadd(kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,JacPDstandard2kxz))),ToReal(-1.),JacPDstandard1kyz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro313 = - kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,kadd(JacPDstandard1kzz,knmsub(kyzL,gamma231,knmsub(kzzL,gamma331,kmsub(kxzL,ksub(gamma333,gamma131),JacPDstandard3kxz)))))); + kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,kmadd(kxzL,gamma333,kmadd(kmadd(kxzL,gamma131,kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,JacPDstandard3kxz))),ToReal(-1.),JacPDstandard1kzz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro321 = - kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kadd(JacPDstandard2kxz,knmsub(kxxL,gamma132,knmsub(kxzL,gamma332,kmsub(kxyL,ksub(gamma131,gamma232),JacPDstandard1kyz)))))); + kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kadd(JacPDstandard2kxz,kmadd(kmadd(kxxL,gamma132,kmadd(kxzL,gamma332,JacPDstandard1kyz)),ToReal(-1.),kmul(kxyL,kmadd(gamma232,ToReal(-1.),gamma131)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro322 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro322 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro323 = - kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,kadd(JacPDstandard2kzz,knmsub(kxzL,gamma132,knmsub(kzzL,gamma332,kmsub(kyzL,ksub(gamma333,gamma232),JacPDstandard3kyz)))))); + kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,kmadd(kyzL,gamma333,kmadd(kmadd(kxzL,gamma132,kmadd(kyzL,gamma232,kmadd(kzzL,gamma332,JacPDstandard3kyz))),ToReal(-1.),JacPDstandard2kzz)))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro331 = - kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,kadd(JacPDstandard3kxz,knmsub(kxxL,gamma133,knmsub(kxyL,gamma233,kmsub(kxzL,ksub(gamma131,gamma333),JacPDstandard1kzz)))))); + kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,kadd(JacPDstandard3kxz,kmadd(kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,JacPDstandard1kzz)),ToReal(-1.),kmul(kxzL,kmadd(gamma333,ToReal(-1.),gamma131)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro332 = - kmadd(kxzL,gamma132,kmadd(kzzL,gamma332,kadd(JacPDstandard3kyz,knmsub(kxyL,gamma133,knmsub(kyyL,gamma233,kmsub(kyzL,ksub(gamma232,gamma333),JacPDstandard2kzz)))))); + kmadd(kxzL,gamma132,kmadd(kzzL,gamma332,kadd(JacPDstandard3kyz,kmadd(kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,JacPDstandard2kzz)),ToReal(-1.),kmul(kyzL,kmadd(gamma333,ToReal(-1.),gamma232)))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro333 = ToReal(0); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro333 = ToReal(0.); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo11 = - kmadd(kadd(gInv23,gInv32),kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv22,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313))))); + kmadd(kadd(gInv23,gInv32),kmadd(kxxL,kyzL,kmadd(kxyL,kmul(kxzL,ToReal(-1.)),R1213)),kmadd(gInv22,kmadd(kxxL,kyyL,kmadd(kmul(kxyL,kxyL),ToReal(-1.),R1212)),kmul(gInv33,kmadd(kxxL,kzzL,kmadd(kmul(kxzL,kxzL),ToReal(-1.),R1313))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo12 = - kmadd(gInv23,R1223,knmsub(gInv21,R1212,knmsub(gInv31,R1213,kmadd(gInv12,kmsub(kxyL,kxyL,kmul(kxxL,kyyL)),kmadd(gInv32,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv13,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))); + kmadd(gInv23,R1223,kmadd(kmadd(gInv21,R1212,kmul(gInv31,R1213)),ToReal(-1.),kmadd(gInv12,kmadd(kxyL,kxyL,kmul(kxxL,kmul(kyyL,ToReal(-1.)))),kmadd(gInv32,kmadd(kxyL,kyzL,kmul(kxzL,kmul(kyyL,ToReal(-1.)))),kmadd(gInv13,kmadd(kxyL,kxzL,kmul(kxxL,kmul(kyzL,ToReal(-1.)))),kmul(gInv33,kmadd(kxyL,kzzL,kmadd(kxzL,kmul(kyzL,ToReal(-1.)),R1323)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo13 = - kmadd(kxzL,kmul(kyzL,gInv23),kmadd(gInv13,kmul(kxzL,kxzL),knmsub(gInv21,R1213,knmsub(gInv31,R1313,knmsub(gInv32,R1323,kmadd(gInv12,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmsub(gInv22,kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmul(kzzL,kmadd(kxyL,gInv23,kmul(kxxL,gInv13)))))))))); + kmadd(kxzL,kmadd(kyyL,gInv22,kmul(kyzL,gInv23)),kmadd(gInv13,kmul(kxzL,kxzL),kmadd(kmadd(kzzL,kmadd(kxxL,gInv13,kmul(kxyL,gInv23)),kmadd(gInv21,R1213,kmadd(gInv22,kmadd(kxyL,kyzL,R1223),kmadd(gInv31,R1313,kmul(gInv32,R1323))))),ToReal(-1.),kmul(gInv12,kmadd(kxyL,kxzL,kmul(kxxL,kmul(kyzL,ToReal(-1.)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo21 = - kmadd(gInv32,R1223,knmsub(gInv12,R1212,knmsub(gInv13,R1213,kmadd(gInv21,kmsub(kxyL,kxyL,kmul(kxxL,kyyL)),kmadd(gInv23,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv31,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))); + kmadd(gInv32,R1223,kmadd(kmadd(gInv12,R1212,kmul(gInv13,R1213)),ToReal(-1.),kmadd(gInv21,kmadd(kxyL,kxyL,kmul(kxxL,kmul(kyyL,ToReal(-1.)))),kmadd(gInv23,kmadd(kxyL,kyzL,kmul(kxzL,kmul(kyyL,ToReal(-1.)))),kmadd(gInv31,kmadd(kxyL,kxzL,kmul(kxxL,kmul(kyzL,ToReal(-1.)))),kmul(gInv33,kmadd(kxyL,kzzL,kmadd(kxzL,kmul(kyzL,ToReal(-1.)),R1323)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo22 = - kmadd(kadd(gInv13,gInv31),kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmadd(gInv11,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))); + kmadd(kadd(gInv13,gInv31),kmadd(kxzL,kyyL,kmul(kmadd(kxyL,kyzL,R1223),ToReal(-1.))),kmadd(gInv11,kmadd(kxxL,kyyL,kmadd(kmul(kxyL,kxyL),ToReal(-1.),R1212)),kmul(gInv33,kmadd(kyyL,kzzL,kmadd(kmul(kyzL,kyzL),ToReal(-1.),R2323))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo23 = - kmadd(gInv12,R1223,knmsub(gInv31,R1323,knmsub(gInv32,R2323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv21,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv13,kmsub(kxzL,kyzL,kmul(kxyL,kzzL)),kmul(gInv23,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))))))))); + kmadd(gInv12,R1223,kmadd(gInv23,kmul(kyzL,kyzL),kmadd(kmadd(kyyL,kmul(kzzL,gInv23),kmadd(gInv31,R1323,kmul(gInv32,R2323))),ToReal(-1.),kmadd(gInv11,kmadd(kxxL,kyzL,kmadd(kxyL,kmul(kxzL,ToReal(-1.)),R1213)),kmadd(gInv21,kmadd(kxyL,kyzL,kmul(kxzL,kmul(kyyL,ToReal(-1.)))),kmul(gInv13,kmadd(kxzL,kyzL,kmul(kxyL,kmul(kzzL,ToReal(-1.)))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo31 = - kmadd(kxzL,kmul(kyzL,gInv32),kmadd(gInv31,kmul(kxzL,kxzL),knmsub(gInv12,R1213,knmsub(gInv13,R1313,knmsub(gInv23,R1323,kmadd(gInv21,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmsub(gInv22,kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmul(kzzL,kmadd(kxyL,gInv32,kmul(kxxL,gInv31)))))))))); + kmadd(kxzL,kmul(kyzL,gInv32),kmadd(gInv31,kmul(kxzL,kxzL),kmadd(kmadd(kzzL,kmadd(kxxL,gInv31,kmul(kxyL,gInv32)),kmadd(gInv12,R1213,kmadd(gInv13,R1313,kmul(gInv23,R1323)))),ToReal(-1.),kmadd(gInv21,kmadd(kxyL,kxzL,kmul(kxxL,kmul(kyzL,ToReal(-1.)))),kmul(gInv22,kmadd(kxzL,kyyL,kmul(kmadd(kxyL,kyzL,R1223),ToReal(-1.)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo32 = - kmadd(gInv21,R1223,knmsub(gInv13,R1323,knmsub(gInv23,R2323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv12,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv31,kmsub(kxzL,kyzL,kmul(kxyL,kzzL)),kmul(gInv32,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))))))))); + kmadd(gInv21,R1223,kmadd(gInv32,kmul(kyzL,kyzL),kmadd(kmadd(kyyL,kmul(kzzL,gInv32),kmadd(gInv13,R1323,kmul(gInv23,R2323))),ToReal(-1.),kmadd(gInv11,kmadd(kxxL,kyzL,kmadd(kxyL,kmul(kxzL,ToReal(-1.)),R1213)),kmadd(gInv12,kmadd(kxyL,kyzL,kmul(kxzL,kmul(kyyL,ToReal(-1.)))),kmul(gInv31,kmadd(kxzL,kyzL,kmul(kxyL,kmul(kzzL,ToReal(-1.)))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo33 = - kmadd(kadd(gInv12,gInv21),kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)),kmadd(gInv11,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)),kmul(gInv22,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))); + kmadd(kadd(gInv12,gInv21),kmadd(kxyL,kzzL,kmadd(kxzL,kmul(kyzL,ToReal(-1.)),R1323)),kmadd(gInv11,kmadd(kxxL,kzzL,kmadd(kmul(kxzL,kxzL),ToReal(-1.),R1313)),kmul(gInv22,kmadd(kyyL,kzzL,kmadd(kmul(kyzL,kyzL),ToReal(-1.),R2323))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi4rL = - kmadd(rmbar1,kmul(rmbar2,kmul(Rojo12,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo13,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo21,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo23,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo31,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo32,kmul(nn,nn))),kmadd(R4p1212,kmul(kmul(n2,n2),kmul(rmbar1,rmbar1)),kmadd(Rojo11,kmul(kmul(nn,nn),kmul(rmbar1,rmbar1)),kmadd(R4p1212,kmul(kmul(n1,n1),kmul(rmbar2,rmbar2)),kmadd(Rojo22,kmul(kmul(nn,nn),kmul(rmbar2,rmbar2)),kmadd(R4p1313,kmul(kmul(n1,n1),kmul(rmbar3,rmbar3)),kmadd(R4p2323,kmul(kmul(n2,n2),kmul(rmbar3,rmbar3)),kmadd(kmadd(n1,kmul(n2,kmul(R4p1212,kmul(rmbar1,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1213,kmul(rmbar1,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1313,kmul(rmbar1,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p1323,kmul(rmbar2,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p2323,kmul(rmbar2,rmbar3))),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(imbar3,imbar3))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(imbar3,imbar3))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(imbar3,imbar3))),kmul(R4p1223,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2)))))))))))),ToReal(-2),kmadd(n1,kmul(n3,kmul(R4p1223,kmul(kmul(rmbar2,rmbar2),ToReal(-2)))),kmadd(R4p1313,kmsub(kmul(n3,n3),kmul(rmbar1,rmbar1),kmul(kmul(imbar3,imbar3),kmul(n1,n1))),kmadd(R4p2323,kmsub(kmul(n3,n3),kmul(rmbar2,rmbar2),kmul(kmul(imbar3,imbar3),kmul(n2,n2))),kmadd(Rojo33,kmsub(kmul(nn,nn),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(nn,nn))),kmadd(kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro112))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro113))),kmadd(n2,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro123))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro132))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro133))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro211))),kmadd(n1,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro213))),kmadd(n2,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro223))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro231))),kmadd(n3,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro233))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro311))),kmadd(n1,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro312))),kmadd(n2,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro322))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro331))),kmadd(n3,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro332))),kmadd(R4p1213,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmul(R4p1323,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3)))))))))))))))))))),ToReal(2),kmadd(n1,kmul(nn,kmul(Ro111,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro121,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro131,kmul(kmul(rmbar1,rmbar1),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro212,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro222,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro232,kmul(kmul(rmbar2,rmbar2),ToReal(2)))),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(rmbar3,rmbar3),ToReal(2)))),knmsub(kmul(imbar2,imbar2),kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmadd(Rojo22,kmul(nn,nn),kmadd(n1,kmul(n3,kmul(R4p1223,ToReal(-2))),kmadd(n1,kmul(nn,kmul(Ro212,ToReal(2))),kmadd(n2,kmul(nn,kmul(Ro222,ToReal(2))),kmul(n3,kmul(nn,kmul(Ro232,ToReal(2)))))))))),kmadd(n2,kmadd(n1,kmul(R4p1213,kmul(rmbar1,kmul(rmbar3,ToReal(-2)))),kmadd(n3,kmul(R4p1323,kmul(rmbar1,kmul(rmbar3,ToReal(-2)))),kmadd(nn,kmul(Ro323,kmul(kmul(imbar3,imbar3),ToReal(-2))),kmadd(n3,kmul(R4p1223,kmul(rmbar1,kmul(rmbar2,ToReal(2)))),kmadd(n1,kmul(R4p1223,kmul(rmbar2,kmul(rmbar3,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro122,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro221,ToReal(2)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro321,ToReal(2)))),kmul(n3,kmul(R4p1213,kmul(kmul(rmbar1,rmbar1),ToReal(2)))))))))))),knmsub(kmul(imbar1,imbar1),kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmadd(n2,kmul(kmadd(n3,R4p1213,kmul(nn,Ro121)),ToReal(2)),kmadd(n3,kmul(nn,kmul(Ro131,ToReal(2))),kmul(nn,kmadd(nn,Rojo11,kmul(n1,kmul(Ro111,ToReal(2))))))))),kmsub(imbar1,kmul(imbar3,kmadd(R4p1223,kmul(kmul(n2,n2),ToReal(2)),kmadd(n1,kmul(kmadd(n2,R4p1213,kmsub(n3,R4p1313,kmul(nn,kadd(Ro113,Ro311)))),ToReal(2)),kmsub(n2,kmul(kmsub(n3,R4p1323,kmul(nn,kadd(Ro123,Ro321))),ToReal(2)),kmul(nn,kmadd(nn,kadd(Rojo13,Rojo31),kmul(n3,kmul(kadd(Ro133,Ro331),ToReal(2))))))))),kmul(imbar2,kmadd(imbar3,kmadd(R4p1213,kmul(kmul(n1,n1),ToReal(2)),kmadd(n1,kmul(kmadd(n2,R4p1223,kmsub(nn,kadd(Ro213,Ro312),kmul(n3,R4p1323))),ToReal(2)),kmadd(n2,kmadd(n3,kmul(R4p2323,ToReal(-2)),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2)))),kmul(nn,kmadd(nn,kadd(Rojo23,Rojo32),kmul(n3,kmul(kadd(Ro233,Ro332),ToReal(2)))))))),kmul(imbar1,kmadd(Rojo12,kmul(nn,nn),kmadd(Rojo21,kmul(nn,nn),kmadd(n1,kmul(ToReal(-2),kmadd(n2,R4p1212,kmsub(n3,R4p1213,kmul(nn,kadd(Ro112,Ro211))))),kmadd(n3,kmul(nn,kmul(Ro132,ToReal(2))),kmadd(n2,kmul(kmadd(n3,R4p1223,kmul(nn,kadd(Ro122,Ro221))),ToReal(2)),kmadd(n3,kmul(nn,kmul(Ro231,ToReal(2))),kmul(R4p1323,kmul(kmul(n3,n3),ToReal(2)))))))))))))))))))))))))))))))))))))))))))); + kmadd(kmadd(kmadd(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),ToReal(-1.)))),kmadd(imbar1,imbar2,kmul(rmbar1,kmul(rmbar2,ToReal(-1.)))),kmadd(kmadd(n1,kmadd(n2,R4p1213,kmul(n3,R4p1313)),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2)))),kmadd(imbar1,imbar3,kmul(rmbar1,kmul(rmbar3,ToReal(-1.)))),kmadd(kmadd(rmbar2,rmbar3,kmul(imbar2,kmul(imbar3,ToReal(-1.)))),kmadd(R4p1213,kmul(n1,n1),kmadd(n2,kmul(n3,kmul(R4p2323,ToReal(-1.))),kmul(n1,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))),kmul(nn,kmadd(kmadd(n1,kadd(Ro112,Ro211),kmadd(n2,kadd(Ro122,Ro221),kmul(n3,kadd(Ro132,Ro231)))),kmadd(rmbar1,rmbar2,kmul(imbar1,kmul(imbar2,ToReal(-1.)))),kmadd(kmadd(n1,kadd(Ro113,Ro311),kmadd(n2,kadd(Ro123,Ro321),kmul(n3,kadd(Ro133,Ro331)))),kmadd(rmbar1,rmbar3,kmul(imbar1,kmul(imbar3,ToReal(-1.)))),kmadd(kmadd(n1,kadd(Ro213,Ro312),kmadd(n2,kadd(Ro223,Ro322),kmul(n3,kadd(Ro233,Ro332)))),kmadd(rmbar2,rmbar3,kmul(imbar2,kmul(imbar3,ToReal(-1.)))),kmadd(kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131))),kmadd(rmbar1,rmbar1,kmul(kmul(imbar1,imbar1),ToReal(-1.))),kmadd(kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232))),kmadd(rmbar2,rmbar2,kmul(kmul(imbar2,imbar2),ToReal(-1.))),kmul(kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))),kmadd(rmbar3,rmbar3,kmul(kmul(imbar3,imbar3),ToReal(-1.))))))))))))),ToReal(2.),kmul(ToReal(-1.),kmadd(kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2.)))))),kmadd(imbar2,imbar2,kmul(kmul(rmbar2,rmbar2),ToReal(-1.))),kmadd(kmul(nn,nn),kmadd(imbar1,kmadd(imbar2,kadd(Rojo12,Rojo21),kmul(imbar3,kadd(Rojo13,Rojo31))),kmadd(imbar2,kmul(imbar3,Rojo32),kmadd(Rojo11,kmul(imbar1,imbar1),kmadd(kmadd(rmbar1,kmadd(rmbar2,kadd(Rojo12,Rojo21),kmul(rmbar3,kadd(Rojo13,Rojo31))),kmadd(rmbar2,kmul(rmbar3,Rojo32),kmul(Rojo11,kmul(rmbar1,rmbar1)))),ToReal(-1.),kmadd(Rojo23,kmadd(imbar2,imbar3,kmul(rmbar2,kmul(rmbar3,ToReal(-1.)))),kmadd(Rojo22,kmadd(imbar2,imbar2,kmul(kmul(rmbar2,rmbar2),ToReal(-1.))),kmul(Rojo33,kmadd(imbar3,imbar3,kmul(kmul(rmbar3,rmbar3),ToReal(-1.)))))))))),kmadd(kmadd(imbar1,imbar1,kmul(kmul(rmbar1,rmbar1),ToReal(-1.))),kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2.)))))),kmul(kmadd(imbar3,imbar3,kmul(kmul(rmbar3,rmbar3),ToReal(-1.))),kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2.)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi4iL = - knmadd(im1,kmadd(kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmul(kmul(n3,n3),ToReal(2)),kmadd(kmul(n2,n2),kmul(kmsub(R4p1212,rm1,kmul(R4p1223,rm3)),ToReal(2)),kmadd(n2,kmul(ToReal(2),knmsub(n1,kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmadd(n3,kmadd(R4p1223,rm2,kmsub(R4p1213,kmul(rm1,ToReal(2)),kmul(R4p1323,rm3))),kmul(nn,kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2))))))))),kmadd(n3,kmul(ToReal(2),kmsub(nn,kmadd(rm2,kadd(Ro132,Ro231),kmadd(rm3,kadd(Ro133,Ro331),kmul(rm1,kmul(Ro131,ToReal(2))))),kmul(n1,kmadd(R4p1213,rm2,kmul(R4p1313,rm3))))),kmul(nn,kmadd(n1,kmul(ToReal(2),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2)))))),kmul(nn,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2)))))))))))),kmadd(im2,kmadd(rm1,kmul(Rojo12,kmul(nn,nn)),kmadd(rm1,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo23,kmul(nn,nn)),kmadd(rm3,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro132,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro231,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro233,ToReal(2)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro332,ToReal(2)))),kmadd(kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmul(kmul(n1,n1),ToReal(2)),kmadd(R4p1323,kmul(rm1,kmul(kmul(n3,n3),ToReal(2))),kmadd(R4p2323,kmul(rm2,kmul(kmul(n3,n3),ToReal(2))),kmadd(rm2,kmul(Rojo22,kmul(kmul(nn,nn),ToReal(2))),kmadd(n1,kmul(ToReal(2),kmadd(n2,kmsub(R4p1223,rm3,kmul(R4p1212,rm1)),kmsub(nn,kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2))))),kmul(n3,kmadd(R4p1213,rm1,kmadd(R4p1323,rm3,kmul(R4p1223,kmul(rm2,ToReal(2))))))))),kmadd(n2,kmul(ToReal(2),kmadd(n3,kmsub(R4p1223,rm1,kmul(R4p2323,rm3)),kmul(nn,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))))))))))))))))))),kmul(im3,kmadd(kmadd(R4p1213,rm2,kmul(R4p1313,rm3)),kmul(kmul(n1,n1),ToReal(2)),kmadd(kmul(n2,n2),kmadd(R4p1223,kmul(rm1,ToReal(-2)),kmul(R4p2323,kmul(rm3,ToReal(2)))),kmadd(n1,kmul(ToReal(2),knmsub(n3,kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmadd(n2,kmadd(R4p1223,rm2,kmsub(R4p1323,kmul(rm3,ToReal(2)),kmul(R4p1213,rm1))),kmul(nn,kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2))))))))),kmadd(n2,kmul(ToReal(2),kmsub(nn,kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2))))),kmul(n3,kmadd(R4p1323,rm1,kmul(R4p2323,rm2))))),kmul(nn,kmadd(n3,kmul(ToReal(2),kmadd(rm1,kadd(Ro133,Ro331),kmadd(rm2,kadd(Ro233,Ro332),kmul(rm3,kmul(Ro333,ToReal(2)))))),kmul(nn,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,ToReal(2))))))))))))))); + kmadd(kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kmadd(n2,R4p1213,kmul(n3,R4p1313)),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2)))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),ToReal(-1.)))),kmul(nn,kmadd(kmadd(im1,kmul(rm1,kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232)))),kmul(im3,kmul(rm3,kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))))))),ToReal(-2.),kmul(kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kadd(Ro112,Ro211),kmadd(n2,kadd(Ro122,Ro221),kmul(n3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kadd(Ro113,Ro311),kmadd(n2,kadd(Ro123,Ro321),kmul(n3,kadd(Ro133,Ro331)))),kmul(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(n1,kadd(Ro213,Ro312),kmadd(n2,kadd(Ro223,Ro322),kmul(n3,kadd(Ro233,Ro332))))))),ToReal(-1.)))))),ToReal(2.),kmadd(ToReal(-2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(n1,n1),kmadd(n2,kmul(n3,kmul(R4p2323,ToReal(-1.))),kmul(n1,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2.))))))),kmul(im3,kmul(rm3,kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2.))))))))))),kmul(kmul(nn,nn),kmul(ToReal(-1.),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,ToReal(2.)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi3rL = - kmadd(ltet2,kmul(R4p1212,kmul(rm2,kmul(n1,n1))),kmadd(ltet3,kmul(R4p1213,kmul(rm2,kmul(n1,n1))),kmadd(ltet2,kmul(R4p1213,kmul(rm3,kmul(n1,n1))),kmadd(ltet3,kmul(R4p1313,kmul(rm3,kmul(n1,n1))),kmadd(ltet1,kmul(R4p1212,kmul(rm1,kmul(n2,n2))),kmadd(ltet3,kmul(R4p2323,kmul(rm3,kmul(n2,n2))),kmadd(ltet1,kmul(R4p1313,kmul(rm1,kmul(n3,n3))),kmadd(ltet2,kmul(R4p1323,kmul(rm1,kmul(n3,n3))),kmadd(ltet1,kmul(R4p1323,kmul(rm2,kmul(n3,n3))),kmadd(ltet2,kmul(R4p2323,kmul(rm2,kmul(n3,n3))),kmadd(n1,kmul(nn,kmul(rm2,kmul(Ro221,kmadd(ltet2,ToReal(-2),n2)))),kmadd(n2,kmul(nn,kmul(rm2,kmul(Ro222,kmadd(ltet2,ToReal(-2),n2)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro223,kmadd(ltet2,ToReal(-2),n2)))),kmadd(n1,kmul(nn,kmul(rm3,kmul(Ro331,kmadd(ltet3,ToReal(-2),n3)))),kmadd(n2,kmul(nn,kmul(rm3,kmul(Ro332,kmadd(ltet3,ToReal(-2),n3)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro333,kmadd(ltet3,ToReal(-2),n3)))),kmadd(rm1,kmul(Rojo11,kmul(kmul(nn,nn),ksub(ltet1,n1))),kmadd(rm2,kmul(Rojo12,kmul(kmul(nn,nn),ksub(ltet1,n1))),kmadd(rm3,kmul(Rojo13,kmul(kmul(nn,nn),ksub(ltet1,n1))),kmadd(rm1,kmul(Rojo21,kmul(kmul(nn,nn),ksub(ltet2,n2))),kmadd(rm2,kmul(Rojo22,kmul(kmul(nn,nn),ksub(ltet2,n2))),kmadd(rm3,kmul(Rojo23,kmul(kmul(nn,nn),ksub(ltet2,n2))),kmadd(rm1,kmul(Rojo31,kmul(kmul(nn,nn),ksub(ltet3,n3))),kmadd(rm2,kmul(Rojo32,kmul(kmul(nn,nn),ksub(ltet3,n3))),kmadd(rm3,kmul(Rojo33,kmul(kmul(nn,nn),ksub(ltet3,n3))),kmadd(n2,kmul(nn,kmul(Ro212,kmsub(rm1,ksub(n2,ltet2),kmul(ltet1,rm2)))),kmadd(n2,kmul(nn,kmul(Ro232,kmsub(rm3,ksub(n2,ltet2),kmul(ltet3,rm2)))),kmadd(n3,kmul(nn,kmul(Ro313,kmsub(rm1,ksub(n3,ltet3),kmul(ltet1,rm3)))),kmadd(n3,kmul(nn,kmul(Ro323,kmsub(rm2,ksub(n3,ltet3),kmul(ltet2,rm3)))),kmadd(ltet3,kmsub(n2,kmul(rm2,kmsub(n1,R4p1223,kmul(n3,R4p2323))),kmul(R4p1223,kmul(rm1,kmul(n2,n2)))),kmadd(ltet1,kmsub(n3,kmsub(n2,kmul(R4p1223,rm2),kmul(n1,kmul(R4p1313,rm3))),kmul(R4p1223,kmul(rm3,kmul(n2,n2)))),kmadd(n3,kmsub(rm1,kmsub(ltet2,kmul(n2,R4p1223),kmul(ltet3,kmul(n1,R4p1313))),kmadd(ltet1,kmadd(n2,kmul(R4p1323,rm3),kmul(n1,kmul(R4p1213,rm2))),kmul(nn,kmul(Ro123,kmadd(ltet2,rm1,kmul(rm2,ksub(ltet1,n1))))))),kmadd(n1,kmsub(ltet2,kmul(n3,kmul(R4p1223,kmul(rm2,ToReal(-2)))),kmadd(ltet1,kmul(n2,kmul(R4p1213,rm3)),kmul(nn,kmul(Ro211,kmadd(ltet1,rm2,kmul(rm1,ksub(ltet2,n2))))))),kmadd(nn,kmsub(n1,kmul(rm1,kmul(Ro111,kmadd(ltet1,ToReal(-2),n1))),kmul(n3,kmul(Ro213,kmadd(ltet1,rm2,kmul(rm1,ksub(ltet2,n2)))))),kmadd(nn,kmsub(n2,kmul(rm1,kmul(Ro112,kmadd(ltet1,ToReal(-2),n1))),kmul(n3,kmul(Ro233,kmadd(ltet3,rm2,kmul(rm3,ksub(ltet2,n2)))))),kmadd(nn,kmsub(n3,kmul(rm1,kmul(Ro113,kmadd(ltet1,ToReal(-2),n1))),kmul(n1,kmul(Ro311,kmadd(ltet1,rm3,kmul(rm1,ksub(ltet3,n3)))))),kmadd(nn,kmsub(n1,kmul(Ro121,kmsub(rm2,ksub(n1,ltet1),kmul(ltet2,rm1))),kmul(n2,kmul(Ro312,kmadd(ltet1,rm3,kmul(rm1,ksub(ltet3,n3)))))),kmadd(n1,kmul(nn,kmsub(Ro131,kmsub(rm3,ksub(n1,ltet1),kmul(ltet3,rm1)),kmul(Ro321,kmadd(ltet2,rm3,kmul(rm2,ksub(ltet3,n3)))))),knmsub(n1,kmadd(ltet2,kmul(n3,kmul(R4p1323,rm3)),kmadd(n2,kmadd(rm1,kmadd(ltet3,R4p1213,kmul(ltet2,R4p1212)),kmul(ltet1,kmul(R4p1212,rm2))),kmul(nn,kmul(Ro231,kmadd(ltet3,rm2,kmul(rm3,ksub(ltet2,n2))))))),kmadd(n2,kmsub(ltet2,kmul(n1,kmul(R4p1223,rm3)),kmul(nn,kmadd(Ro322,kmadd(ltet2,rm3,kmul(rm2,ksub(ltet3,n3))),kmul(Ro122,kmadd(ltet2,rm1,kmul(rm2,ksub(ltet1,n1))))))),kmadd(n2,kmsub(ltet3,kmul(n1,kmul(R4p1323,kmul(rm3,ToReal(2)))),kmul(nn,kmul(Ro132,kmadd(ltet3,rm1,kmul(rm3,ksub(ltet1,n1)))))),kmul(n3,knmsub(ltet2,kmul(n2,kmul(R4p2323,rm3)),knmsub(ltet3,kmul(R4p1323,kmadd(n1,rm2,kmul(n2,rm1))),kmsub(R4p1213,kmul(rm1,kmsub(ltet1,kmul(n2,ToReal(2)),kmul(ltet2,n1))),kmul(nn,kmul(Ro133,kmadd(ltet3,rm1,kmul(rm3,ksub(ltet1,n1)))))))))))))))))))))))))))))))))))))))))))))))))); + kmadd(kmul(nn,nn),kmul(ToReal(-1.),kmadd(kmadd(rm1,Rojo11,kmadd(rm2,Rojo12,kmul(rm3,Rojo13))),kmadd(ltet1,ToReal(-1.),n1),kmadd(kmadd(rm1,Rojo21,kmadd(rm2,Rojo22,kmul(rm3,Rojo23))),kmadd(ltet2,ToReal(-1.),n2),kmul(kmadd(rm1,Rojo31,kmadd(rm2,Rojo32,kmul(rm3,Rojo33))),kmadd(ltet3,ToReal(-1.),n3))))),kmadd(n3,kmadd(ltet1,kmadd(kmadd(n2,R4p1213,kmul(n3,R4p1313)),rm1,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),rm2)),kmadd(ltet2,kmadd(n3,kmadd(R4p1323,rm1,kmul(R4p2323,rm2)),kmul(n1,kmul(kmadd(R4p1213,rm1,kmul(R4p1223,rm2)),ToReal(-1.)))),kmadd(nn,kmadd(rm1,kmul(Ro113,kmadd(ltet1,ToReal(-2.),n1)),kmadd(rm2,kmul(Ro223,kmadd(ltet2,ToReal(-2.),n2)),kmadd(rm3,kmul(Ro333,kmadd(ltet3,ToReal(-2.),n3)),kmadd(Ro313,kmadd(ltet1,kmul(rm3,ToReal(-1.)),kmul(rm1,kmadd(ltet3,ToReal(-1.),n3))),kmadd(Ro323,kmadd(ltet2,kmul(rm3,ToReal(-1.)),kmul(rm2,kmadd(ltet3,ToReal(-1.),n3))),kmul(Ro133,kmul(ToReal(-1.),kmadd(ltet3,rm1,kmul(rm3,kmadd(n1,ToReal(-1.),ltet1)))))))))),kmul(ToReal(-1.),kmadd(ltet3,kmadd(kmadd(n1,R4p1313,kmul(n2,R4p1323)),rm1,kmul(kmadd(n1,R4p1323,kmul(n2,R4p2323)),rm2)),kmul(nn,kmadd(Ro123,kmadd(ltet2,rm1,kmul(rm2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro213,kmadd(ltet1,rm2,kmul(rm1,kmadd(n2,ToReal(-1.),ltet2))),kmul(Ro233,kmadd(ltet3,rm2,kmul(rm3,kmadd(n2,ToReal(-1.),ltet2)))))))))))),kmadd(n1,kmadd(ltet3,kmadd(kmadd(n1,R4p1213,kmul(n2,R4p1223)),rm2,kmul(kmadd(n1,R4p1313,kmul(n2,R4p1323)),rm3)),kmadd(ltet2,kmadd(n1,kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmul(n3,kmul(kmadd(R4p1223,rm2,kmul(R4p1323,rm3)),ToReal(-1.)))),kmadd(nn,kmadd(rm1,kmul(Ro111,kmadd(ltet1,ToReal(-2.),n1)),kmadd(rm2,kmul(Ro221,kmadd(ltet2,ToReal(-2.),n2)),kmadd(rm3,kmul(Ro331,kmadd(ltet3,ToReal(-2.),n3)),kmadd(Ro121,kmadd(ltet2,kmul(rm1,ToReal(-1.)),kmul(rm2,kmadd(ltet1,ToReal(-1.),n1))),kmadd(Ro131,kmadd(ltet3,kmul(rm1,ToReal(-1.)),kmul(rm3,kmadd(ltet1,ToReal(-1.),n1))),kmul(Ro311,kmul(ToReal(-1.),kmadd(ltet1,rm3,kmul(rm1,kmadd(n3,ToReal(-1.),ltet3)))))))))),kmul(ToReal(-1.),kmadd(ltet1,kmadd(kmadd(n2,R4p1212,kmul(n3,R4p1213)),rm2,kmul(kmadd(n2,R4p1213,kmul(n3,R4p1313)),rm3)),kmul(nn,kmadd(Ro211,kmadd(ltet1,rm2,kmul(rm1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro231,kmadd(ltet3,rm2,kmul(rm3,kmadd(n2,ToReal(-1.),ltet2))),kmul(Ro321,kmadd(ltet2,rm3,kmul(rm2,kmadd(n3,ToReal(-1.),ltet3)))))))))))),kmul(n2,kmadd(ltet3,kmadd(kmadd(n1,R4p1323,kmul(n2,R4p2323)),rm3,kmul(kmadd(n1,R4p1213,kmul(n2,R4p1223)),kmul(rm1,ToReal(-1.)))),kmadd(ltet1,kmadd(kmadd(n2,R4p1212,kmul(n3,R4p1213)),rm1,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),kmul(rm3,ToReal(-1.)))),kmadd(ltet2,kmadd(rm1,kmadd(n3,R4p1223,kmul(n1,kmul(R4p1212,ToReal(-1.)))),kmul(rm3,kmadd(n1,R4p1223,kmul(n3,kmul(R4p2323,ToReal(-1.)))))),kmul(nn,kmadd(rm1,kmul(Ro112,kmadd(ltet1,ToReal(-2.),n1)),kmadd(rm2,kmul(Ro222,kmadd(ltet2,ToReal(-2.),n2)),kmadd(rm3,kmul(Ro332,kmadd(ltet3,ToReal(-2.),n3)),kmadd(Ro212,kmadd(ltet1,kmul(rm2,ToReal(-1.)),kmul(rm1,kmadd(ltet2,ToReal(-1.),n2))),kmadd(Ro232,kmadd(ltet3,kmul(rm2,ToReal(-1.)),kmul(rm3,kmadd(ltet2,ToReal(-1.),n2))),kmul(ToReal(-1.),kmadd(Ro122,kmadd(ltet2,rm1,kmul(rm2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro132,kmadd(ltet3,rm1,kmul(rm3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro312,kmadd(ltet1,rm3,kmul(rm1,kmadd(n3,ToReal(-1.),ltet3))),kmul(Ro322,kmadd(ltet2,rm3,kmul(rm2,kmadd(n3,ToReal(-1.),ltet3))))))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi3iL = - kmadd(im2,kmadd(n1,kmul(Rojo12,kmul(nn,nn)),kmadd(n2,kmul(Rojo22,kmul(nn,nn)),kmadd(n3,kmul(Rojo32,kmul(nn,nn)),knmsub(ltet2,kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmadd(n1,kmul(kmadd(n3,R4p1223,kmul(nn,Ro221)),ToReal(-2)),kmadd(n3,kmul(nn,kmul(Ro223,ToReal(-2))),kmul(nn,kmadd(nn,Rojo22,kmul(n2,kmul(Ro222,ToReal(-2))))))))),kmadd(nn,kmsub(ltet3,kmul(n3,Ro233),kmul(Ro121,kmul(n1,n1))),kmadd(nn,kmsub(ltet3,kmul(n1,Ro321),kmul(Ro222,kmul(n2,n2))),kmadd(nn,kmsub(ltet3,kmul(n2,Ro322),kmul(Ro323,kmul(n3,n3))),kmadd(ltet3,kmsub(n3,kmul(nn,Ro323),kmul(Rojo32,kmul(nn,nn))),kmadd(ltet3,kmadd(n2,kmul(n3,R4p2323),kmadd(n1,kmul(nn,Ro231),kmsub(n1,kmsub(n3,R4p1323,kmul(n2,R4p1223)),kmul(R4p1213,kmul(n1,n1))))),kmadd(ltet1,kmadd(n3,kmul(nn,Ro123),kmadd(n1,kmadd(n2,R4p1212,kmadd(n3,R4p1213,kmul(nn,kadd(Ro121,Ro211)))),kmadd(n2,kmul(nn,Ro212),kmadd(n3,kmul(nn,Ro213),knmsub(R4p1323,kmul(n3,n3),kmsub(n2,kmsub(nn,Ro122,kmul(n3,R4p1223)),kmul(Rojo12,kmul(nn,nn)))))))),kmsub(n2,kmul(nn,kmsub(ltet3,Ro232,kmul(n3,kadd(Ro322,Ro223)))),kmul(n1,kmul(nn,kmadd(n3,kadd(Ro321,Ro123),kmul(n2,kadd(Ro221,Ro122)))))))))))))))),kmadd(im1,kmadd(ltet3,kmul(R4p1223,kmul(n2,n2)),kmadd(n1,kmul(Rojo11,kmul(nn,nn)),kmadd(n2,kmul(Rojo21,kmul(nn,nn)),kmadd(n3,kmul(Rojo31,kmul(nn,nn)),kmadd(n3,kmul(nn,kmsub(ltet3,Ro133,kmul(n1,Ro311))),kmadd(nn,kmsub(ltet3,kmul(n1,Ro311),kmul(Ro111,kmul(n1,n1))),kmadd(nn,kmsub(ltet3,kmul(n2,Ro312),kmul(Ro212,kmul(n2,n2))),kmadd(nn,kmul(Ro313,kmsub(ltet3,n3,kmul(n3,n3))),kmadd(ltet3,kmadd(n1,kmul(n2,R4p1213),kmadd(n1,kmul(n3,R4p1313),kmadd(n2,kmul(n3,R4p1323),kmsub(n1,kmul(nn,Ro131),kmul(Rojo31,kmul(nn,nn)))))),kmadd(ltet2,kmadd(n3,kmul(nn,Ro123),kmadd(n1,kmadd(n2,R4p1212,kmadd(n3,R4p1213,kmul(nn,kadd(Ro121,Ro211)))),kmadd(n2,kmul(nn,Ro212),kmadd(n3,kmul(nn,Ro213),knmsub(R4p1323,kmul(n3,n3),kmsub(n2,kmsub(nn,Ro122,kmul(n3,R4p1223)),kmul(Rojo21,kmul(nn,nn)))))))),knmsub(n1,kmul(nn,kmadd(n2,kadd(Ro211,Ro112),kmul(n3,Ro113))),kmsub(n2,kmul(nn,kmsub(ltet3,Ro132,kmul(n3,kadd(Ro312,Ro213)))),kmul(ltet1,kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmadd(Rojo11,kmul(nn,nn),kmadd(n1,kmul(nn,kmul(Ro111,ToReal(-2))),kmadd(n2,kmul(nn,kmul(Ro112,ToReal(-2))),kmadd(n3,kmul(nn,kmul(Ro113,ToReal(-2))),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2))))))))))))))))))))))),kmul(im3,kmadd(n1,kmul(Rojo13,kmul(nn,nn)),kmadd(n2,kmul(Rojo23,kmul(nn,nn)),kmadd(n3,kmul(Rojo33,kmul(nn,nn)),knmsub(n3,kmul(nn,kmadd(n2,kadd(Ro332,Ro233),kmul(n1,Ro331))),kmadd(ltet1,kmadd(n1,kmadd(n2,R4p1213,kmadd(n3,R4p1313,kmul(nn,kadd(Ro131,Ro311)))),kmadd(n2,kmadd(n3,R4p1323,kmul(nn,kadd(Ro132,Ro312))),kmadd(R4p1223,kmul(n2,n2),kmul(nn,kmsub(n3,kadd(Ro133,Ro313),kmul(nn,Rojo13)))))),kmadd(ltet2,kmadd(n2,kmadd(n3,R4p2323,kmul(nn,kadd(Ro232,Ro322))),knmsub(R4p1213,kmul(n1,n1),kmadd(n1,kmadd(n3,R4p1323,kmsub(nn,kadd(Ro231,Ro321),kmul(n2,R4p1223))),kmul(nn,kmsub(n3,kadd(Ro233,Ro323),kmul(nn,Rojo23)))))),kmadd(ltet3,knmsub(R4p1313,kmul(n1,n1),knmsub(R4p2323,kmul(n2,n2),kmsub(n2,kmul(nn,kmul(Ro332,ToReal(2))),kmul(Rojo33,kmul(nn,nn))))),kmadd(nn,knmsub(Ro131,kmul(n1,n1),kmsub(Ro333,kmsub(ltet3,kmul(n3,ToReal(2)),kmul(n3,n3)),kmul(Ro232,kmul(n2,n2)))),kmul(n1,kmadd(n2,kmsub(ltet3,kmul(R4p1323,ToReal(-2)),kmul(nn,kadd(Ro231,Ro132))),kmul(nn,kmsub(ltet3,kmul(Ro331,ToReal(2)),kmul(n3,Ro133)))))))))))))))); + kmadd(kmul(nn,nn),kmadd(kmadd(im1,Rojo11,kmadd(im2,Rojo12,kmul(im3,Rojo13))),kmadd(ltet1,ToReal(-1.),n1),kmadd(kmadd(im1,Rojo21,kmadd(im2,Rojo22,kmul(im3,Rojo23))),kmadd(ltet2,ToReal(-1.),n2),kmul(kmadd(im1,Rojo31,kmadd(im2,Rojo32,kmul(im3,Rojo33))),kmadd(ltet3,ToReal(-1.),n3)))),kmadd(ToReal(-1.),kmadd(n3,kmadd(im1,kmadd(R4p1213,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1313,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p1323,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.))))))),kmul(im2,kmadd(R4p1223,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1323,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p2323,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.))))))))),kmadd(n1,kmadd(im2,kmadd(R4p1212,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1213,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p1223,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.))))))),kmul(im3,kmadd(R4p1213,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1313,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p1323,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.))))))))),kmul(n2,kmadd(im1,kmadd(R4p1212,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1213,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p1223,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.))))))),kmul(im3,kmadd(R4p1223,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1323,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p2323,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.)))))))))))),kmul(nn,kmadd(kmadd(n1,Ro121,kmadd(n2,Ro122,kmul(n3,Ro123))),kmadd(im1,ltet2,kmul(im2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(kmadd(n1,Ro131,kmadd(n2,Ro132,kmul(n3,Ro133))),kmadd(im1,ltet3,kmul(im3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(kmadd(n1,Ro211,kmadd(n2,Ro212,kmul(n3,Ro213))),kmadd(im2,ltet1,kmul(im1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(kmadd(n1,Ro231,kmadd(n2,Ro232,kmul(n3,Ro233))),kmadd(im2,ltet3,kmul(im3,kmadd(n2,ToReal(-1.),ltet2))),kmadd(kmadd(n1,Ro311,kmadd(n2,Ro312,kmul(n3,Ro313))),kmadd(im3,ltet1,kmul(im1,kmadd(n3,ToReal(-1.),ltet3))),kmadd(kmadd(n1,Ro321,kmadd(n2,Ro322,kmul(n3,Ro323))),kmadd(im3,ltet2,kmul(im2,kmadd(n3,ToReal(-1.),ltet3))),kmadd(im1,kmul(kmadd(n1,Ro111,kmadd(n2,Ro112,kmul(n3,Ro113))),kmadd(n1,ToReal(-1.),kmul(ltet1,ToReal(2.)))),kmadd(im2,kmul(kmadd(n1,Ro221,kmadd(n2,Ro222,kmul(n3,Ro223))),kmadd(n2,ToReal(-1.),kmul(ltet2,ToReal(2.)))),kmul(im3,kmul(kmadd(n1,Ro331,kmadd(n2,Ro332,kmul(n3,Ro333))),kmadd(n3,ToReal(-1.),kmul(ltet3,ToReal(2.)))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi2rL = - kmadd(ltet1,kmul(n3,kmul(R4p1223,kmadd(im2,im2,kmul(rm2,rm2)))),kmadd(nn,kmul(kmadd(im1,im2,kmul(rm1,rm2)),kmul(Ro211,ksub(n1,ltet1))),kmadd(nn,kmul(kmadd(im1,im3,kmul(rm1,rm3)),kmul(Ro311,ksub(n1,ltet1))),kmadd(nn,kmul(kmadd(im2,im3,kmul(rm2,rm3)),kmul(Ro322,ksub(n2,ltet2))),kmadd(nn,kmul(Ro222,kmul(kmadd(im2,im2,kmul(rm2,rm2)),ksub(n2,ltet2))),kmadd(nn,kmul(kmadd(im1,im3,kmul(rm1,rm3)),kmul(Ro133,ksub(n3,ltet3))),kmadd(nn,kmul(kmadd(im2,im3,kmul(rm2,rm3)),kmul(Ro233,ksub(n3,ltet3))),kmadd(nn,kmul(Ro112,kmsub(n2,kmadd(im1,im1,kmul(rm1,rm1)),kmul(ltet1,kmadd(im1,im2,kmul(rm1,rm2))))),kmadd(nn,kmul(Ro221,kmsub(n1,kmadd(im2,im2,kmul(rm2,rm2)),kmul(ltet2,kmadd(im1,im2,kmul(rm1,rm2))))),kmadd(nn,kmul(Ro113,kmsub(n3,kmadd(im1,im1,kmul(rm1,rm1)),kmul(ltet1,kmadd(im1,im3,kmul(rm1,rm3))))),kmadd(nn,kmul(Ro331,kmsub(n1,kmadd(im3,im3,kmul(rm3,rm3)),kmul(ltet3,kmadd(im1,im3,kmul(rm1,rm3))))),kmadd(nn,kmul(Ro223,kmsub(n3,kmadd(im2,im2,kmul(rm2,rm2)),kmul(ltet2,kmadd(im2,im3,kmul(rm2,rm3))))),kmadd(nn,kmul(Ro332,kmsub(n2,kmadd(im3,im3,kmul(rm3,rm3)),kmul(ltet3,kmadd(im2,im3,kmul(rm2,rm3))))),kmadd(kmadd(im1,im2,kmul(rm1,rm2)),kmsub(ltet3,kmul(n3,kmul(R4p1323,ToReal(-2))),kmul(Rojo12,kmul(nn,nn))),kmadd(kmadd(im1,im3,kmul(rm1,rm3)),kmsub(ltet2,kmul(n1,R4p1213),kmul(Rojo13,kmul(nn,nn))),kmadd(kmadd(im2,im2,kmul(rm2,rm2)),kmsub(ltet3,kmul(n1,R4p1223),kmul(Rojo22,kmul(nn,nn))),kmadd(kmadd(im2,im3,kmul(rm2,rm3)),kmsub(ltet1,kmul(n1,kmul(R4p1213,ToReal(-2))),kmul(Rojo23,kmul(nn,nn))),kmadd(kmadd(im2,im3,kmul(rm2,rm3)),kmsub(ltet3,kmul(n1,R4p1323),kmul(Rojo32,kmul(nn,nn))),kmadd(nn,kmul(Ro121,kmsub(n1,kmadd(im1,im2,kmul(rm1,rm2)),kmul(ltet2,kmadd(im1,im1,kmul(rm1,rm1))))),kmadd(nn,kmul(Ro131,kmsub(n1,kmadd(im1,im3,kmul(rm1,rm3)),kmul(ltet3,kmadd(im1,im1,kmul(rm1,rm1))))),kmadd(ltet3,kmul(R4p1313,kmsub(n1,kmadd(im1,im3,kmul(rm1,rm3)),kmul(n3,kmadd(im1,im1,kmul(rm1,rm1))))),kmadd(n3,kmsub(ltet1,kmul(R4p1313,kmadd(im1,im3,kmul(rm1,rm3))),kmul(ltet2,kmul(R4p1213,kmadd(im1,im1,kmul(rm1,rm1))))),kmadd(nn,kmul(Ro212,kmsub(n2,kmadd(im1,im2,kmul(rm1,rm2)),kmul(ltet1,kmadd(im2,im2,kmul(rm2,rm2))))),kmadd(nn,kmul(Ro232,kmsub(n2,kmadd(im2,im3,kmul(rm2,rm3)),kmul(ltet3,kmadd(im2,im2,kmul(rm2,rm2))))),kmadd(ltet1,kmsub(n3,kmul(R4p1323,kmadd(im2,im3,kmul(rm2,rm3))),kmul(n1,kmul(R4p1212,kmadd(im2,im2,kmul(rm2,rm2))))),kmadd(ltet3,kmsub(n2,kmul(R4p1323,kmadd(im1,im3,kmul(rm1,rm3))),kmul(n3,kmul(R4p2323,kmadd(im2,im2,kmul(rm2,rm2))))),kmadd(nn,kmul(Ro313,kmsub(n3,kmadd(im1,im3,kmul(rm1,rm3)),kmul(ltet1,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(nn,kmul(Ro323,kmsub(n3,kmadd(im2,im3,kmul(rm2,rm3)),kmul(ltet2,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(ltet2,kmul(R4p1323,kmsub(n3,kmadd(im1,im3,kmul(rm1,rm3)),kmul(n1,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(ltet2,kmul(R4p2323,kmsub(n3,kmadd(im2,im3,kmul(rm2,rm3)),kmul(n2,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(n2,kmsub(ltet3,kmul(R4p2323,kmadd(im2,im3,kmul(rm2,rm3))),kmul(ltet1,kmul(R4p1323,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(kmadd(im1,im1,kmul(rm1,rm1)),kmsub(nn,kmul(Ro111,ksub(n1,ltet1)),kmul(Rojo11,kmul(nn,nn))),kmadd(kmadd(im1,im2,kmul(rm1,rm2)),kmsub(nn,kmul(Ro122,ksub(n2,ltet2)),kmul(Rojo21,kmul(nn,nn))),kmadd(kmadd(im3,im3,kmul(rm3,rm3)),kmsub(nn,kmul(Ro333,ksub(n3,ltet3)),kmul(Rojo33,kmul(nn,nn))),knmsub(R4p1223,kmadd(kmadd(im1,im2,kmul(rm1,rm2)),kmadd(ltet2,n3,kmul(ltet3,n2)),kmul(ltet2,kmul(n1,kmadd(im2,im3,kmul(rm2,rm3))))),kmadd(nn,kmul(Ro321,kmadd(im3,kmsub(im2,n1,kmul(im1,ltet2)),kmul(rm3,kmsub(n1,rm2,kmul(ltet2,rm1))))),kmadd(nn,kmul(Ro231,kmadd(im2,kmsub(im3,n1,kmul(im1,ltet3)),kmul(rm2,kmsub(n1,rm3,kmul(ltet3,rm1))))),kmadd(nn,kmul(Ro312,kmadd(im3,kmsub(im1,n2,kmul(im2,ltet1)),kmul(rm3,kmsub(n2,rm1,kmul(ltet1,rm2))))),kmadd(nn,kmul(Ro132,kmadd(im1,kmsub(im3,n2,kmul(im2,ltet3)),kmul(rm1,kmsub(n2,rm3,kmul(ltet3,rm2))))),kmadd(nn,kmul(Ro213,kmadd(im2,kmsub(im1,n3,kmul(im3,ltet1)),kmul(rm2,kmsub(n3,rm1,kmul(ltet1,rm3))))),kmadd(nn,kmul(Ro123,kmadd(im1,kmsub(im2,n3,kmul(im3,ltet2)),kmul(rm1,kmsub(n3,rm2,kmul(ltet2,rm3))))),kmadd(ltet1,kmsub(n2,kmsub(R4p1213,kmadd(im1,im3,kmul(rm1,rm3)),kmul(R4p1223,kmadd(im2,im3,kmul(rm2,rm3)))),kmul(n1,kmul(R4p1313,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(R4p1212,kmadd(ltet1,kmul(n2,kmadd(im1,im2,kmul(rm1,rm2))),kmul(ltet2,kmsub(n1,kmadd(im1,im2,kmul(rm1,rm2)),kmul(n2,kmadd(im1,im1,kmul(rm1,rm1)))))),kmadd(R4p1213,kmadd(ltet1,kmul(n3,kmadd(im1,im2,kmul(rm1,rm2))),kmul(ltet3,kmsub(n1,kmadd(im1,im2,kmul(rm1,rm2)),kmul(n2,kmadd(im1,im1,kmul(rm1,rm1)))))),kmul(kmadd(im1,im3,kmul(rm1,rm3)),kmsub(ltet2,kmul(n2,kmul(R4p1223,ToReal(2))),kmul(Rojo31,kmul(nn,nn)))))))))))))))))))))))))))))))))))))))))))))))); + kmadd(kmadd(im1,im2,kmul(rm1,rm2)),kmadd(n1,kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),kmadd(ltet1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(kmadd(n3,kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),kmul(ltet3,kmadd(n2,R4p1223,kmul(n3,R4p1323)))),ToReal(-1.)))),kmadd(ToReal(-1.),kmadd(kmadd(n2,kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),kmul(n3,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)))),kmadd(im1,im1,kmul(rm1,rm1)),kmadd(kmadd(n1,kmadd(ltet1,R4p1313,kmul(ltet2,R4p1323)),kmul(n2,kmadd(ltet1,R4p1323,kmul(ltet2,R4p2323)))),kmadd(im3,im3,kmul(rm3,rm3)),kmadd(kmul(nn,nn),kmadd(im1,kmadd(im2,kadd(Rojo12,Rojo21),kmul(im3,kadd(Rojo13,Rojo31))),kmadd(rm1,kmadd(rm2,kadd(Rojo12,Rojo21),kmul(rm3,kadd(Rojo13,Rojo31))),kmadd(kmadd(im2,im3,kmul(rm2,rm3)),kadd(Rojo23,Rojo32),kmadd(Rojo11,kmadd(im1,im1,kmul(rm1,rm1)),kmadd(Rojo22,kmadd(im2,im2,kmul(rm2,rm2)),kmul(Rojo33,kmadd(im3,im3,kmul(rm3,rm3)))))))),kmadd(kmadd(im2,im2,kmul(rm2,rm2)),kmadd(n3,kmadd(ltet3,R4p2323,kmul(ltet1,kmul(R4p1223,ToReal(-1.)))),kmul(n1,kmadd(ltet1,R4p1212,kmul(ltet3,kmul(R4p1223,ToReal(-1.)))))),kmul(kmadd(im2,im3,kmul(rm2,rm3)),kmadd(n1,kmadd(ltet1,R4p1213,kmul(ltet2,R4p1223)),kmadd(kmadd(ltet2,kmul(n3,R4p2323),kmul(ltet3,kmadd(n1,R4p1323,kmul(n2,R4p2323)))),ToReal(-1.),kmul(ltet1,kmadd(n1,R4p1213,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))))))))),kmadd(nn,kmadd(kmadd(kmadd(im1,im3,kmul(rm1,rm3)),Ro311,kmul(Ro111,kmadd(im1,im1,kmul(rm1,rm1)))),kmadd(ltet1,ToReal(-1.),n1),kmadd(kmadd(kmadd(im2,im3,kmul(rm2,rm3)),Ro322,kmul(Ro222,kmadd(im2,im2,kmul(rm2,rm2)))),kmadd(ltet2,ToReal(-1.),n2),kmadd(kmadd(kmadd(im1,im3,kmul(rm1,rm3)),Ro133,kmadd(kmadd(im2,im3,kmul(rm2,rm3)),Ro233,kmul(Ro333,kmadd(im3,im3,kmul(rm3,rm3))))),kmadd(ltet3,ToReal(-1.),n3),kmadd(Ro112,kmadd(n2,kmadd(im1,im1,kmul(rm1,rm1)),kmul(ltet1,kmul(kmadd(im1,im2,kmul(rm1,rm2)),ToReal(-1.)))),kmadd(Ro221,kmadd(n1,kmadd(im2,im2,kmul(rm2,rm2)),kmul(ltet2,kmul(kmadd(im1,im2,kmul(rm1,rm2)),ToReal(-1.)))),kmadd(Ro113,kmadd(n3,kmadd(im1,im1,kmul(rm1,rm1)),kmul(ltet1,kmul(kmadd(im1,im3,kmul(rm1,rm3)),ToReal(-1.)))),kmadd(Ro331,kmadd(n1,kmadd(im3,im3,kmul(rm3,rm3)),kmul(ltet3,kmul(kmadd(im1,im3,kmul(rm1,rm3)),ToReal(-1.)))),kmadd(Ro223,kmadd(n3,kmadd(im2,im2,kmul(rm2,rm2)),kmul(ltet2,kmul(kmadd(im2,im3,kmul(rm2,rm3)),ToReal(-1.)))),kmadd(Ro332,kmadd(n2,kmadd(im3,im3,kmul(rm3,rm3)),kmul(ltet3,kmul(kmadd(im2,im3,kmul(rm2,rm3)),ToReal(-1.)))),kmadd(Ro121,kmadd(n1,kmadd(im1,im2,kmul(rm1,rm2)),kmul(ltet2,kmul(kmadd(im1,im1,kmul(rm1,rm1)),ToReal(-1.)))),kmadd(Ro131,kmadd(n1,kmadd(im1,im3,kmul(rm1,rm3)),kmul(ltet3,kmul(kmadd(im1,im1,kmul(rm1,rm1)),ToReal(-1.)))),kmadd(Ro212,kmadd(n2,kmadd(im1,im2,kmul(rm1,rm2)),kmul(ltet1,kmul(kmadd(im2,im2,kmul(rm2,rm2)),ToReal(-1.)))),kmadd(Ro232,kmadd(n2,kmadd(im2,im3,kmul(rm2,rm3)),kmul(ltet3,kmul(kmadd(im2,im2,kmul(rm2,rm2)),ToReal(-1.)))),kmadd(Ro313,kmadd(n3,kmadd(im1,im3,kmul(rm1,rm3)),kmul(ltet1,kmul(kmadd(im3,im3,kmul(rm3,rm3)),ToReal(-1.)))),kmadd(Ro323,kmadd(n3,kmadd(im2,im3,kmul(rm2,rm3)),kmul(ltet2,kmul(kmadd(im3,im3,kmul(rm3,rm3)),ToReal(-1.)))),kmadd(kmadd(im1,im2,kmul(rm1,rm2)),kmadd(Ro211,kmadd(ltet1,ToReal(-1.),n1),kmul(Ro122,kmadd(ltet2,ToReal(-1.),n2))),kmadd(Ro321,kmadd(im3,kmadd(im2,n1,kmul(im1,kmul(ltet2,ToReal(-1.)))),kmul(rm3,kmadd(n1,rm2,kmul(ltet2,kmul(rm1,ToReal(-1.)))))),kmadd(Ro231,kmadd(im2,kmadd(im3,n1,kmul(im1,kmul(ltet3,ToReal(-1.)))),kmul(rm2,kmadd(n1,rm3,kmul(ltet3,kmul(rm1,ToReal(-1.)))))),kmadd(Ro312,kmadd(im3,kmadd(im1,n2,kmul(im2,kmul(ltet1,ToReal(-1.)))),kmul(rm3,kmadd(n2,rm1,kmul(ltet1,kmul(rm2,ToReal(-1.)))))),kmadd(Ro132,kmadd(im1,kmadd(im3,n2,kmul(im2,kmul(ltet3,ToReal(-1.)))),kmul(rm1,kmadd(n2,rm3,kmul(ltet3,kmul(rm2,ToReal(-1.)))))),kmadd(Ro213,kmadd(im2,kmadd(im1,n3,kmul(im3,kmul(ltet1,ToReal(-1.)))),kmul(rm2,kmadd(n3,rm1,kmul(ltet1,kmul(rm3,ToReal(-1.)))))),kmul(Ro123,kmadd(im1,kmadd(im2,n3,kmul(im3,kmul(ltet2,ToReal(-1.)))),kmul(rm1,kmadd(n3,rm2,kmul(ltet2,kmul(rm3,ToReal(-1.)))))))))))))))))))))))))))),kmul(kmadd(im1,im3,kmul(rm1,rm3)),kmadd(kmadd(ltet2,n1,kmul(ltet1,n2)),R4p1213,kmadd(kmadd(ltet3,n1,kmul(ltet1,n3)),R4p1313,kmadd(ltet2,kmul(n3,R4p1323),kmul(n2,kmadd(ltet3,R4p1323,kmul(ltet2,kmul(R4p1223,ToReal(2)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi2iL = - kmadd(im3,kmadd(n2,kmul(nn,kmul(rm1,Ro312)),kmadd(n3,kmul(nn,kmul(rm1,Ro313)),kmadd(n2,kmul(nn,kmul(rm2,Ro322)),kmadd(n3,kmul(nn,kmul(rm2,Ro323)),kmadd(ltet3,kmul(nn,kmul(rm1,Ro331)),kmadd(ltet3,kmul(nn,kmul(rm2,Ro332)),kmadd(rm1,kmul(Rojo13,kmul(nn,nn)),kmadd(rm2,kmul(Rojo23,kmul(nn,nn)),kmadd(rm1,kmsub(n1,kmul(nn,Ro311),kmul(Rojo31,kmul(nn,nn))),kmadd(rm2,kmsub(n1,kmul(nn,Ro321),kmul(Rojo32,kmul(nn,nn))),kmadd(ltet1,kmadd(n3,kmul(R4p1313,rm1),kmadd(n3,kmul(R4p1323,rm2),kmadd(nn,kmul(rm1,Ro311),kmadd(nn,kmul(rm2,Ro312),kmadd(rm1,kmsub(n2,R4p1213,kmul(nn,Ro113)),kmul(rm2,kmsub(n2,R4p1223,kmul(nn,Ro213)))))))),kmsub(ltet2,kmadd(nn,kmul(rm1,Ro321),kmadd(nn,kmul(rm2,Ro322),knmsub(n1,kmadd(R4p1213,rm1,kmul(R4p1223,rm2)),kmadd(rm1,kmsub(n3,R4p1323,kmul(nn,Ro123)),kmul(rm2,kmsub(n3,R4p2323,kmul(nn,Ro223))))))),kmadd(ltet3,kmadd(rm2,kmadd(n1,R4p1323,kmadd(nn,Ro233,kmul(n2,R4p2323))),kmul(rm1,kmadd(n1,R4p1313,kmadd(nn,Ro133,kmul(n2,R4p1323))))),kmul(nn,kmadd(rm2,kmadd(n1,Ro231,kmadd(n3,Ro233,kmul(n2,Ro232))),kmul(rm1,kmadd(n1,Ro131,kmadd(n3,Ro133,kmul(n2,Ro132))))))))))))))))))),kmadd(im2,kmadd(n3,kmul(nn,kmul(rm1,Ro213)),kmadd(ltet3,kmul(nn,kmul(rm1,Ro231)),kmadd(n2,kmul(nn,kmul(rm3,Ro232)),kmadd(ltet3,kmul(nn,kmul(rm3,Ro233)),kmadd(n3,kmul(nn,kmul(rm3,Ro233)),kmadd(rm1,kmul(Rojo12,kmul(nn,nn)),kmadd(rm3,kmul(Rojo32,kmul(nn,nn)),kmadd(rm1,kmsub(n2,kmul(nn,Ro212),kmul(Rojo21,kmul(nn,nn))),kmadd(rm3,kmsub(n1,kmul(nn,Ro231),kmul(Rojo23,kmul(nn,nn))),kmadd(nn,kmsub(rm1,kmsub(n1,Ro211,kmul(n3,Ro123)),kmul(n3,kmul(rm3,Ro323))),knmsub(ltet3,kmadd(rm1,kmadd(n1,R4p1213,kmadd(nn,Ro132,kmul(n2,R4p1223))),kmul(nn,kmul(rm3,Ro332))),kmadd(n1,kmsub(rm3,kmsub(ltet3,R4p1323,kmul(nn,Ro321)),kmul(nn,kmul(rm1,Ro121))),kmadd(n2,kmsub(rm3,kmsub(ltet3,R4p2323,kmul(nn,Ro322)),kmul(nn,kmul(rm1,Ro122))),kmadd(ltet2,kmadd(nn,kmul(rm1,Ro221),kmadd(nn,kmul(rm3,Ro223),kmadd(rm1,kmsub(n3,R4p1223,kmadd(nn,Ro122,kmul(n1,R4p1212))),kmul(rm3,kmsub(n1,R4p1223,kmadd(nn,Ro322,kmul(n3,R4p2323))))))),kmul(ltet1,kmadd(nn,kmul(rm3,Ro213),kmadd(n2,kmsub(R4p1212,rm1,kmul(R4p1223,rm3)),kmadd(n3,kmsub(R4p1213,rm1,kmul(R4p1323,rm3)),kmul(nn,kmsub(rm1,ksub(Ro211,Ro112),kmul(rm3,Ro312))))))))))))))))))))),kmul(im1,kmadd(rm2,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo31,kmul(nn,nn)),knmsub(kmul(nn,nn),kmadd(rm3,Rojo13,kmul(rm2,Rojo12)),knmsub(ltet1,kmadd(n3,kmul(R4p1213,rm2),kmadd(n3,kmul(R4p1313,rm3),kmadd(nn,kmul(rm2,Ro211),kmadd(nn,kmul(rm3,Ro311),kmadd(rm2,kmsub(n2,R4p1212,kmul(nn,Ro112)),kmul(rm3,kmsub(n2,R4p1213,kmul(nn,Ro113)))))))),kmadd(nn,kmadd(n1,kmul(rm2,ksub(Ro121,Ro211)),kmadd(n2,kmul(rm2,ksub(Ro122,Ro212)),kmadd(n3,kmul(rm2,ksub(Ro123,Ro213)),kmadd(n1,kmul(rm3,ksub(Ro131,Ro311)),kmadd(n2,kmul(rm3,ksub(Ro132,Ro312)),kmul(n3,kmul(rm3,ksub(Ro133,Ro313)))))))),kmadd(ltet3,kmadd(n2,kmul(R4p1223,rm2),kmadd(n2,kmul(R4p1323,rm3),kmadd(nn,kmul(rm2,Ro132),kmadd(nn,kmul(rm3,Ro133),kmadd(rm2,kmsub(n1,R4p1213,kmul(nn,Ro231)),kmul(rm3,kmsub(n1,R4p1313,kmul(nn,Ro331)))))))),kmul(ltet2,kmadd(n1,kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmsub(nn,kmadd(rm2,ksub(Ro122,Ro221),kmul(rm3,ksub(Ro123,Ro321))),kmul(n3,kmadd(R4p1223,rm2,kmul(R4p1323,rm3))))))))))))))); + kmadd(kmadd(n2,kmadd(ltet1,R4p1213,kmul(ltet2,R4p1223)),kmadd(n3,kmadd(ltet1,R4p1313,kmul(ltet2,R4p1323)),kmul(kmadd(ltet2,kmadd(n1,R4p1213,kmul(n2,R4p1223)),kmul(ltet3,kmadd(n1,R4p1313,kmul(n2,R4p1323)))),ToReal(-1.)))),kmadd(im3,rm1,kmul(im1,kmul(rm3,ToReal(-1.)))),kmadd(kmadd(im2,rm1,kmul(im1,kmul(rm2,ToReal(-1.)))),kmadd(ltet1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(ToReal(-1.),kmadd(n1,kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),kmadd(ltet3,kmadd(n2,R4p1223,kmul(n3,R4p1323)),kmul(n3,kmul(kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),ToReal(-1.))))))),kmadd(kmadd(im3,rm2,kmul(im2,kmul(rm3,ToReal(-1.)))),kmadd(ltet1,kmadd(n1,R4p1213,kmul(n2,R4p1223)),kmul(ToReal(-1.),kmadd(n1,kmadd(ltet1,R4p1213,kmul(ltet2,R4p1223)),kmadd(ltet3,kmadd(n1,R4p1323,kmul(n2,R4p2323)),kmul(n3,kmul(kmadd(ltet1,R4p1323,kmul(ltet2,R4p2323)),ToReal(-1.))))))),kmadd(kmul(nn,nn),kmadd(im1,kmadd(rm2,Rojo21,kmadd(rm3,Rojo31,kmul(kmadd(rm2,Rojo12,kmul(rm3,Rojo13)),ToReal(-1.)))),kmadd(im2,kmadd(rm1,kmadd(Rojo21,ToReal(-1.),Rojo12),kmul(rm3,kmadd(Rojo23,ToReal(-1.),Rojo32))),kmul(im3,kmadd(rm1,kmadd(Rojo31,ToReal(-1.),Rojo13),kmul(rm2,kmadd(Rojo32,ToReal(-1.),Rojo23)))))),kmul(nn,kmul(ToReal(-1.),kmadd(im1,kmadd(rm2,kmadd(n2,Ro212,kmadd(n3,Ro213,kmadd(ltet2,Ro221,kmul(ltet3,Ro231)))),kmadd(rm3,kmadd(n2,Ro312,kmadd(n3,Ro313,kmadd(ltet2,Ro321,kmul(ltet3,Ro331)))),kmadd(kmadd(rm2,kmadd(kadd(ltet2,n2),Ro122,kmadd(n3,Ro123,kmul(ltet3,Ro132))),kmul(rm3,kmadd(ltet2,Ro123,kmadd(n2,Ro132,kmul(kadd(ltet3,n3),Ro133))))),ToReal(-1.),kmadd(n1,kmadd(rm2,Ro211,kmadd(rm3,Ro311,kmul(kmadd(rm2,Ro121,kmul(rm3,Ro131)),ToReal(-1.)))),kmul(ltet1,kmadd(rm2,kmadd(Ro112,ToReal(-1.),Ro211),kmul(rm3,kmadd(Ro113,ToReal(-1.),Ro311)))))))),kmadd(im2,kmadd(rm3,kmadd(kadd(ltet2,n2),Ro322,kmadd(n3,Ro323,kmadd(ltet3,Ro332,kmul(kmadd(ltet2,Ro223,kmadd(n2,Ro232,kmul(kadd(ltet3,n3),Ro233))),ToReal(-1.))))),kmadd(ltet1,kmadd(rm1,kmadd(Ro211,ToReal(-1.),Ro112),kmul(rm3,kmadd(Ro213,ToReal(-1.),Ro312))),kmadd(rm1,kmadd(kadd(ltet2,n2),Ro122,kmadd(kmadd(n2,Ro212,kmul(ltet2,Ro221)),ToReal(-1.),kmadd(n3,kmadd(Ro213,ToReal(-1.),Ro123),kmul(ltet3,kmadd(Ro231,ToReal(-1.),Ro132))))),kmul(n1,kmadd(rm1,kmadd(Ro211,ToReal(-1.),Ro121),kmul(rm3,kmadd(Ro231,ToReal(-1.),Ro321))))))),kmul(im3,kmadd(rm1,kmadd(kadd(ltet3,n3),Ro133,kmadd(kmadd(n3,Ro313,kmul(ltet3,Ro331)),ToReal(-1.),kmadd(n1,kmadd(Ro311,ToReal(-1.),Ro131),kmul(n2,kmadd(Ro312,ToReal(-1.),Ro132))))),kmadd(ltet1,kmadd(rm1,kmadd(Ro311,ToReal(-1.),Ro113),kmul(rm2,kmadd(Ro312,ToReal(-1.),Ro213))),kmadd(ltet2,kmadd(rm1,kmadd(Ro321,ToReal(-1.),Ro123),kmul(rm2,kmadd(Ro322,ToReal(-1.),Ro223))),kmul(rm2,kmadd(kadd(ltet3,n3),Ro233,kmadd(kmadd(n3,Ro323,kmul(ltet3,Ro332)),ToReal(-1.),kmadd(n1,kmadd(Ro321,ToReal(-1.),Ro231),kmul(n2,kmadd(Ro322,ToReal(-1.),Ro232)))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi1rL = - kmadd(kmadd(n2,kmul(R4p1212,rm2),kmadd(n3,kmul(R4p1213,rm2),kmadd(n2,kmul(R4p1213,rm3),kmadd(n3,kmul(R4p1313,rm3),kmadd(nn,kmul(rm1,Ro111),kmadd(nn,kmul(rm2,Ro121),kmul(nn,kmul(rm3,Ro131)))))))),kmul(ltet1,ltet1),kmadd(n1,kmul(R4p1313,kmul(rm1,kmul(ltet3,ltet3))),kmadd(n2,kmul(R4p1323,kmul(rm1,kmul(ltet3,ltet3))),kmadd(n2,kmul(R4p2323,kmul(rm2,kmul(ltet3,ltet3))),kmadd(nn,kmul(rm1,kmul(Ro313,kmul(ltet3,ltet3))),kmadd(nn,kmul(rm2,kmul(Ro323,kmul(ltet3,ltet3))),kmadd(n1,kmul(rm1,kmul(Rojo11,kmul(nn,nn))),kmadd(n1,kmul(rm2,kmul(Rojo12,kmul(nn,nn))),kmadd(n1,kmul(rm3,kmul(Rojo13,kmul(nn,nn))),kmadd(n2,kmul(rm1,kmul(Rojo21,kmul(nn,nn))),kmadd(n2,kmul(rm2,kmul(Rojo22,kmul(nn,nn))),kmadd(n2,kmul(rm3,kmul(Rojo23,kmul(nn,nn))),kmadd(n3,kmul(rm1,kmul(Rojo31,kmul(nn,nn))),kmadd(n3,kmul(rm2,kmul(Rojo32,kmul(nn,nn))),kmadd(n3,kmul(rm3,kmul(Rojo33,kmul(nn,nn))),kmadd(rm2,kmsub(n1,kmul(R4p1323,kmul(ltet3,ltet3)),kmul(ltet3,kmul(Rojo32,kmul(nn,nn)))),kmadd(rm3,kmsub(nn,kmul(Ro333,kmul(ltet3,ltet3)),kmul(ltet3,kmul(Rojo33,kmul(nn,nn)))),kmadd(kmul(ltet2,ltet2),kmadd(nn,kmul(rm1,Ro212),kmadd(nn,kmul(rm2,Ro222),kmadd(nn,kmul(rm3,Ro232),kmadd(rm3,kmsub(n3,R4p2323,kmul(n1,R4p1223)),kmul(rm1,kmsub(n1,R4p1212,kmul(n3,R4p1223))))))),kmsub(ltet3,kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro333,ToReal(-2)))),knmsub(rm1,kmul(Rojo31,kmul(nn,nn)),kmadd(n2,kmul(nn,kmsub(rm2,kmul(Ro223,ToReal(-2)),kmul(rm3,kadd(Ro323,Ro233)))),kmul(nn,knmsub(n3,kmul(rm2,kadd(Ro323,Ro233)),kmsub(rm1,kmsub(n1,kmul(Ro113,ToReal(-2)),kmadd(n3,kadd(Ro313,Ro133),kmul(n2,kadd(Ro213,Ro123)))),kmul(n1,kmadd(rm3,kadd(Ro313,Ro133),kmul(rm2,kadd(Ro213,Ro123)))))))))),kmadd(ltet2,kmadd(ltet1,kmadd(n3,kmsub(R4p1213,rm1,kmul(R4p1223,rm2)),kmadd(n2,kmsub(R4p1212,rm1,kmul(R4p1223,rm3)),kmadd(rm3,kmsub(n1,R4p1213,kmul(nn,Ro132)),knmsub(nn,kmul(rm1,kadd(Ro211,Ro112)),kmadd(rm3,kmsub(n3,kmul(R4p1323,ToReal(-2)),kmul(nn,Ro231)),kmul(rm2,kmsub(n1,R4p1212,kmul(nn,kadd(Ro221,Ro122))))))))),kmsub(nn,kmadd(n3,kmul(rm1,Ro132),kmadd(n3,kmul(rm2,Ro232),kmadd(n3,kmul(rm1,Ro312),kmadd(n3,kmul(rm2,Ro322),kmadd(nn,kmul(rm1,Rojo21),kmadd(nn,kmul(rm2,Rojo22),kmadd(nn,kmul(rm3,Rojo23),kmadd(n3,kmul(rm3,kmul(Ro332,ToReal(2))),kmadd(n1,kmadd(rm2,kadd(Ro122,Ro212),kmadd(rm3,kadd(Ro132,Ro312),kmul(rm1,kmul(Ro112,ToReal(2))))),kmul(n2,kmadd(rm1,kadd(Ro122,Ro212),kmadd(rm3,kadd(Ro232,Ro322),kmul(rm2,kmul(Ro222,ToReal(2))))))))))))))),kmul(ltet3,kmadd(nn,kmul(rm1,Ro213),kmadd(nn,kmul(rm2,Ro223),kmadd(nn,kmul(rm3,Ro233),kmadd(nn,kmul(rm1,Ro312),kmadd(nn,kmul(rm2,Ro322),kmadd(nn,kmul(rm3,Ro332),kmadd(rm1,kmsub(n2,R4p1223,kmul(n3,R4p1323)),kmsub(n1,kmadd(R4p1223,rm2,kmsub(R4p1213,kmul(rm1,ToReal(2)),kmul(R4p1323,rm3))),kmul(R4p2323,kmadd(n2,rm3,kmul(n3,rm2)))))))))))))),kmul(ltet1,kmadd(ltet3,kmadd(rm2,kmsub(n1,R4p1213,kmul(nn,Ro123)),kmadd(rm2,kmsub(n3,R4p1323,kmul(nn,Ro321)),kmadd(rm1,kmsub(n3,R4p1313,kmul(nn,kadd(Ro311,Ro113))),kmadd(rm3,kmsub(n1,R4p1313,kmul(nn,kadd(Ro331,Ro133))),kmul(n2,kmadd(R4p1213,rm1,kmadd(R4p1323,rm3,kmul(R4p1223,kmul(rm2,ToReal(2)))))))))),kmul(nn,kmadd(n3,kmul(rm1,Ro131),kmadd(n3,kmul(rm2,Ro231),kmadd(n3,kmul(rm1,Ro311),kmadd(n3,kmul(rm2,Ro321),kmadd(nn,kmul(rm1,Rojo11),kmadd(nn,kmul(rm2,Rojo12),kmadd(nn,kmul(rm3,Rojo13),kmadd(n3,kmul(rm3,kmul(Ro331,ToReal(2))),kmadd(n1,kmadd(rm2,kadd(Ro121,Ro211),kmadd(rm3,kadd(Ro131,Ro311),kmul(rm1,kmul(Ro111,ToReal(2))))),kmul(n2,kmadd(rm1,kadd(Ro121,Ro211),kmadd(rm3,kadd(Ro231,Ro321),kmul(rm2,kmul(Ro221,ToReal(2)))))))))))))))))))))))))))))))))))))); + kmadd(kmul(nn,nn),kmul(ToReal(-1.),kmadd(kmadd(rm1,Rojo11,kmadd(rm2,Rojo12,kmul(rm3,Rojo13))),kmadd(n1,ToReal(-1.),ltet1),kmadd(rm1,kmadd(Rojo21,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo31,kmadd(n3,ToReal(-1.),ltet3))),kmadd(rm2,kmadd(Rojo22,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo32,kmadd(n3,ToReal(-1.),ltet3))),kmul(rm3,kmadd(Rojo23,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo33,kmadd(n3,ToReal(-1.),ltet3)))))))),kmadd(ltet1,kmadd(ltet1,kmadd(kmadd(n2,R4p1212,kmul(n3,R4p1213)),rm2,kmul(kmadd(n2,R4p1213,kmul(n3,R4p1313)),rm3)),kmadd(ltet2,kmadd(n3,kmadd(R4p1223,rm2,kmul(R4p1323,rm3)),kmul(n1,kmul(kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),ToReal(-1.)))),kmadd(ToReal(-1.),kmadd(ltet3,kmadd(kmadd(n1,R4p1213,kmul(n2,R4p1223)),rm2,kmul(kmadd(n1,R4p1313,kmul(n2,R4p1323)),rm3)),kmul(nn,kmul(Ro231,kmadd(n3,rm2,kmul(rm3,kmadd(ltet2,ToReal(-1.),n2)))))),kmul(nn,kmadd(rm1,kmul(Ro111,kmadd(n1,ToReal(-2.),ltet1)),kmadd(rm2,kmul(Ro221,kmadd(n2,ToReal(-2.),ltet2)),kmadd(rm3,kmul(Ro331,kmadd(n3,ToReal(-2.),ltet3)),kmadd(Ro121,kmadd(n2,kmul(rm1,ToReal(-1.)),kmul(rm2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro131,kmadd(n3,kmul(rm1,ToReal(-1.)),kmul(rm3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro211,kmadd(n1,kmul(rm2,ToReal(-1.)),kmul(rm1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro311,kmadd(n1,kmul(rm3,ToReal(-1.)),kmul(rm1,kmadd(n3,ToReal(-1.),ltet3))),kmul(Ro321,kmadd(n2,kmul(rm3,ToReal(-1.)),kmul(rm2,kmadd(n3,ToReal(-1.),ltet3))))))))))))))),kmadd(ltet2,kmadd(ltet1,kmadd(kmadd(n2,R4p1223,kmul(n3,R4p1323)),rm3,kmul(kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(rm1,ToReal(-1.)))),kmadd(ltet3,kmadd(kmadd(n1,R4p1213,kmul(n2,R4p1223)),rm1,kmul(kmadd(n1,R4p1323,kmul(n2,R4p2323)),kmul(rm3,ToReal(-1.)))),kmadd(ltet2,kmadd(rm3,kmadd(n3,R4p2323,kmul(n1,kmul(R4p1223,ToReal(-1.)))),kmul(rm1,kmadd(n1,R4p1212,kmul(n3,kmul(R4p1223,ToReal(-1.)))))),kmul(nn,kmadd(rm1,kmul(Ro112,kmadd(n1,ToReal(-2.),ltet1)),kmadd(rm2,kmul(Ro222,kmadd(n2,ToReal(-2.),ltet2)),kmadd(rm3,kmul(Ro332,kmadd(n3,ToReal(-2.),ltet3)),kmadd(Ro212,kmadd(n1,kmul(rm2,ToReal(-1.)),kmul(rm1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro232,kmadd(n3,kmul(rm2,ToReal(-1.)),kmul(rm3,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro312,kmadd(n1,kmul(rm3,ToReal(-1.)),kmul(rm1,kmadd(n3,ToReal(-1.),ltet3))),kmadd(Ro322,kmadd(n2,kmul(rm3,ToReal(-1.)),kmul(rm2,kmadd(n3,ToReal(-1.),ltet3))),kmul(ToReal(-1.),kmadd(Ro122,kmadd(n2,rm1,kmul(rm2,kmadd(ltet1,ToReal(-1.),n1))),kmul(Ro132,kmadd(n3,rm1,kmul(rm3,kmadd(ltet1,ToReal(-1.),n1))))))))))))))))),kmul(ltet3,kmadd(ltet3,kmadd(kmadd(n1,R4p1313,kmul(n2,R4p1323)),rm1,kmul(kmadd(n1,R4p1323,kmul(n2,R4p2323)),rm2)),kmadd(ltet2,kmadd(n1,kmadd(R4p1213,rm1,kmul(R4p1223,rm2)),kmul(n3,kmul(kmadd(R4p1323,rm1,kmul(R4p2323,rm2)),ToReal(-1.)))),kmadd(nn,kmadd(rm1,kmul(Ro113,kmadd(n1,ToReal(-2.),ltet1)),kmadd(rm2,kmul(Ro223,kmadd(n2,ToReal(-2.),ltet2)),kmadd(rm3,kmul(Ro333,kmadd(n3,ToReal(-2.),ltet3)),kmadd(Ro133,kmul(ToReal(-1.),kmadd(n3,rm1,kmul(rm3,kmadd(ltet1,ToReal(-1.),n1)))),kmadd(Ro213,kmadd(n1,kmul(rm2,ToReal(-1.)),kmul(rm1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro313,kmadd(n1,kmul(rm3,ToReal(-1.)),kmul(rm1,kmadd(n3,ToReal(-1.),ltet3))),kmul(Ro323,kmadd(n2,kmul(rm3,ToReal(-1.)),kmul(rm2,kmadd(n3,ToReal(-1.),ltet3)))))))))),kmul(ToReal(-1.),kmadd(ltet1,kmadd(kmadd(n2,R4p1213,kmul(n3,R4p1313)),rm1,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),rm2)),kmul(nn,kmadd(Ro123,kmadd(n2,rm1,kmul(rm2,kmadd(ltet1,ToReal(-1.),n1))),kmul(Ro233,kmadd(n3,rm2,kmul(rm3,kmadd(ltet2,ToReal(-1.),n2))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi1iL = - kmadd(im1,kmadd(nn,kmul(Ro111,kmul(ltet1,ltet1)),kmadd(n1,kmul(R4p1313,kmul(ltet3,ltet3)),kmadd(n2,kmul(R4p1323,kmul(ltet3,ltet3)),kmadd(nn,kmul(Ro313,kmul(ltet3,ltet3)),kmadd(n1,kmul(Rojo11,kmul(nn,nn)),kmadd(n2,kmul(Rojo21,kmul(nn,nn)),kmadd(n3,kmul(Rojo31,kmul(nn,nn)),kmadd(kmul(ltet2,ltet2),kmadd(n1,R4p1212,kmsub(nn,Ro212,kmul(n3,R4p1223))),kmsub(ltet3,kmsub(nn,kmsub(n1,kmul(Ro113,ToReal(-2)),kmadd(n3,kadd(Ro313,Ro133),kmul(n2,kadd(Ro213,Ro123)))),kmul(Rojo31,kmul(nn,nn))),kmadd(ltet2,kmadd(ltet1,kmadd(n2,R4p1212,kmsub(n3,R4p1213,kmul(nn,kadd(Ro112,Ro211)))),kmsub(nn,kmadd(n3,Ro132,kmadd(n2,kadd(Ro122,Ro212),kmadd(n3,Ro312,kmadd(nn,Rojo21,kmul(n1,kmul(Ro112,ToReal(2))))))),kmul(ltet3,kmadd(n2,R4p1223,kmadd(nn,Ro213,kmadd(nn,Ro312,kmsub(n1,kmul(R4p1213,ToReal(2)),kmul(n3,R4p1323)))))))),kmul(ltet1,kmadd(ltet3,kmadd(n2,R4p1213,kmsub(n3,R4p1313,kmul(nn,kadd(Ro113,Ro311)))),kmul(nn,kmadd(n3,Ro131,kmadd(n2,kadd(Ro121,Ro211),kmadd(n3,Ro311,kmadd(nn,Rojo11,kmul(n1,kmul(Ro111,ToReal(2)))))))))))))))))))),kmadd(im2,kmadd(kmadd(n2,R4p1212,kmadd(n3,R4p1213,kmul(nn,Ro121))),kmul(ltet1,ltet1),kmadd(nn,kmul(Ro222,kmul(ltet2,ltet2)),kmadd(n1,kmul(R4p1323,kmul(ltet3,ltet3)),kmadd(n2,kmul(R4p2323,kmul(ltet3,ltet3)),kmadd(nn,kmul(Ro323,kmul(ltet3,ltet3)),kmadd(n1,kmul(Rojo12,kmul(nn,nn)),kmadd(n2,kmul(Rojo22,kmul(nn,nn)),kmadd(n3,kmul(Rojo32,kmul(nn,nn)),kmadd(ltet3,kmul(nn,kmsub(n2,kmul(Ro223,ToReal(-2)),kmul(n3,Ro323))),knmsub(ltet3,kmadd(nn,kmadd(n1,kadd(Ro213,Ro123),kmul(n3,Ro233)),kmul(Rojo32,kmul(nn,nn))),kmsub(ltet2,kmsub(ltet3,kmadd(n1,R4p1223,kmsub(nn,kadd(Ro223,Ro322),kmul(n3,R4p2323))),kmul(nn,kmadd(n1,kadd(Ro122,Ro212),kmadd(n3,Ro232,kmadd(n3,Ro322,kmadd(nn,Rojo22,kmul(n2,kmul(Ro222,ToReal(2))))))))),kmul(ltet1,kmadd(ltet2,kmsub(n1,R4p1212,kmadd(nn,kadd(Ro122,Ro221),kmul(n3,R4p1223))),kmadd(ltet3,kmadd(n1,R4p1213,kmadd(n3,R4p1323,kmsub(n2,kmul(R4p1223,ToReal(2)),kmul(nn,kadd(Ro321,Ro123))))),kmul(nn,kmadd(n1,kadd(Ro121,Ro211),kmadd(n3,Ro231,kmadd(n3,Ro321,kmadd(nn,Rojo12,kmul(n2,kmul(Ro221,ToReal(2)))))))))))))))))))))),kmul(im3,kmadd(kmadd(n2,R4p1213,kmadd(n3,R4p1313,kmul(nn,Ro131))),kmul(ltet1,ltet1),kmsub(kmul(ltet2,ltet2),kmadd(n3,R4p2323,kmsub(nn,Ro232,kmul(n1,R4p1223))),kmadd(ltet1,kmadd(ltet3,kmadd(n1,R4p1313,kmsub(n2,R4p1323,kmul(nn,kadd(Ro133,Ro331)))),kmadd(ltet2,kmadd(n1,R4p1213,kmsub(n3,kmul(R4p1323,ToReal(-2)),kmadd(nn,kadd(Ro231,Ro132),kmul(n2,R4p1223)))),kmul(nn,kmadd(n1,kadd(Ro131,Ro311),kmadd(n2,kadd(Ro231,Ro321),kmadd(nn,Rojo13,kmul(n3,kmul(Ro331,ToReal(2))))))))),kmadd(nn,knmsub(nn,kmadd(n1,Rojo13,kmadd(n2,Rojo23,kmul(n3,Rojo33))),kmsub(ltet3,kmadd(n1,kadd(Ro133,Ro313),kmadd(n2,kadd(Ro233,Ro323),kmadd(nn,Rojo33,kmul(n3,kmul(Ro333,ToReal(2)))))),kmul(Ro333,kmul(ltet3,ltet3)))),kmul(ltet2,kmadd(ltet3,kmadd(n1,R4p1323,kmsub(n2,R4p2323,kmul(nn,kadd(Ro233,Ro332)))),kmul(nn,kmadd(n1,kadd(Ro132,Ro312),kmadd(n2,kadd(Ro232,Ro322),kmadd(nn,Rojo23,kmul(n3,kmul(Ro332,ToReal(2)))))))))))))))); + kmadd(kmul(nn,nn),kmul(ToReal(-1.),kmadd(im1,kmadd(Rojo11,kmadd(n1,ToReal(-1.),ltet1),kmadd(Rojo21,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo31,kmadd(n3,ToReal(-1.),ltet3)))),kmadd(im2,kmadd(Rojo12,kmadd(n1,ToReal(-1.),ltet1),kmadd(Rojo22,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo32,kmadd(n3,ToReal(-1.),ltet3)))),kmul(im3,kmadd(Rojo13,kmadd(n1,ToReal(-1.),ltet1),kmadd(Rojo23,kmadd(n2,ToReal(-1.),ltet2),kmul(Rojo33,kmadd(n3,ToReal(-1.),ltet3)))))))),kmadd(ltet3,kmadd(im1,kmadd(nn,kmul(Ro113,kmadd(n1,ToReal(-2.),ltet1)),kmadd(R4p1213,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1313,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p1323,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.)))))))),kmadd(im2,kmadd(nn,kmul(Ro223,kmadd(n2,ToReal(-2.),ltet2)),kmadd(R4p1223,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1323,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p2323,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.)))))))),kmul(nn,kmadd(im3,kmul(Ro333,kmadd(n3,ToReal(-2.),ltet3)),kmadd(Ro213,kmul(ToReal(-1.),kmadd(im2,n1,kmul(im1,kmadd(ltet2,ToReal(-1.),n2)))),kmadd(Ro123,kmadd(im1,kmul(n2,ToReal(-1.)),kmul(im2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro133,kmadd(im1,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro233,kmadd(im2,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro313,kmadd(im3,kmul(n1,ToReal(-1.)),kmul(im1,kmadd(n3,ToReal(-1.),ltet3))),kmul(Ro323,kmadd(im3,kmul(n2,ToReal(-1.)),kmul(im2,kmadd(n3,ToReal(-1.),ltet3))))))))))))),kmadd(ltet1,kmadd(im2,kmadd(nn,kmul(Ro221,kmadd(n2,ToReal(-2.),ltet2)),kmadd(R4p1212,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1213,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p1223,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.)))))))),kmadd(im3,kmadd(nn,kmul(Ro331,kmadd(n3,ToReal(-2.),ltet3)),kmadd(R4p1213,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1313,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p1323,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.)))))))),kmul(nn,kmadd(im1,kmul(Ro111,kmadd(n1,ToReal(-2.),ltet1)),kmadd(Ro121,kmadd(im1,kmul(n2,ToReal(-1.)),kmul(im2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro131,kmadd(im1,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro231,kmadd(im2,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n2,ToReal(-1.),ltet2))),kmul(ToReal(-1.),kmadd(Ro211,kmadd(im2,n1,kmul(im1,kmadd(ltet2,ToReal(-1.),n2))),kmadd(Ro311,kmadd(im3,n1,kmul(im1,kmadd(ltet3,ToReal(-1.),n3))),kmul(Ro321,kmadd(im3,n2,kmul(im2,kmadd(ltet3,ToReal(-1.),n3)))))))))))))),kmul(ltet2,kmadd(im3,kmadd(nn,kmul(Ro332,kmadd(n3,ToReal(-2.),ltet3)),kmadd(R4p1223,kmadd(ltet1,n2,kmul(ltet2,kmul(n1,ToReal(-1.)))),kmadd(R4p1323,kmadd(ltet1,n3,kmul(ltet3,kmul(n1,ToReal(-1.)))),kmul(R4p2323,kmadd(ltet2,n3,kmul(ltet3,kmul(n2,ToReal(-1.)))))))),kmadd(im1,kmadd(nn,kmul(Ro112,kmadd(n1,ToReal(-2.),ltet1)),kmadd(R4p1212,kmadd(ltet2,n1,kmul(ltet1,kmul(n2,ToReal(-1.)))),kmadd(R4p1213,kmadd(ltet3,n1,kmul(ltet1,kmul(n3,ToReal(-1.)))),kmul(R4p1223,kmadd(ltet3,n2,kmul(ltet2,kmul(n3,ToReal(-1.)))))))),kmul(nn,kmadd(im2,kmul(Ro222,kmadd(n2,ToReal(-2.),ltet2)),kmadd(Ro122,kmadd(im1,kmul(n2,ToReal(-1.)),kmul(im2,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro132,kmadd(im1,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n1,ToReal(-1.),ltet1))),kmadd(Ro212,kmadd(im2,kmul(n1,ToReal(-1.)),kmul(im1,kmadd(n2,ToReal(-1.),ltet2))),kmadd(Ro232,kmadd(im2,kmul(n3,ToReal(-1.)),kmul(im3,kmadd(n2,ToReal(-1.),ltet2))),kmul(ToReal(-1.),kmadd(Ro312,kmadd(im3,n1,kmul(im1,kmadd(ltet3,ToReal(-1.),n3))),kmul(Ro322,kmadd(im3,n2,kmul(im2,kmadd(ltet3,ToReal(-1.),n3)))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi0rL = - kmadd(rm1,kmul(rm2,kmul(Rojo12,kmul(nn,nn))),kmadd(rm1,kmul(rm3,kmul(Rojo13,kmul(nn,nn))),kmadd(rm1,kmul(rm2,kmul(Rojo21,kmul(nn,nn))),kmadd(rm2,kmul(rm3,kmul(Rojo23,kmul(nn,nn))),kmadd(rm1,kmul(rm3,kmul(Rojo31,kmul(nn,nn))),kmadd(rm2,kmul(rm3,kmul(Rojo32,kmul(nn,nn))),kmadd(R4p1212,kmul(kmul(ltet2,ltet2),kmul(rm1,rm1)),kmadd(Rojo11,kmul(kmul(nn,nn),kmul(rm1,rm1)),kmadd(R4p1212,kmul(kmul(ltet1,ltet1),kmul(rm2,rm2)),kmadd(Rojo22,kmul(kmul(nn,nn),kmul(rm2,rm2)),kmadd(R4p1313,kmul(kmul(ltet1,ltet1),kmul(rm3,rm3)),kmadd(R4p2323,kmul(kmul(ltet2,ltet2),kmul(rm3,rm3)),kmadd(kmadd(ltet1,kmul(ltet2,kmul(R4p1212,kmul(rm1,rm2))),kmadd(ltet1,kmul(ltet3,kmul(R4p1213,kmul(rm1,rm2))),kmadd(ltet1,kmul(ltet3,kmul(R4p1313,kmul(rm1,rm3))),kmadd(ltet1,kmul(ltet3,kmul(R4p1323,kmul(rm2,rm3))),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,kmul(rm2,rm3))),kmadd(ltet1,kmul(ltet2,kmul(R4p1323,kmul(im3,im3))),kmadd(ltet1,kmul(nn,kmul(Ro313,kmul(im3,im3))),kmadd(ltet3,kmul(nn,kmul(Ro333,kmul(im3,im3))),kmul(R4p1223,kmul(rm1,kmul(rm3,kmul(ltet2,ltet2)))))))))))),ToReal(-2),kmadd(ltet1,kmul(ltet3,kmul(R4p1223,kmul(kmul(rm2,rm2),ToReal(-2)))),kmadd(R4p1313,kmsub(kmul(ltet3,ltet3),kmul(rm1,rm1),kmul(kmul(im3,im3),kmul(ltet1,ltet1))),kmadd(R4p2323,kmsub(kmul(ltet3,ltet3),kmul(rm2,rm2),kmul(kmul(im3,im3),kmul(ltet2,ltet2))),kmadd(Rojo33,kmsub(kmul(nn,nn),kmul(rm3,rm3),kmul(kmul(im3,im3),kmul(nn,nn))),kmadd(kmadd(ltet1,kmul(nn,kmul(rm1,kmul(rm2,Ro112))),kmadd(ltet1,kmul(nn,kmul(rm1,kmul(rm3,Ro113))),kmadd(ltet2,kmul(nn,kmul(rm1,kmul(rm3,Ro123))),kmadd(ltet3,kmul(nn,kmul(rm1,kmul(rm2,Ro132))),kmadd(ltet3,kmul(nn,kmul(rm1,kmul(rm3,Ro133))),kmadd(ltet1,kmul(nn,kmul(rm1,kmul(rm2,Ro211))),kmadd(ltet1,kmul(nn,kmul(rm2,kmul(rm3,Ro213))),kmadd(ltet2,kmul(nn,kmul(rm2,kmul(rm3,Ro223))),kmadd(ltet3,kmul(nn,kmul(rm1,kmul(rm2,Ro231))),kmadd(ltet3,kmul(nn,kmul(rm2,kmul(rm3,Ro233))),kmadd(ltet1,kmul(nn,kmul(rm1,kmul(rm3,Ro311))),kmadd(ltet1,kmul(nn,kmul(rm2,kmul(rm3,Ro312))),kmadd(ltet2,kmul(nn,kmul(rm2,kmul(rm3,Ro322))),kmadd(ltet3,kmul(nn,kmul(rm1,kmul(rm3,Ro331))),kmadd(ltet3,kmul(nn,kmul(rm2,kmul(rm3,Ro332))),kmadd(R4p1213,kmul(rm2,kmul(rm3,kmul(ltet1,ltet1))),kmul(R4p1323,kmul(rm1,kmul(rm2,kmul(ltet3,ltet3)))))))))))))))))))),ToReal(2),kmadd(ltet1,kmul(nn,kmul(Ro111,kmul(kmul(rm1,rm1),ToReal(2)))),kmadd(ltet2,kmul(nn,kmul(Ro121,kmul(kmul(rm1,rm1),ToReal(2)))),kmadd(ltet3,kmul(nn,kmul(Ro131,kmul(kmul(rm1,rm1),ToReal(2)))),kmadd(ltet1,kmul(nn,kmul(Ro212,kmul(kmul(rm2,rm2),ToReal(2)))),kmadd(ltet2,kmul(nn,kmul(Ro222,kmul(kmul(rm2,rm2),ToReal(2)))),kmadd(ltet3,kmul(nn,kmul(Ro232,kmul(kmul(rm2,rm2),ToReal(2)))),kmadd(ltet1,kmul(ltet2,kmul(R4p1323,kmul(kmul(rm3,rm3),ToReal(2)))),kmadd(ltet1,kmul(nn,kmul(Ro313,kmul(kmul(rm3,rm3),ToReal(2)))),kmadd(ltet2,kmul(nn,kmul(Ro323,kmul(kmul(rm3,rm3),ToReal(2)))),kmadd(ltet3,kmul(nn,kmul(Ro333,kmul(kmul(rm3,rm3),ToReal(2)))),knmsub(kmul(im2,im2),kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmadd(Rojo22,kmul(nn,nn),kmadd(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2))),kmadd(ltet1,kmul(nn,kmul(Ro212,ToReal(2))),kmadd(ltet2,kmul(nn,kmul(Ro222,ToReal(2))),kmul(ltet3,kmul(nn,kmul(Ro232,ToReal(2)))))))))),kmadd(ltet2,kmadd(ltet1,kmul(R4p1213,kmul(rm1,kmul(rm3,ToReal(-2)))),kmadd(ltet3,kmul(R4p1323,kmul(rm1,kmul(rm3,ToReal(-2)))),kmadd(nn,kmul(Ro323,kmul(kmul(im3,im3),ToReal(-2))),kmadd(ltet3,kmul(R4p1223,kmul(rm1,kmul(rm2,ToReal(2)))),kmadd(ltet1,kmul(R4p1223,kmul(rm2,kmul(rm3,ToReal(2)))),kmadd(nn,kmul(rm1,kmul(rm2,kmul(Ro122,ToReal(2)))),kmadd(nn,kmul(rm1,kmul(rm2,kmul(Ro221,ToReal(2)))),kmadd(nn,kmul(rm1,kmul(rm3,kmul(Ro321,ToReal(2)))),kmul(ltet3,kmul(R4p1213,kmul(kmul(rm1,rm1),ToReal(2)))))))))))),knmsub(kmul(im1,im1),kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmadd(ltet2,kmul(kmadd(ltet3,R4p1213,kmul(nn,Ro121)),ToReal(2)),kmadd(ltet3,kmul(nn,kmul(Ro131,ToReal(2))),kmul(nn,kmadd(nn,Rojo11,kmul(ltet1,kmul(Ro111,ToReal(2))))))))),kmsub(im1,kmul(im3,kmadd(R4p1223,kmul(kmul(ltet2,ltet2),ToReal(2)),kmadd(ltet1,kmul(kmadd(ltet2,R4p1213,kmsub(ltet3,R4p1313,kmul(nn,kadd(Ro113,Ro311)))),ToReal(2)),kmsub(ltet2,kmul(kmsub(ltet3,R4p1323,kmul(nn,kadd(Ro123,Ro321))),ToReal(2)),kmul(nn,kmadd(nn,kadd(Rojo13,Rojo31),kmul(ltet3,kmul(kadd(Ro133,Ro331),ToReal(2))))))))),kmul(im2,kmadd(im3,kmadd(R4p1213,kmul(kmul(ltet1,ltet1),ToReal(2)),kmadd(ltet1,kmul(kmadd(ltet2,R4p1223,kmsub(nn,kadd(Ro213,Ro312),kmul(ltet3,R4p1323))),ToReal(2)),kmadd(ltet2,kmadd(ltet3,kmul(R4p2323,ToReal(-2)),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2)))),kmul(nn,kmadd(nn,kadd(Rojo23,Rojo32),kmul(ltet3,kmul(kadd(Ro233,Ro332),ToReal(2)))))))),kmul(im1,kmadd(Rojo12,kmul(nn,nn),kmadd(Rojo21,kmul(nn,nn),kmadd(ltet1,kmul(ToReal(-2),kmadd(ltet2,R4p1212,kmsub(ltet3,R4p1213,kmul(nn,kadd(Ro112,Ro211))))),kmadd(ltet3,kmul(nn,kmul(Ro132,ToReal(2))),kmadd(ltet2,kmul(kmadd(ltet3,R4p1223,kmul(nn,kadd(Ro122,Ro221))),ToReal(2)),kmadd(ltet3,kmul(nn,kmul(Ro231,ToReal(2))),kmul(R4p1323,kmul(kmul(ltet3,ltet3),ToReal(2)))))))))))))))))))))))))))))))))))))))))))); + kmadd(kmadd(kmadd(ltet1,kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),kmul(ltet3,kmul(kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),ToReal(-1.)))),kmadd(im1,im2,kmul(rm1,kmul(rm2,ToReal(-1.)))),kmadd(kmadd(ltet1,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)),kmadd(ltet2,kmul(ltet3,R4p1323),kmul(R4p1223,kmul(ltet2,ltet2)))),kmadd(im1,im3,kmul(rm1,kmul(rm3,ToReal(-1.)))),kmadd(kmadd(rm2,rm3,kmul(im2,kmul(im3,ToReal(-1.)))),kmadd(R4p1213,kmul(ltet1,ltet1),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,ToReal(-1.))),kmul(ltet1,kmadd(ltet2,R4p1223,kmul(ltet3,kmul(R4p1323,ToReal(-1.))))))),kmul(nn,kmadd(kmadd(ltet1,kadd(Ro112,Ro211),kmadd(ltet2,kadd(Ro122,Ro221),kmul(ltet3,kadd(Ro132,Ro231)))),kmadd(rm1,rm2,kmul(im1,kmul(im2,ToReal(-1.)))),kmadd(kmadd(ltet1,kadd(Ro113,Ro311),kmadd(ltet2,kadd(Ro123,Ro321),kmul(ltet3,kadd(Ro133,Ro331)))),kmadd(rm1,rm3,kmul(im1,kmul(im3,ToReal(-1.)))),kmadd(kmadd(ltet1,kadd(Ro213,Ro312),kmadd(ltet2,kadd(Ro223,Ro322),kmul(ltet3,kadd(Ro233,Ro332)))),kmadd(rm2,rm3,kmul(im2,kmul(im3,ToReal(-1.)))),kmadd(kmadd(ltet1,Ro111,kmadd(ltet2,Ro121,kmul(ltet3,Ro131))),kmadd(rm1,rm1,kmul(kmul(im1,im1),ToReal(-1.))),kmadd(kmadd(ltet1,Ro212,kmadd(ltet2,Ro222,kmul(ltet3,Ro232))),kmadd(rm2,rm2,kmul(kmul(im2,im2),ToReal(-1.))),kmul(kmadd(ltet1,Ro313,kmadd(ltet2,Ro323,kmul(ltet3,Ro333))),kmadd(rm3,rm3,kmul(kmul(im3,im3),ToReal(-1.))))))))))))),ToReal(2.),kmul(ToReal(-1.),kmadd(kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmul(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2.)))))),kmadd(im2,im2,kmul(kmul(rm2,rm2),ToReal(-1.))),kmadd(kmul(nn,nn),kmadd(im1,kmadd(im2,kadd(Rojo12,Rojo21),kmul(im3,kadd(Rojo13,Rojo31))),kmadd(im2,kmul(im3,Rojo32),kmadd(Rojo11,kmul(im1,im1),kmadd(kmadd(rm1,kmadd(rm2,kadd(Rojo12,Rojo21),kmul(rm3,kadd(Rojo13,Rojo31))),kmadd(rm2,kmul(rm3,Rojo32),kmul(Rojo11,kmul(rm1,rm1)))),ToReal(-1.),kmadd(Rojo23,kmadd(im2,im3,kmul(rm2,kmul(rm3,ToReal(-1.)))),kmadd(Rojo22,kmadd(im2,im2,kmul(kmul(rm2,rm2),ToReal(-1.))),kmul(Rojo33,kmadd(im3,im3,kmul(kmul(rm3,rm3),ToReal(-1.)))))))))),kmadd(kmadd(im1,im1,kmul(kmul(rm1,rm1),ToReal(-1.))),kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmul(ltet2,kmul(ltet3,kmul(R4p1213,ToReal(2.)))))),kmul(kmadd(im3,im3,kmul(kmul(rm3,rm3),ToReal(-1.))),kmadd(R4p1313,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet2,ltet2),kmul(ltet1,kmul(ltet2,kmul(R4p1323,ToReal(2.)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi0iL = - kmadd(im1,kmadd(nn,kmadd(ltet1,kmul(ToReal(2),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2)))))),kmul(nn,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2))))))),kmul(ToReal(2),kmadd(kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmul(ltet3,ltet3),kmadd(kmul(ltet2,ltet2),kmsub(R4p1212,rm1,kmul(R4p1223,rm3)),kmadd(ltet2,knmsub(ltet1,kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmadd(ltet3,kmadd(R4p1223,rm2,kmsub(R4p1213,kmul(rm1,ToReal(2)),kmul(R4p1323,rm3))),kmul(nn,kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2)))))))),kmul(ltet3,kmsub(nn,kmadd(rm2,kadd(Ro132,Ro231),kmadd(rm3,kadd(Ro133,Ro331),kmul(rm1,kmul(Ro131,ToReal(2))))),kmul(ltet1,kmadd(R4p1213,rm2,kmul(R4p1313,rm3)))))))))),kmadd(im3,kmadd(kmul(ltet2,ltet2),kmadd(R4p1223,kmul(rm1,ToReal(-2)),kmul(R4p2323,kmul(rm3,ToReal(2)))),kmadd(nn,kmadd(ltet3,kmul(ToReal(2),kmadd(rm1,kadd(Ro133,Ro331),kmadd(rm2,kadd(Ro233,Ro332),kmul(rm3,kmul(Ro333,ToReal(2)))))),kmul(nn,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,ToReal(2))))))),kmul(ToReal(2),kmadd(kmadd(R4p1213,rm2,kmul(R4p1313,rm3)),kmul(ltet1,ltet1),kmadd(ltet1,knmsub(ltet3,kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmadd(ltet2,kmadd(R4p1223,rm2,kmsub(R4p1323,kmul(rm3,ToReal(2)),kmul(R4p1213,rm1))),kmul(nn,kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2)))))))),kmul(ltet2,kmsub(nn,kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2))))),kmul(ltet3,kmadd(R4p1323,rm1,kmul(R4p2323,rm2)))))))))),kmul(im2,kmadd(kmadd(rm1,kadd(Rojo12,Rojo21),kmul(rm3,kadd(Rojo23,Rojo32))),kmul(nn,nn),kmadd(ToReal(2),kmadd(ltet3,kmul(nn,kmadd(rm1,kadd(Ro132,Ro231),kmul(rm3,Ro233))),kmadd(kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmul(ltet1,ltet1),kmadd(kmadd(R4p1323,rm1,kmul(R4p2323,rm2)),kmul(ltet3,ltet3),kmadd(rm2,kmul(Rojo22,kmul(nn,nn)),kmadd(ltet1,kmadd(ltet2,kmsub(R4p1223,rm3,kmul(R4p1212,rm1)),kmsub(nn,kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2))))),kmul(ltet3,kmadd(R4p1213,rm1,kmadd(R4p1323,rm3,kmul(R4p1223,kmul(rm2,ToReal(2)))))))),kmul(ltet2,kmadd(ltet3,kmsub(R4p1223,rm1,kmul(R4p2323,rm3)),kmul(nn,kmadd(rm1,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro223,Ro322),kmul(rm2,kmul(Ro222,ToReal(2))))))))))))),kmul(ltet3,kmul(nn,kmadd(rm3,kmul(Ro332,ToReal(2)),kmul(rm2,kmul(Ro232,ToReal(4))))))))))); + kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmul(kmadd(ltet1,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)),kmadd(ltet2,kmul(ltet3,R4p1323),kmul(R4p1223,kmul(ltet2,ltet2)))),ToReal(-2.)),kmadd(ToReal(2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmul(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet3,kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),kmul(ltet1,kmul(kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),ToReal(-1.)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(ltet1,ltet1),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,ToReal(-1.))),kmul(ltet1,kmadd(ltet2,R4p1223,kmul(ltet3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmul(ltet2,kmul(ltet3,kmul(R4p1213,ToReal(2.))))))),kmadd(im3,kmul(rm3,kmadd(R4p1313,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet2,ltet2),kmul(ltet1,kmul(ltet2,kmul(R4p1323,ToReal(2.))))))),kmul(nn,kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet1,kadd(Ro112,Ro211),kmadd(ltet2,kadd(Ro122,Ro221),kmul(ltet3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(ltet1,kadd(Ro113,Ro311),kmadd(ltet2,kadd(Ro123,Ro321),kmul(ltet3,kadd(Ro133,Ro331)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(ltet1,kadd(Ro213,Ro312),kmadd(ltet2,kadd(Ro223,Ro322),kmul(ltet3,kadd(Ro233,Ro332)))),kmul(kmadd(im1,kmul(rm1,kmadd(ltet1,Ro111,kmadd(ltet2,Ro121,kmul(ltet3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(ltet1,Ro212,kmadd(ltet2,Ro222,kmul(ltet3,Ro232)))),kmul(im3,kmul(rm3,kmadd(ltet1,Ro313,kmadd(ltet2,Ro323,kmul(ltet3,Ro333))))))),ToReal(2.))))))))))),kmul(kmul(nn,nn),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,ToReal(2.))))))))))); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); @@ -1367,7 +1379,7 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int vec_store_nta_partial(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } - LC_ENDLOOP3VEC(WeylScal4_psis_calc_Nth); + CCTK_ENDLOOP3STR(WeylScal4_psis_calc_Nth); } extern "C" void WeylScal4_psis_calc_Nth(CCTK_ARGUMENTS) |