diff options
Diffstat (limited to 'src/WeylScal4_psi4_calc_2nd.cc')
| -rw-r--r-- | src/WeylScal4_psi4_calc_2nd.cc | 745 |
1 files changed, 392 insertions, 353 deletions
diff --git a/src/WeylScal4_psi4_calc_2nd.cc b/src/WeylScal4_psi4_calc_2nd.cc index ddf5a89..eed89f8 100644 --- a/src/WeylScal4_psi4_calc_2nd.cc +++ b/src/WeylScal4_psi4_calc_2nd.cc @@ -32,7 +32,9 @@ extern "C" void WeylScal4_psi4_calc_2nd_SelectBCs(CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; - CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; + if (cctk_iteration % WeylScal4_psi4_calc_2nd_calc_every != WeylScal4_psi4_calc_2nd_calc_offset) + return; + CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); @@ -42,7 +44,7 @@ extern "C" void WeylScal4_psi4_calc_2nd_SelectBCs(CCTK_ARGUMENTS) return; } -static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const imin[3], int const imax[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) +static void WeylScal4_psi4_calc_2nd_Body(const cGH* restrict const cctkGH, const int dir, const int face, const CCTK_REAL normal[3], const CCTK_REAL tangentA[3], const CCTK_REAL tangentB[3], const int imin[3], const int imax[3], const int n_subblock_gfs, CCTK_REAL* restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; @@ -51,111 +53,133 @@ static void WeylScal4_psi4_calc_2nd_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); + const ptrdiff_t di CCTK_ATTRIBUTE_UNUSED = 1; + const ptrdiff_t dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); + const ptrdiff_t dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); + const ptrdiff_t cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di; + const ptrdiff_t cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj; + const ptrdiff_t cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk; + const CCTK_REAL_VEC dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0)); + const CCTK_REAL_VEC dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1)); + const CCTK_REAL_VEC dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2)); + const CCTK_REAL_VEC dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME); + const CCTK_REAL_VEC t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time); + const CCTK_REAL_VEC dxi CCTK_ATTRIBUTE_UNUSED = INV(dx); + const CCTK_REAL_VEC dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); + const CCTK_REAL_VEC dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); + const CCTK_REAL_VEC khalf CCTK_ATTRIBUTE_UNUSED = ToReal(0.5); + const CCTK_REAL_VEC kthird CCTK_ATTRIBUTE_UNUSED = ToReal(1.0/3.0); + const CCTK_REAL_VEC ktwothird CCTK_ATTRIBUTE_UNUSED = ToReal(2.0/3.0); + const CCTK_REAL_VEC kfourthird CCTK_ATTRIBUTE_UNUSED = ToReal(4.0/3.0); + const CCTK_REAL_VEC keightthird CCTK_ATTRIBUTE_UNUSED = ToReal(8.0/3.0); + const CCTK_REAL_VEC hdxi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dxi); + const CCTK_REAL_VEC hdyi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dyi); + const CCTK_REAL_VEC 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)); + const CCTK_REAL_VEC p1o12dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dx); + const CCTK_REAL_VEC p1o12dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dy); + const CCTK_REAL_VEC p1o12dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dz); + const CCTK_REAL_VEC p1o144dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dy,dx)); + const CCTK_REAL_VEC p1o144dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dx)); + const CCTK_REAL_VEC p1o144dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dy)); + const CCTK_REAL_VEC p1o180dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dx,dx)); + const CCTK_REAL_VEC p1o180dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dy,dy)); + const CCTK_REAL_VEC p1o180dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dz,dz)); + const CCTK_REAL_VEC p1o2dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dx); + const CCTK_REAL_VEC p1o2dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dy); + const CCTK_REAL_VEC p1o2dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.5),dz); + const CCTK_REAL_VEC p1o3600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dy,dx)); + const CCTK_REAL_VEC p1o3600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dx)); + const CCTK_REAL_VEC p1o3600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dy)); + const CCTK_REAL_VEC p1o4dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dy,dx)); + const CCTK_REAL_VEC p1o4dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dx)); + const CCTK_REAL_VEC p1o4dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dy)); + const CCTK_REAL_VEC p1o5040dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dx,dx)); + const CCTK_REAL_VEC p1o5040dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dy,dy)); + const CCTK_REAL_VEC p1o5040dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dz,dz)); + const CCTK_REAL_VEC p1o60dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dx); + const CCTK_REAL_VEC p1o60dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dy); + const CCTK_REAL_VEC p1o60dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0166666666666666666666666666667),dz); + const CCTK_REAL_VEC p1o705600dxdy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dy,dx)); + const CCTK_REAL_VEC p1o705600dxdz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dx)); + const CCTK_REAL_VEC p1o705600dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dy)); + const CCTK_REAL_VEC p1o840dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dx); + const CCTK_REAL_VEC p1o840dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dy); + const CCTK_REAL_VEC p1o840dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00119047619047619047619047619048),dz); + const CCTK_REAL_VEC p1odx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dx,dx)); + const CCTK_REAL_VEC p1ody2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dy,dy)); + const CCTK_REAL_VEC p1odz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dz,dz)); + const CCTK_REAL_VEC pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx)); + const CCTK_REAL_VEC pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy)); + const CCTK_REAL_VEC pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz)); /* Jacobian variable pointers */ - 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) + const bool use_jacobian1 = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) + && strlen(jacobian_group) > 0; + const bool use_jacobian = assume_use_jacobian>=0 ? assume_use_jacobian : use_jacobian1; + const bool usejacobian CCTK_ATTRIBUTE_UNUSED = use_jacobian; + if (use_jacobian && (strlen(jacobian_determinant_group) == 0 || strlen(jacobian_inverse_group) == 0 || strlen(jacobian_derivative_group) == 0)) { - CCTK_WARN (1, "GenericFD::jacobian_group and GenericFD::jacobian_derivative_group must both be set to valid group names"); + CCTK_WARN(1, "GenericFD::jacobian_group, GenericFD::jacobian_determinant_group, GenericFD::jacobian_inverse_group, and GenericFD::jacobian_derivative_group must all be set to valid group names"); } - CCTK_REAL const *restrict jacobian_ptrs[9]; + const CCTK_REAL* restrict jacobian_ptrs[9]; if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_group, 9, jacobian_ptrs); - CCTK_REAL const *restrict const J11 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[0] : 0; - CCTK_REAL const *restrict const J12 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[1] : 0; - CCTK_REAL const *restrict const J13 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[2] : 0; - CCTK_REAL const *restrict const J21 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[3] : 0; - CCTK_REAL const *restrict const J22 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[4] : 0; - CCTK_REAL const *restrict const J23 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[5] : 0; - CCTK_REAL const *restrict const J31 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[6] : 0; - CCTK_REAL const *restrict const J32 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[7] : 0; - CCTK_REAL const *restrict const J33 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[8] : 0; + const CCTK_REAL* restrict const J11 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[0] : 0; + const CCTK_REAL* restrict const J12 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[1] : 0; + const CCTK_REAL* restrict const J13 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[2] : 0; + const CCTK_REAL* restrict const J21 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[3] : 0; + const CCTK_REAL* restrict const J22 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[4] : 0; + const CCTK_REAL* restrict const J23 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[5] : 0; + const CCTK_REAL* restrict const J31 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[6] : 0; + const CCTK_REAL* restrict const J32 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[7] : 0; + const CCTK_REAL* restrict const J33 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[8] : 0; + + const CCTK_REAL* restrict jacobian_determinant_ptrs[1] CCTK_ATTRIBUTE_UNUSED; + if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_determinant_group, + 1, jacobian_determinant_ptrs); + + const CCTK_REAL* restrict const detJ CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_ptrs[0] : 0; - CCTK_REAL const *restrict jacobian_derivative_ptrs[18] CCTK_ATTRIBUTE_UNUSED; + const CCTK_REAL* restrict jacobian_inverse_ptrs[9] CCTK_ATTRIBUTE_UNUSED; + if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_inverse_group, + 9, jacobian_inverse_ptrs); + + const CCTK_REAL* restrict const iJ11 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[0] : 0; + const CCTK_REAL* restrict const iJ12 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[1] : 0; + const CCTK_REAL* restrict const iJ13 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[2] : 0; + const CCTK_REAL* restrict const iJ21 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[3] : 0; + const CCTK_REAL* restrict const iJ22 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[4] : 0; + const CCTK_REAL* restrict const iJ23 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[5] : 0; + const CCTK_REAL* restrict const iJ31 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[6] : 0; + const CCTK_REAL* restrict const iJ32 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[7] : 0; + const CCTK_REAL* restrict const iJ33 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_inverse_ptrs[8] : 0; + + const CCTK_REAL* restrict jacobian_derivative_ptrs[18] CCTK_ATTRIBUTE_UNUSED; if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_derivative_group, 18, jacobian_derivative_ptrs); - CCTK_REAL const *restrict const dJ111 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[0] : 0; - CCTK_REAL const *restrict const dJ112 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[1] : 0; - CCTK_REAL const *restrict const dJ113 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[2] : 0; - CCTK_REAL const *restrict const dJ122 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[3] : 0; - CCTK_REAL const *restrict const dJ123 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[4] : 0; - CCTK_REAL const *restrict const dJ133 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[5] : 0; - CCTK_REAL const *restrict const dJ211 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[6] : 0; - CCTK_REAL const *restrict const dJ212 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[7] : 0; - CCTK_REAL const *restrict const dJ213 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[8] : 0; - CCTK_REAL const *restrict const dJ222 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[9] : 0; - CCTK_REAL const *restrict const dJ223 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[10] : 0; - CCTK_REAL const *restrict const dJ233 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[11] : 0; - CCTK_REAL const *restrict const dJ311 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[12] : 0; - CCTK_REAL const *restrict const dJ312 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[13] : 0; - CCTK_REAL const *restrict const dJ313 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[14] : 0; - CCTK_REAL const *restrict const dJ322 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[15] : 0; - CCTK_REAL const *restrict const dJ323 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[16] : 0; - CCTK_REAL const *restrict const dJ333 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[17] : 0; + const CCTK_REAL* restrict const dJ111 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[0] : 0; + const CCTK_REAL* restrict const dJ112 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[1] : 0; + const CCTK_REAL* restrict const dJ113 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[2] : 0; + const CCTK_REAL* restrict const dJ122 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[3] : 0; + const CCTK_REAL* restrict const dJ123 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[4] : 0; + const CCTK_REAL* restrict const dJ133 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[5] : 0; + const CCTK_REAL* restrict const dJ211 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[6] : 0; + const CCTK_REAL* restrict const dJ212 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[7] : 0; + const CCTK_REAL* restrict const dJ213 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[8] : 0; + const CCTK_REAL* restrict const dJ222 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[9] : 0; + const CCTK_REAL* restrict const dJ223 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[10] : 0; + const CCTK_REAL* restrict const dJ233 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[11] : 0; + const CCTK_REAL* restrict const dJ311 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[12] : 0; + const CCTK_REAL* restrict const dJ312 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[13] : 0; + const CCTK_REAL* restrict const dJ313 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[14] : 0; + const CCTK_REAL* restrict const dJ322 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[15] : 0; + const CCTK_REAL* restrict const dJ323 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[16] : 0; + const CCTK_REAL* restrict const dJ333 CCTK_ATTRIBUTE_UNUSED = use_jacobian ? jacobian_derivative_ptrs[17] : 0; /* Assign local copies of arrays functions */ @@ -166,13 +190,20 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int /* Copy local copies back to grid functions */ /* Loop over the grid points */ - #pragma omp parallel + const int imin0=imin[0]; + const int imin1=imin[1]; + const int imin2=imin[2]; + const int imax0=imax[0]; + const int imax1=imax[1]; + const int imax2=imax[2]; + #pragma omp parallel // reduction(+: vec_iter_counter, vec_op_counter, vec_mem_counter) CCTK_LOOP3STR(WeylScal4_psi4_calc_2nd, - i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], + i,j,k, imin0,imin1,imin2, imax0,imax1,imax2, cctk_ash[0],cctk_ash[1],cctk_ash[2], - kimin,kimax, CCTK_REAL_VEC_SIZE) + vecimin,vecimax, CCTK_REAL_VEC_SIZE) { - ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; + const ptrdiff_t index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; + // vec_iter_counter+=CCTK_REAL_VEC_SIZE; /* Assign local copies of grid functions */ @@ -229,82 +260,83 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int /* Include user supplied include files */ /* Precompute derivatives */ - CCTK_REAL_VEC PDstandard2nd1gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd11gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd22gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd33gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd12gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd13gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd23gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd1kzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd2kzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC PDstandard2nd3kzz CCTK_ATTRIBUTE_UNUSED ; - - switch(fdOrder) + CCTK_REAL_VEC PDstandard2nd1gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd11gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd22gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd33gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd12gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd13gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd23gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd1kzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd2kzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC PDstandard2nd3kzz CCTK_ATTRIBUTE_UNUSED; + + switch (fdOrder) { case 2: + { PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); @@ -378,8 +410,10 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; + } case 4: + { PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); @@ -453,8 +487,10 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; + } case 6: + { PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); @@ -528,8 +564,10 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; + } case 8: + { PDstandard2nd1gxx = PDstandard2nd1(&gxx[index]); PDstandard2nd2gxx = PDstandard2nd2(&gxx[index]); PDstandard2nd3gxx = PDstandard2nd3(&gxx[index]); @@ -603,66 +641,69 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int PDstandard2nd2kzz = PDstandard2nd2(&kzz[index]); PDstandard2nd3kzz = PDstandard2nd3(&kzz[index]); break; + } + default: + CCTK_BUILTIN_UNREACHABLE(); } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC JacPDstandard2nd11gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd11gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd11gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd12gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd12gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd12gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd12gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd13gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1kyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1kyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd1kzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd21gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd22gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd22gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd22gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd23gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd23gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd23gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd23gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2kxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2kyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd2kzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd31gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd31gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd31gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd31gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd32gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd33gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd33gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd33gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gyz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3gzz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3kxx CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3kxy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3kxz CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3kyy CCTK_ATTRIBUTE_UNUSED ; - CCTK_REAL_VEC JacPDstandard2nd3kyz CCTK_ATTRIBUTE_UNUSED ; + CCTK_REAL_VEC JacPDstandard2nd11gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd11gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd11gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd12gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd12gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd12gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd12gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd13gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1kyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1kyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd1kzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd21gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd22gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd22gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd22gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd23gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd23gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd23gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd23gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2kxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2kyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd2kzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd31gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd31gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd31gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd31gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd32gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd33gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd33gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd33gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gyz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3gzz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3kxx CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3kxy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3kxz CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3kyy CCTK_ATTRIBUTE_UNUSED; + CCTK_REAL_VEC JacPDstandard2nd3kyz CCTK_ATTRIBUTE_UNUSED; if (use_jacobian) { @@ -766,31 +807,31 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int kmadd(J13L,PDstandard2nd1kyz,kmadd(J23L,PDstandard2nd2kyz,kmul(J33L,PDstandard2nd3kyz))); JacPDstandard2nd11gyy = - kmadd(dJ111L,PDstandard2nd1gyy,kmadd(dJ211L,PDstandard2nd2gyy,kmadd(dJ311L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,kmul(J11L,J11L),kmadd(PDstandard2nd22gyy,kmul(J21L,J21L),kmadd(PDstandard2nd33gyy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy)),kmul(J21L,kmul(J31L,PDstandard2nd23gyy))),ToReal(2.)))))))); + kmadd(dJ111L,PDstandard2nd1gyy,kmadd(dJ211L,PDstandard2nd2gyy,kmadd(dJ311L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,kmul(J11L,J11L),kmadd(PDstandard2nd22gyy,kmul(J21L,J21L),kmadd(PDstandard2nd33gyy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyy,kmul(J31L,PDstandard2nd13gyy)),kmul(J21L,kmul(J31L,PDstandard2nd23gyy))),ToReal(2)))))))); JacPDstandard2nd11gyz = - kmadd(dJ111L,PDstandard2nd1gyz,kmadd(dJ211L,PDstandard2nd2gyz,kmadd(dJ311L,PDstandard2nd3gyz,kmadd(PDstandard2nd11gyz,kmul(J11L,J11L),kmadd(PDstandard2nd22gyz,kmul(J21L,J21L),kmadd(PDstandard2nd33gyz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz)),kmul(J21L,kmul(J31L,PDstandard2nd23gyz))),ToReal(2.)))))))); + kmadd(dJ111L,PDstandard2nd1gyz,kmadd(dJ211L,PDstandard2nd2gyz,kmadd(dJ311L,PDstandard2nd3gyz,kmadd(PDstandard2nd11gyz,kmul(J11L,J11L),kmadd(PDstandard2nd22gyz,kmul(J21L,J21L),kmadd(PDstandard2nd33gyz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gyz,kmul(J31L,PDstandard2nd13gyz)),kmul(J21L,kmul(J31L,PDstandard2nd23gyz))),ToReal(2)))))))); JacPDstandard2nd11gzz = - kmadd(dJ111L,PDstandard2nd1gzz,kmadd(dJ211L,PDstandard2nd2gzz,kmadd(dJ311L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,kmul(J11L,J11L),kmadd(PDstandard2nd22gzz,kmul(J21L,J21L),kmadd(PDstandard2nd33gzz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz)),kmul(J21L,kmul(J31L,PDstandard2nd23gzz))),ToReal(2.)))))))); + kmadd(dJ111L,PDstandard2nd1gzz,kmadd(dJ211L,PDstandard2nd2gzz,kmadd(dJ311L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,kmul(J11L,J11L),kmadd(PDstandard2nd22gzz,kmul(J21L,J21L),kmadd(PDstandard2nd33gzz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard2nd12gzz,kmul(J31L,PDstandard2nd13gzz)),kmul(J21L,kmul(J31L,PDstandard2nd23gzz))),ToReal(2)))))))); JacPDstandard2nd22gxx = - kmadd(dJ122L,PDstandard2nd1gxx,kmadd(dJ222L,PDstandard2nd2gxx,kmadd(dJ322L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,kmul(J12L,J12L),kmadd(PDstandard2nd22gxx,kmul(J22L,J22L),kmadd(PDstandard2nd33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx)),kmul(J22L,kmul(J32L,PDstandard2nd23gxx))),ToReal(2.)))))))); + kmadd(dJ122L,PDstandard2nd1gxx,kmadd(dJ222L,PDstandard2nd2gxx,kmadd(dJ322L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,kmul(J12L,J12L),kmadd(PDstandard2nd22gxx,kmul(J22L,J22L),kmadd(PDstandard2nd33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxx,kmul(J32L,PDstandard2nd13gxx)),kmul(J22L,kmul(J32L,PDstandard2nd23gxx))),ToReal(2)))))))); JacPDstandard2nd22gxz = - kmadd(dJ122L,PDstandard2nd1gxz,kmadd(dJ222L,PDstandard2nd2gxz,kmadd(dJ322L,PDstandard2nd3gxz,kmadd(PDstandard2nd11gxz,kmul(J12L,J12L),kmadd(PDstandard2nd22gxz,kmul(J22L,J22L),kmadd(PDstandard2nd33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz)),kmul(J22L,kmul(J32L,PDstandard2nd23gxz))),ToReal(2.)))))))); + kmadd(dJ122L,PDstandard2nd1gxz,kmadd(dJ222L,PDstandard2nd2gxz,kmadd(dJ322L,PDstandard2nd3gxz,kmadd(PDstandard2nd11gxz,kmul(J12L,J12L),kmadd(PDstandard2nd22gxz,kmul(J22L,J22L),kmadd(PDstandard2nd33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gxz,kmul(J32L,PDstandard2nd13gxz)),kmul(J22L,kmul(J32L,PDstandard2nd23gxz))),ToReal(2)))))))); JacPDstandard2nd22gzz = - kmadd(dJ122L,PDstandard2nd1gzz,kmadd(dJ222L,PDstandard2nd2gzz,kmadd(dJ322L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,kmul(J12L,J12L),kmadd(PDstandard2nd22gzz,kmul(J22L,J22L),kmadd(PDstandard2nd33gzz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz)),kmul(J22L,kmul(J32L,PDstandard2nd23gzz))),ToReal(2.)))))))); + kmadd(dJ122L,PDstandard2nd1gzz,kmadd(dJ222L,PDstandard2nd2gzz,kmadd(dJ322L,PDstandard2nd3gzz,kmadd(PDstandard2nd11gzz,kmul(J12L,J12L),kmadd(PDstandard2nd22gzz,kmul(J22L,J22L),kmadd(PDstandard2nd33gzz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard2nd12gzz,kmul(J32L,PDstandard2nd13gzz)),kmul(J22L,kmul(J32L,PDstandard2nd23gzz))),ToReal(2)))))))); JacPDstandard2nd33gxx = - kmadd(dJ133L,PDstandard2nd1gxx,kmadd(dJ233L,PDstandard2nd2gxx,kmadd(dJ333L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,kmul(J13L,J13L),kmadd(PDstandard2nd22gxx,kmul(J23L,J23L),kmadd(PDstandard2nd33gxx,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmul(J23L,kmul(J33L,PDstandard2nd23gxx))),ToReal(2.)))))))); + kmadd(dJ133L,PDstandard2nd1gxx,kmadd(dJ233L,PDstandard2nd2gxx,kmadd(dJ333L,PDstandard2nd3gxx,kmadd(PDstandard2nd11gxx,kmul(J13L,J13L),kmadd(PDstandard2nd22gxx,kmul(J23L,J23L),kmadd(PDstandard2nd33gxx,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxx,kmul(J33L,PDstandard2nd13gxx)),kmul(J23L,kmul(J33L,PDstandard2nd23gxx))),ToReal(2)))))))); JacPDstandard2nd33gxy = - kmadd(dJ133L,PDstandard2nd1gxy,kmadd(dJ233L,PDstandard2nd2gxy,kmadd(dJ333L,PDstandard2nd3gxy,kmadd(PDstandard2nd11gxy,kmul(J13L,J13L),kmadd(PDstandard2nd22gxy,kmul(J23L,J23L),kmadd(PDstandard2nd33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmul(J23L,kmul(J33L,PDstandard2nd23gxy))),ToReal(2.)))))))); + kmadd(dJ133L,PDstandard2nd1gxy,kmadd(dJ233L,PDstandard2nd2gxy,kmadd(dJ333L,PDstandard2nd3gxy,kmadd(PDstandard2nd11gxy,kmul(J13L,J13L),kmadd(PDstandard2nd22gxy,kmul(J23L,J23L),kmadd(PDstandard2nd33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gxy,kmul(J33L,PDstandard2nd13gxy)),kmul(J23L,kmul(J33L,PDstandard2nd23gxy))),ToReal(2)))))))); JacPDstandard2nd33gyy = - kmadd(dJ133L,PDstandard2nd1gyy,kmadd(dJ233L,PDstandard2nd2gyy,kmadd(dJ333L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,kmul(J13L,J13L),kmadd(PDstandard2nd22gyy,kmul(J23L,J23L),kmadd(PDstandard2nd33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmul(J23L,kmul(J33L,PDstandard2nd23gyy))),ToReal(2.)))))))); + kmadd(dJ133L,PDstandard2nd1gyy,kmadd(dJ233L,PDstandard2nd2gyy,kmadd(dJ333L,PDstandard2nd3gyy,kmadd(PDstandard2nd11gyy,kmul(J13L,J13L),kmadd(PDstandard2nd22gyy,kmul(J23L,J23L),kmadd(PDstandard2nd33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard2nd12gyy,kmul(J33L,PDstandard2nd13gyy)),kmul(J23L,kmul(J33L,PDstandard2nd23gyy))),ToReal(2)))))))); JacPDstandard2nd12gxy = kmadd(J12L,kmadd(J11L,PDstandard2nd11gxy,kmadd(J21L,PDstandard2nd12gxy,kmul(J31L,PDstandard2nd13gxy))),kmadd(J11L,kmadd(J22L,PDstandard2nd12gxy,kmul(J32L,PDstandard2nd13gxy)),kmadd(dJ112L,PDstandard2nd1gxy,kmadd(J22L,kmadd(J21L,PDstandard2nd22gxy,kmul(J31L,PDstandard2nd23gxy)),kmadd(dJ212L,PDstandard2nd2gxy,kmadd(J32L,kmadd(J21L,PDstandard2nd23gxy,kmul(J31L,PDstandard2nd33gxy)),kmul(dJ312L,PDstandard2nd3gxy))))))); @@ -955,105 +996,102 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int } CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED detg = - 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.)))))); + knmsub(gyyL,kmul(gxzL,gxzL),knmsub(gxxL,kmul(gyzL,gyzL),kmadd(gzzL,kmsub(gxxL,gyyL,kmul(gxyL,gxyL)),kmul(gxyL,kmul(gxzL,kmul(gyzL,ToReal(2))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED invdetg = kdiv(ToReal(1),detg); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv11 = - kmul(invdetg,kmadd(gyyL,gzzL,kmul(kmul(gyzL,gyzL),ToReal(-1.)))); + kmul(invdetg,kmsub(gyyL,gzzL,kmul(gyzL,gyzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv12 = - kmul(invdetg,kmadd(gxzL,gyzL,kmul(gxyL,kmul(gzzL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv13 = - kmul(invdetg,kmadd(gxyL,gyzL,kmul(gxzL,kmul(gyyL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv21 = - kmul(invdetg,kmadd(gxzL,gyzL,kmul(gxyL,kmul(gzzL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxzL,gyzL,kmul(gxyL,gzzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv22 = - kmul(invdetg,kmadd(gxxL,gzzL,kmul(kmul(gxzL,gxzL),ToReal(-1.)))); + kmul(invdetg,kmsub(gxxL,gzzL,kmul(gxzL,gxzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv23 = - kmul(invdetg,kmadd(gxyL,gxzL,kmul(gxxL,kmul(gyzL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv31 = - kmul(invdetg,kmadd(gxyL,gyzL,kmul(gxzL,kmul(gyyL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxyL,gyzL,kmul(gxzL,gyyL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv32 = - kmul(invdetg,kmadd(gxyL,gxzL,kmul(gxxL,kmul(gyzL,ToReal(-1.))))); + kmul(invdetg,kmsub(gxyL,gxzL,kmul(gxxL,gyzL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gInv33 = - kmul(invdetg,kmadd(gxxL,gyyL,kmul(kmul(gxyL,gxyL),ToReal(-1.)))); + kmul(invdetg,kmsub(gxxL,gyyL,kmul(gxyL,gxyL))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma111 = - kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard2nd1gxx,kmadd(kmadd(gInv12,JacPDstandard2nd2gxx,kmul(gInv13,JacPDstandard2nd3gxx)),ToReal(-1.),kmul(kmadd(gInv12,JacPDstandard2nd1gxy,kmul(gInv13,JacPDstandard2nd1gxz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv11,JacPDstandard2nd1gxx,knmsub(gInv12,JacPDstandard2nd2gxx,kmsub(kmadd(gInv12,JacPDstandard2nd1gxy,kmul(gInv13,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv13,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma211 = - kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard2nd1gxx,kmadd(kmadd(gInv22,JacPDstandard2nd2gxx,kmul(gInv23,JacPDstandard2nd3gxx)),ToReal(-1.),kmul(kmadd(gInv22,JacPDstandard2nd1gxy,kmul(gInv23,JacPDstandard2nd1gxz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv21,JacPDstandard2nd1gxx,knmsub(gInv22,JacPDstandard2nd2gxx,kmsub(kmadd(gInv22,JacPDstandard2nd1gxy,kmul(gInv23,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv23,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma311 = - kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard2nd1gxx,kmadd(kmadd(gInv32,JacPDstandard2nd2gxx,kmul(gInv33,JacPDstandard2nd3gxx)),ToReal(-1.),kmul(kmadd(gInv32,JacPDstandard2nd1gxy,kmul(gInv33,JacPDstandard2nd1gxz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv31,JacPDstandard2nd1gxx,knmsub(gInv32,JacPDstandard2nd2gxx,kmsub(kmadd(gInv32,JacPDstandard2nd1gxy,kmul(gInv33,JacPDstandard2nd1gxz)),ToReal(2),kmul(gInv33,JacPDstandard2nd3gxx))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma121 = - kmul(kmadd(gInv12,JacPDstandard2nd1gyy,kmadd(gInv11,JacPDstandard2nd2gxx,kmul(gInv13,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd3gxy,ToReal(-1.),JacPDstandard2nd2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv12,JacPDstandard2nd1gyy,kmadd(gInv11,JacPDstandard2nd2gxx,kmul(gInv13,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma221 = - kmul(kmadd(gInv22,JacPDstandard2nd1gyy,kmadd(gInv21,JacPDstandard2nd2gxx,kmul(gInv23,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd3gxy,ToReal(-1.),JacPDstandard2nd2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv22,JacPDstandard2nd1gyy,kmadd(gInv21,JacPDstandard2nd2gxx,kmul(gInv23,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma321 = - kmul(kmadd(gInv32,JacPDstandard2nd1gyy,kmadd(gInv31,JacPDstandard2nd2gxx,kmul(gInv33,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd3gxy,ToReal(-1.),JacPDstandard2nd2gxz))))),ToReal(0.5)); + kmul(kmadd(gInv32,JacPDstandard2nd1gyy,kmadd(gInv31,JacPDstandard2nd2gxx,kmul(gInv33,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd2gxz,JacPDstandard2nd3gxy))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma131 = - kmul(kmadd(gInv13,JacPDstandard2nd1gzz,kmadd(gInv11,JacPDstandard2nd3gxx,kmul(gInv12,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd2gxz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv13,JacPDstandard2nd1gzz,kmadd(gInv11,JacPDstandard2nd3gxx,kmul(gInv12,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma231 = - kmul(kmadd(gInv23,JacPDstandard2nd1gzz,kmadd(gInv21,JacPDstandard2nd3gxx,kmul(gInv22,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd2gxz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv23,JacPDstandard2nd1gzz,kmadd(gInv21,JacPDstandard2nd3gxx,kmul(gInv22,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma331 = - kmul(kmadd(gInv33,JacPDstandard2nd1gzz,kmadd(gInv31,JacPDstandard2nd3gxx,kmul(gInv32,kadd(JacPDstandard2nd1gyz,kmadd(JacPDstandard2nd2gxz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv33,JacPDstandard2nd1gzz,kmadd(gInv31,JacPDstandard2nd3gxx,kmul(gInv32,kadd(JacPDstandard2nd1gyz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd2gxz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma122 = - kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2nd2gyy,kmadd(gInv11,kmadd(JacPDstandard2nd1gyy,ToReal(-1.),kmul(JacPDstandard2nd2gxy,ToReal(2.))),kmul(gInv13,kmadd(JacPDstandard2nd3gyy,ToReal(-1.),kmul(JacPDstandard2nd2gyz,ToReal(2.))))))); + kmul(ToReal(0.5),kmadd(gInv12,JacPDstandard2nd2gyy,kmadd(gInv11,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv13,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma222 = - kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2nd2gyy,kmadd(gInv21,kmadd(JacPDstandard2nd1gyy,ToReal(-1.),kmul(JacPDstandard2nd2gxy,ToReal(2.))),kmul(gInv23,kmadd(JacPDstandard2nd3gyy,ToReal(-1.),kmul(JacPDstandard2nd2gyz,ToReal(2.))))))); + kmul(ToReal(0.5),kmadd(gInv22,JacPDstandard2nd2gyy,kmadd(gInv21,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv23,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma322 = - kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2nd2gyy,kmadd(gInv31,kmadd(JacPDstandard2nd1gyy,ToReal(-1.),kmul(JacPDstandard2nd2gxy,ToReal(2.))),kmul(gInv33,kmadd(JacPDstandard2nd3gyy,ToReal(-1.),kmul(JacPDstandard2nd2gyz,ToReal(2.))))))); + kmul(ToReal(0.5),kmadd(gInv32,JacPDstandard2nd2gyy,kmadd(gInv31,kmsub(JacPDstandard2nd2gxy,ToReal(2),JacPDstandard2nd1gyy),kmul(gInv33,kmsub(JacPDstandard2nd2gyz,ToReal(2),JacPDstandard2nd3gyy))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma132 = - kmul(kmadd(gInv13,JacPDstandard2nd2gzz,kmadd(gInv12,JacPDstandard2nd3gyy,kmul(gInv11,kadd(JacPDstandard2nd2gxz,kmadd(JacPDstandard2nd1gyz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv13,JacPDstandard2nd2gzz,kmadd(gInv12,JacPDstandard2nd3gyy,kmul(gInv11,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma232 = - kmul(kmadd(gInv23,JacPDstandard2nd2gzz,kmadd(gInv22,JacPDstandard2nd3gyy,kmul(gInv21,kadd(JacPDstandard2nd2gxz,kmadd(JacPDstandard2nd1gyz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv23,JacPDstandard2nd2gzz,kmadd(gInv22,JacPDstandard2nd3gyy,kmul(gInv21,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma332 = - kmul(kmadd(gInv33,JacPDstandard2nd2gzz,kmadd(gInv32,JacPDstandard2nd3gyy,kmul(gInv31,kadd(JacPDstandard2nd2gxz,kmadd(JacPDstandard2nd1gyz,ToReal(-1.),JacPDstandard2nd3gxy))))),ToReal(0.5)); + kmul(kmadd(gInv33,JacPDstandard2nd2gzz,kmadd(gInv32,JacPDstandard2nd3gyy,kmul(gInv31,kadd(JacPDstandard2nd2gxz,ksub(JacPDstandard2nd3gxy,JacPDstandard2nd1gyz))))),ToReal(0.5)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma133 = - kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard2nd3gzz,kmadd(kmadd(gInv11,JacPDstandard2nd1gzz,kmul(gInv12,JacPDstandard2nd2gzz)),ToReal(-1.),kmul(kmadd(gInv11,JacPDstandard2nd3gxz,kmul(gInv12,JacPDstandard2nd3gyz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv13,JacPDstandard2nd3gzz,kmadd(gInv11,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv12,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma233 = - kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard2nd3gzz,kmadd(kmadd(gInv21,JacPDstandard2nd1gzz,kmul(gInv22,JacPDstandard2nd2gzz)),ToReal(-1.),kmul(kmadd(gInv21,JacPDstandard2nd3gxz,kmul(gInv22,JacPDstandard2nd3gyz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv23,JacPDstandard2nd3gzz,kmadd(gInv21,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv22,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gamma333 = - kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard2nd3gzz,kmadd(kmadd(gInv31,JacPDstandard2nd1gzz,kmul(gInv32,JacPDstandard2nd2gzz)),ToReal(-1.),kmul(kmadd(gInv31,JacPDstandard2nd3gxz,kmul(gInv32,JacPDstandard2nd3gyz)),ToReal(2.))))); + kmul(ToReal(0.5),kmadd(gInv33,JacPDstandard2nd3gzz,kmadd(gInv31,kmsub(JacPDstandard2nd3gxz,ToReal(2),JacPDstandard2nd1gzz),kmul(gInv32,kmsub(JacPDstandard2nd3gyz,ToReal(2),JacPDstandard2nd2gzz))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED xmoved = - kadd(xL,ToReal(-1.*xorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED xmoved = ksub(xL,ToReal(xorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ymoved = - kadd(yL,ToReal(-1.*yorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED ymoved = ksub(yL,ToReal(yorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED zmoved = - kadd(zL,ToReal(-1.*zorig)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED zmoved = ksub(zL,ToReal(zorig)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va1 = kmul(ymoved,ToReal(-1.)); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED va1 = kneg(ymoved); 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)); @@ -1062,13 +1100,13 @@ static void WeylScal4_psi4_calc_2nd_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,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.)))))))); + kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv11,va2,kmul(gInv12,va1)),kmadd(vb1,kmsub(gInv12,va3,kmul(gInv13,va2)),kmul(vb2,kmsub(gInv13,va1,kmul(gInv11,va3)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED vc2 = - 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.)))))))); + kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv21,va2,kmul(gInv22,va1)),kmadd(vb1,kmsub(gInv22,va3,kmul(gInv23,va2)),kmul(vb2,kmsub(gInv23,va1,kmul(gInv21,va3)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED vc3 = - 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.)))))))); + kmul(ksqrt(detg),kmadd(vb3,kmsub(gInv31,va2,kmul(gInv32,va1)),kmadd(vb1,kmsub(gInv32,va3,kmul(gInv33,va2)),kmul(vb2,kmsub(gInv33,va1,kmul(gInv31,va3)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wa1 = va1; @@ -1077,7 +1115,7 @@ static void WeylScal4_psi4_calc_2nd_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)); @@ -1088,17 +1126,14 @@ static void WeylScal4_psi4_calc_2nd_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 = - kmadd(ea1,kmul(omega12,ToReal(-1.)),vb1); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb1 = knmsub(ea1,omega12,vb1); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb2 = - kmadd(ea2,kmul(omega12,ToReal(-1.)),vb2); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb2 = knmsub(ea2,omega12,vb2); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb3 = - kmadd(ea3,kmul(omega12,ToReal(-1.)),vb3); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wb3 = knmsub(ea3,omega12,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)); @@ -1113,16 +1148,16 @@ static void WeylScal4_psi4_calc_2nd_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 = - kmadd(kmadd(ea1,omega13,kmul(eb1,omega23)),ToReal(-1.),vc1); + ksub(vc1,kmadd(eb1,omega23,kmul(ea1,omega13))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wc2 = - kmadd(kmadd(ea2,omega13,kmul(eb2,omega23)),ToReal(-1.),vc2); + ksub(vc2,kmadd(eb2,omega23,kmul(ea2,omega13))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED wc3 = - kmadd(kmadd(ea3,omega13,kmul(eb3,omega23)),ToReal(-1.),vc3); + ksub(vc3,kmadd(eb3,omega23,kmul(ea3,omega13))); 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)); @@ -1133,14 +1168,11 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED isqrt2 = ToReal(0.707106781186547524); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n1 = - kmul(eb1,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n1 = kneg(kmul(eb1,isqrt2)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n2 = - kmul(eb2,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n2 = kneg(kmul(eb2,isqrt2)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n3 = - kmul(eb3,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED n3 = kneg(kmul(eb3,isqrt2)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rm1 = kmul(ec1,isqrt2); @@ -1160,160 +1192,157 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rmbar3 = kmul(ec3,isqrt2); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar1 = - kmul(ea1,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar1 = kneg(kmul(ea1,isqrt2)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar2 = - kmul(ea2,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar2 = kneg(kmul(ea2,isqrt2)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar3 = - kmul(ea3,kmul(isqrt2,ToReal(-1.))); + CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED imbar3 = kneg(kmul(ea3,isqrt2)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED nn = isqrt2; CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1212 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxy,kadd(JacPDstandard2nd21gxy,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(JacPDstandard2nd11gyy,JacPDstandard2nd22gxx),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxy,kadd(JacPDstandard2nd21gxy,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),ksub(kmsub(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),JacPDstandard2nd22gxx),JacPDstandard2nd11gyy))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1213 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxz,kadd(JacPDstandard2nd31gxy,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(JacPDstandard2nd11gyz,JacPDstandard2nd23gxx),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd12gxz,kadd(JacPDstandard2nd31gxy,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),ksub(kmsub(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),JacPDstandard2nd23gxx),JacPDstandard2nd11gyz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1223 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd22gxz,kadd(JacPDstandard2nd31gyy,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(JacPDstandard2nd12gyz,JacPDstandard2nd23gxy),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd22gxz,kadd(JacPDstandard2nd31gyy,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),ksub(kmsub(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),JacPDstandard2nd23gxy),JacPDstandard2nd12gyz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1313 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxz,kadd(JacPDstandard2nd31gxz,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(JacPDstandard2nd11gzz,JacPDstandard2nd33gxx),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd13gxz,kadd(JacPDstandard2nd31gxz,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),ksub(kmsub(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),JacPDstandard2nd33gxx),JacPDstandard2nd11gzz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R1323 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd31gyz,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(JacPDstandard2nd12gzz,JacPDstandard2nd33gxy),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd23gxz,kadd(JacPDstandard2nd31gyz,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),ksub(kmsub(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),JacPDstandard2nd33gxy),JacPDstandard2nd12gzz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R2323 = - kmul(ToReal(0.5),kadd(JacPDstandard2nd23gyz,kadd(JacPDstandard2nd32gyz,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(JacPDstandard2nd22gzz,JacPDstandard2nd33gyy),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.))))))); + kmul(ToReal(0.5),kadd(JacPDstandard2nd23gyz,kadd(JacPDstandard2nd32gyz,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),ksub(kmsub(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),JacPDstandard2nd33gyy),JacPDstandard2nd22gzz))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1212 = - kmadd(kxxL,kyyL,kmadd(kmul(kxyL,kxyL),ToReal(-1.),R1212)); + kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1213 = - kmadd(kxxL,kyzL,kmadd(kxyL,kmul(kxzL,ToReal(-1.)),R1213)); + kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1223 = - kmadd(kxyL,kyzL,kmadd(kxzL,kmul(kyyL,ToReal(-1.)),R1223)); + kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1223)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1313 = - kmadd(kxxL,kzzL,kmadd(kmul(kxzL,kxzL),ToReal(-1.),R1313)); + kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p1323 = - kmadd(kxyL,kzzL,kmadd(kxzL,kmul(kyzL,ToReal(-1.)),R1323)); + kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R4p2323 = - kmadd(kyyL,kzzL,kmadd(kmul(kyzL,kyzL),ToReal(-1.),R2323)); + kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,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(kxyL,gamma221,kmadd(kxzL,gamma321,kmadd(kmadd(kxyL,gamma111,kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,JacPDstandard2nd2kxx))),ToReal(-1.),JacPDstandard2nd1kxy)))); + kmadd(kxxL,gamma121,kmadd(kxzL,gamma321,kadd(JacPDstandard2nd1kxy,knmsub(kyyL,gamma211,knmsub(kyzL,gamma311,kmsub(kxyL,ksub(gamma221,gamma111),JacPDstandard2nd2kxx)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro113 = - kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,kmadd(kxzL,gamma331,kmadd(kmadd(kxzL,gamma111,kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,JacPDstandard2nd3kxx))),ToReal(-1.),JacPDstandard2nd1kxz)))); + kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,kadd(JacPDstandard2nd1kxz,knmsub(kyzL,gamma211,knmsub(kzzL,gamma311,kmsub(kxzL,ksub(gamma331,gamma111),JacPDstandard2nd3kxx)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro121 = - kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,kadd(JacPDstandard2nd2kxx,kmadd(kmadd(kxxL,gamma121,kmadd(kxzL,gamma321,JacPDstandard2nd1kxy)),ToReal(-1.),kmul(kxyL,kmadd(gamma221,ToReal(-1.),gamma111)))))); + kmadd(kyyL,gamma211,kmadd(kyzL,gamma311,kadd(JacPDstandard2nd2kxx,knmsub(kxxL,gamma121,knmsub(kxzL,gamma321,kmsub(kxyL,ksub(gamma111,gamma221),JacPDstandard2nd1kxy)))))); - 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,kmadd(kyzL,gamma331,kmadd(kmadd(kxzL,gamma121,kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,JacPDstandard2nd3kxy))),ToReal(-1.),JacPDstandard2nd2kxz)))); + kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kadd(JacPDstandard2nd2kxz,knmsub(kxzL,gamma121,knmsub(kzzL,gamma321,kmsub(kyzL,ksub(gamma331,gamma221),JacPDstandard2nd3kxy)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro131 = - kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,kadd(JacPDstandard2nd3kxx,kmadd(kmadd(kxxL,gamma131,kmadd(kxyL,gamma231,JacPDstandard2nd1kxz)),ToReal(-1.),kmul(kxzL,kmadd(gamma331,ToReal(-1.),gamma111)))))); + kmadd(kyzL,gamma211,kmadd(kzzL,gamma311,kadd(JacPDstandard2nd3kxx,knmsub(kxxL,gamma131,knmsub(kxyL,gamma231,kmsub(kxzL,ksub(gamma111,gamma331),JacPDstandard2nd1kxz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro132 = - kmadd(kxzL,gamma121,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,kmadd(kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,JacPDstandard2nd2kxz)),ToReal(-1.),kmul(kyzL,kmadd(gamma331,ToReal(-1.),gamma221)))))); + kmadd(kxzL,gamma121,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,knmsub(kxyL,gamma131,knmsub(kyyL,gamma231,kmsub(kyzL,ksub(gamma221,gamma331),JacPDstandard2nd2kxz)))))); - 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(kxyL,gamma222,kmadd(kxzL,gamma322,kmadd(kmadd(kxyL,gamma121,kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,JacPDstandard2nd2kxy))),ToReal(-1.),JacPDstandard2nd1kyy)))); + kmadd(kxxL,gamma122,kmadd(kxzL,gamma322,kadd(JacPDstandard2nd1kyy,knmsub(kyyL,gamma221,knmsub(kyzL,gamma321,kmsub(kxyL,ksub(gamma222,gamma121),JacPDstandard2nd2kxy)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro213 = - kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kmadd(kxzL,gamma332,kmadd(kmadd(kxzL,gamma121,kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,JacPDstandard2nd3kxy))),ToReal(-1.),JacPDstandard2nd1kyz)))); + kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,kadd(JacPDstandard2nd1kyz,knmsub(kyzL,gamma221,knmsub(kzzL,gamma321,kmsub(kxzL,ksub(gamma332,gamma121),JacPDstandard2nd3kxy)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro221 = - kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,kadd(JacPDstandard2nd2kxy,kmadd(kmadd(kxxL,gamma122,kmadd(kxzL,gamma322,JacPDstandard2nd1kyy)),ToReal(-1.),kmul(kxyL,kmadd(gamma222,ToReal(-1.),gamma121)))))); + kmadd(kyyL,gamma221,kmadd(kyzL,gamma321,kadd(JacPDstandard2nd2kxy,knmsub(kxxL,gamma122,knmsub(kxzL,gamma322,kmsub(kxyL,ksub(gamma121,gamma222),JacPDstandard2nd1kyy)))))); - 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,kmadd(kyzL,gamma332,kmadd(kmadd(kxzL,gamma122,kmadd(kyzL,gamma222,kmadd(kzzL,gamma322,JacPDstandard2nd3kyy))),ToReal(-1.),JacPDstandard2nd2kyz)))); + kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,kadd(JacPDstandard2nd2kyz,knmsub(kxzL,gamma122,knmsub(kzzL,gamma322,kmsub(kyzL,ksub(gamma332,gamma222),JacPDstandard2nd3kyy)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro231 = - kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,kmadd(kmadd(kxxL,gamma132,kmadd(kxyL,gamma232,JacPDstandard2nd1kyz)),ToReal(-1.),kmul(kxzL,kmadd(gamma332,ToReal(-1.),gamma121)))))); + kmadd(kyzL,gamma221,kmadd(kzzL,gamma321,kadd(JacPDstandard2nd3kxy,knmsub(kxxL,gamma132,knmsub(kxyL,gamma232,kmsub(kxzL,ksub(gamma121,gamma332),JacPDstandard2nd1kyz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro232 = - kmadd(kxzL,gamma122,kmadd(kzzL,gamma322,kadd(JacPDstandard2nd3kyy,kmadd(kmadd(kxyL,gamma132,kmadd(kyyL,gamma232,JacPDstandard2nd2kyz)),ToReal(-1.),kmul(kyzL,kmadd(gamma332,ToReal(-1.),gamma222)))))); + kmadd(kxzL,gamma122,kmadd(kzzL,gamma322,kadd(JacPDstandard2nd3kyy,knmsub(kxyL,gamma132,knmsub(kyyL,gamma232,kmsub(kyzL,ksub(gamma222,gamma332),JacPDstandard2nd2kyz)))))); - 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(kxyL,gamma232,kmadd(kxzL,gamma332,kmadd(kmadd(kxyL,gamma131,kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,JacPDstandard2nd2kxz))),ToReal(-1.),JacPDstandard2nd1kyz)))); + kmadd(kxxL,gamma132,kmadd(kxzL,gamma332,kadd(JacPDstandard2nd1kyz,knmsub(kyyL,gamma231,knmsub(kyzL,gamma331,kmsub(kxyL,ksub(gamma232,gamma131),JacPDstandard2nd2kxz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro313 = - kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,kmadd(kxzL,gamma333,kmadd(kmadd(kxzL,gamma131,kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,JacPDstandard2nd3kxz))),ToReal(-1.),JacPDstandard2nd1kzz)))); + kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,kadd(JacPDstandard2nd1kzz,knmsub(kyzL,gamma231,knmsub(kzzL,gamma331,kmsub(kxzL,ksub(gamma333,gamma131),JacPDstandard2nd3kxz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro321 = - kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kadd(JacPDstandard2nd2kxz,kmadd(kmadd(kxxL,gamma132,kmadd(kxzL,gamma332,JacPDstandard2nd1kyz)),ToReal(-1.),kmul(kxyL,kmadd(gamma232,ToReal(-1.),gamma131)))))); + kmadd(kyyL,gamma231,kmadd(kyzL,gamma331,kadd(JacPDstandard2nd2kxz,knmsub(kxxL,gamma132,knmsub(kxzL,gamma332,kmsub(kxyL,ksub(gamma131,gamma232),JacPDstandard2nd1kyz)))))); - 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,kmadd(kyzL,gamma333,kmadd(kmadd(kxzL,gamma132,kmadd(kyzL,gamma232,kmadd(kzzL,gamma332,JacPDstandard2nd3kyz))),ToReal(-1.),JacPDstandard2nd2kzz)))); + kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,kadd(JacPDstandard2nd2kzz,knmsub(kxzL,gamma132,knmsub(kzzL,gamma332,kmsub(kyzL,ksub(gamma333,gamma232),JacPDstandard2nd3kyz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro331 = - kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,kadd(JacPDstandard2nd3kxz,kmadd(kmadd(kxxL,gamma133,kmadd(kxyL,gamma233,JacPDstandard2nd1kzz)),ToReal(-1.),kmul(kxzL,kmadd(gamma333,ToReal(-1.),gamma131)))))); + kmadd(kyzL,gamma231,kmadd(kzzL,gamma331,kadd(JacPDstandard2nd3kxz,knmsub(kxxL,gamma133,knmsub(kxyL,gamma233,kmsub(kxzL,ksub(gamma131,gamma333),JacPDstandard2nd1kzz)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ro332 = - kmadd(kxzL,gamma132,kmadd(kzzL,gamma332,kadd(JacPDstandard2nd3kyz,kmadd(kmadd(kxyL,gamma133,kmadd(kyyL,gamma233,JacPDstandard2nd2kzz)),ToReal(-1.),kmul(kyzL,kmadd(gamma333,ToReal(-1.),gamma232)))))); + kmadd(kxzL,gamma132,kmadd(kzzL,gamma332,kadd(JacPDstandard2nd3kyz,knmsub(kxyL,gamma133,knmsub(kyyL,gamma233,kmsub(kyzL,ksub(gamma232,gamma333),JacPDstandard2nd2kzz)))))); - 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,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))))); + 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))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo12 = - 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)))))))); + 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))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo13 = - 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.)))))))); + 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)))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo21 = - 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)))))))); + 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))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo22 = - 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))))); + 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))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo23 = - 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.)))))))))); + 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))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo31 = - 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.)))))))); + 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)))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo32 = - 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.)))))))))); + 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))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Rojo33 = - 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))))); + 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))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi4rL = - 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.)))))))))))); + knmsub(kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2)))))),kmsub(imbar2,imbar2,kmul(rmbar2,rmbar2)),knmsub(kmul(nn,nn),kmadd(imbar1,kmadd(imbar2,kadd(Rojo12,Rojo21),kmul(imbar3,kadd(Rojo13,Rojo31))),knmsub(rmbar1,kmadd(rmbar2,kadd(Rojo12,Rojo21),kmul(rmbar3,kadd(Rojo13,Rojo31))),kmadd(Rojo23,kmsub(imbar2,imbar3,kmul(rmbar2,rmbar3)),kmadd(Rojo32,kmsub(imbar2,imbar3,kmul(rmbar2,rmbar3)),kmadd(Rojo11,kmsub(imbar1,imbar1,kmul(rmbar1,rmbar1)),kmadd(Rojo22,kmsub(imbar2,imbar2,kmul(rmbar2,rmbar2)),kmul(Rojo33,kmsub(imbar3,imbar3,kmul(rmbar3,rmbar3))))))))),kmsub(kmadd(kmsub(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmadd(n2,R4p1223,kmul(n3,R4p1323)))),kmsub(imbar1,imbar2,kmul(rmbar1,rmbar2)),kmadd(kmadd(n1,kmul(n2,R4p1213),kmadd(n1,kmul(n3,R4p1313),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2))))),kmsub(imbar1,imbar3,kmul(rmbar1,rmbar3)),kmadd(kmsub(rmbar2,rmbar3,kmul(imbar2,imbar3)),kmadd(R4p1213,kmul(n1,n1),kmsub(n1,kmsub(n2,R4p1223,kmul(n3,R4p1323)),kmul(n2,kmul(n3,R4p2323)))),kmul(nn,kmadd(kmadd(n1,Ro112,kmadd(n2,Ro122,kmul(n3,Ro132))),kmsub(rmbar1,rmbar2,kmul(imbar1,imbar2)),kmadd(kmadd(n1,Ro211,kmadd(n2,Ro221,kmul(n3,Ro231))),kmsub(rmbar1,rmbar2,kmul(imbar1,imbar2)),kmadd(kmadd(n1,Ro113,kmadd(n2,Ro123,kmul(n3,Ro133))),kmsub(rmbar1,rmbar3,kmul(imbar1,imbar3)),kmadd(kmadd(n1,Ro311,kmadd(n2,Ro321,kmul(n3,Ro331))),kmsub(rmbar1,rmbar3,kmul(imbar1,imbar3)),kmadd(kmadd(n1,Ro213,kmadd(n2,Ro223,kmul(n3,Ro233))),kmsub(rmbar2,rmbar3,kmul(imbar2,imbar3)),kmadd(kmadd(n1,Ro312,kmadd(n2,Ro322,kmul(n3,Ro332))),kmsub(rmbar2,rmbar3,kmul(imbar2,imbar3)),kmadd(kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131))),kmsub(rmbar1,rmbar1,kmul(imbar1,imbar1)),kmadd(kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232))),kmsub(rmbar2,rmbar2,kmul(imbar2,imbar2)),kmul(kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))),kmsub(rmbar3,rmbar3,kmul(imbar3,imbar3))))))))))))))),ToReal(2),kmadd(kmsub(imbar3,imbar3,kmul(rmbar3,rmbar3)),kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2)))))),kmul(kmsub(imbar1,imbar1,kmul(rmbar1,rmbar1)),kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Psi4iL = - 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.)))))))))))); + 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)),kmsub(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmadd(n2,R4p1223,kmul(n3,R4p1323)))),kmul(nn,kmadd(kmadd(im1,kmul(rm1,kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131)))),kmul(im2,kmul(rm2,kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232)))))),ToReal(-2),kmsub(im3,kmul(rm3,kmul(kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))),ToReal(-2))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kadd(Ro211,Ro112),kmadd(n3,kadd(Ro231,Ro132),kmul(n2,kadd(Ro221,Ro122)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(n1,kadd(Ro312,Ro213),kmadd(n3,kadd(Ro332,Ro233),kmul(n2,kadd(Ro322,Ro223)))),kmul(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kadd(Ro311,Ro113),kmadd(n3,kadd(Ro331,Ro133),kmul(n2,kadd(Ro321,Ro123)))))))))))),ToReal(2),kmsub(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),kmsub(n1,kmsub(n2,R4p1223,kmul(n3,R4p1323)),kmul(n2,kmul(n3,R4p2323)))),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),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,kimin,kimax); + vec_store_partial_prepare(i,vecimin,vecimax); vec_store_nta_partial(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } @@ -1336,7 +1365,7 @@ extern "C" void WeylScal4_psi4_calc_2nd(CCTK_ARGUMENTS) return; } - const char *const groups[] = { + const char* const groups[] = { "admbase::curv", "admbase::metric", "grid::coordinates", @@ -1344,23 +1373,33 @@ extern "C" void WeylScal4_psi4_calc_2nd(CCTK_ARGUMENTS) "WeylScal4::Psi4r_group"}; GenericFD_AssertGroupStorage(cctkGH, "WeylScal4_psi4_calc_2nd", 5, groups); - switch(fdOrder) + switch (fdOrder) { case 2: + { GenericFD_EnsureStencilFits(cctkGH, "WeylScal4_psi4_calc_2nd", 1, 1, 1); break; + } case 4: + { GenericFD_EnsureStencilFits(cctkGH, "WeylScal4_psi4_calc_2nd", 1, 1, 1); break; + } case 6: + { GenericFD_EnsureStencilFits(cctkGH, "WeylScal4_psi4_calc_2nd", 1, 1, 1); break; + } case 8: + { GenericFD_EnsureStencilFits(cctkGH, "WeylScal4_psi4_calc_2nd", 1, 1, 1); break; + } + default: + CCTK_BUILTIN_UNREACHABLE(); } GenericFD_LoopOverInterior(cctkGH, WeylScal4_psi4_calc_2nd_Body); |