/* File produced by Kranc */ #define KRANC_C #include #include #include #include #include #include "cctk.h" #include "cctk_Arguments.h" #include "cctk_Parameters.h" #include "GenericFD.h" #include "Differencing.h" #include "cctk_Loop.h" #include "loopcontrol.h" #include "vectors.h" /* Define macros used in calculations */ #define INITVALUE (42) #define ScalarINV(x) ((CCTK_REAL)1.0 / (x)) #define ScalarSQR(x) ((x) * (x)) #define ScalarCUB(x) ((x) * ScalarSQR(x)) #define ScalarQAD(x) (ScalarSQR(ScalarSQR(x))) #define INV(x) (kdiv(ToReal(1.0),x)) #define SQR(x) (kmul(x,x)) #define CUB(x) (kmul(x,SQR(x))) #define QAD(x) (SQR(SQR(x))) extern "C" void ML_CCZ4_Dissipation_SelectBCs(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; if (cctk_iteration % ML_CCZ4_Dissipation_calc_every != ML_CCZ4_Dissipation_calc_offset) return; CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_curvrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_curvrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_dtlapserhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_dtlapserhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_dtshiftrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_dtshiftrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_Gammarhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_Gammarhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_lapserhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_lapserhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_log_confacrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_log_confacrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_metricrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_metricrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_shiftrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_shiftrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_Thetarhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_Thetarhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_CCZ4::ML_trace_curvrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_CCZ4::ML_trace_curvrhs."); return; } static void ML_CCZ4_Dissipation_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; /* Include user-supplied include files */ /* Initialise finite differencing variables */ 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 */ const CCTK_REAL_VEC p1o1024dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0009765625),dx); const CCTK_REAL_VEC p1o1024dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0009765625),dy); const CCTK_REAL_VEC p1o1024dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0009765625),dz); const CCTK_REAL_VEC p1o120dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00833333333333333333333333333333),dx); const CCTK_REAL_VEC p1o120dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00833333333333333333333333333333),dy); const CCTK_REAL_VEC p1o120dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00833333333333333333333333333333),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 p1o1680dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000595238095238095238095238095238),dx); const CCTK_REAL_VEC p1o1680dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000595238095238095238095238095238),dy); const CCTK_REAL_VEC p1o1680dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.000595238095238095238095238095238),dz); 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 p1o24dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0416666666666666666666666666667),dx); const CCTK_REAL_VEC p1o24dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0416666666666666666666666666667),dy); const CCTK_REAL_VEC p1o24dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0416666666666666666666666666667),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 p1o4dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),dx); 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 p1o4dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),dy); const CCTK_REAL_VEC p1o4dydz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),kmul(dz,dy)); const CCTK_REAL_VEC p1o4dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.25),dz); 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 p1o560dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00178571428571428571428571428571),dx); const CCTK_REAL_VEC p1o560dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00178571428571428571428571428571),dy); const CCTK_REAL_VEC p1o560dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.00178571428571428571428571428571),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 p1o64dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.015625),dx); const CCTK_REAL_VEC p1o64dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.015625),dy); const CCTK_REAL_VEC p1o64dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.015625),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 p1odx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),dx); const CCTK_REAL_VEC p1odx2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dx,dx)); const CCTK_REAL_VEC p1ody CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),dy); const CCTK_REAL_VEC p1ody2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dy,dy)); const CCTK_REAL_VEC p1odz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),dz); const CCTK_REAL_VEC p1odz2 CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(1),kmul(dz,dz)); const CCTK_REAL_VEC pm1o120dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00833333333333333333333333333333),dx); const CCTK_REAL_VEC pm1o120dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00833333333333333333333333333333),dy); const CCTK_REAL_VEC pm1o120dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00833333333333333333333333333333),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)); const CCTK_REAL_VEC pm1o16dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0625),dx); const CCTK_REAL_VEC pm1o16dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0625),dy); const CCTK_REAL_VEC pm1o16dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0625),dz); const CCTK_REAL_VEC pm1o256dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00390625),dx); const CCTK_REAL_VEC pm1o256dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00390625),dy); const CCTK_REAL_VEC pm1o256dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.00390625),dz); const CCTK_REAL_VEC pm1o2dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.5),dx); const CCTK_REAL_VEC pm1o2dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.5),dy); const CCTK_REAL_VEC pm1o2dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.5),dz); const CCTK_REAL_VEC pm1o4dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.25),dx); const CCTK_REAL_VEC pm1o4dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.25),dy); const CCTK_REAL_VEC pm1o4dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.25),dz); const CCTK_REAL_VEC pm1o60dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0166666666666666666666666666667),dx); const CCTK_REAL_VEC pm1o60dy CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0166666666666666666666666666667),dy); const CCTK_REAL_VEC pm1o60dz CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(-0.0166666666666666666666666666667),dz); /* Jacobian variable pointers */ 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, GenericFD::jacobian_determinant_group, GenericFD::jacobian_inverse_group, and GenericFD::jacobian_derivative_group must all be set to valid group names"); } const CCTK_REAL* restrict jacobian_ptrs[9]; if (use_jacobian) GenericFD_GroupDataPointers(cctkGH, jacobian_group, 9, jacobian_ptrs); 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; 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); 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 */ /* Calculate temporaries and arrays functions */ /* Copy local copies back to grid functions */ /* Loop over the grid points */ 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(ML_CCZ4_Dissipation, i,j,k, imin0,imin1,imin2, imax0,imax1,imax2, cctk_ash[0],cctk_ash[1],cctk_ash[2], vecimin,vecimax, CCTK_REAL_VEC_SIZE) { 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 */ CCTK_REAL_VEC AL CCTK_ATTRIBUTE_UNUSED = vec_load(A[index]); CCTK_REAL_VEC alphaL CCTK_ATTRIBUTE_UNUSED = vec_load(alpha[index]); CCTK_REAL_VEC alpharhsL CCTK_ATTRIBUTE_UNUSED = vec_load(alpharhs[index]); CCTK_REAL_VEC ArhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Arhs[index]); CCTK_REAL_VEC At11L CCTK_ATTRIBUTE_UNUSED = vec_load(At11[index]); CCTK_REAL_VEC At11rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At11rhs[index]); CCTK_REAL_VEC At12L CCTK_ATTRIBUTE_UNUSED = vec_load(At12[index]); CCTK_REAL_VEC At12rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At12rhs[index]); CCTK_REAL_VEC At13L CCTK_ATTRIBUTE_UNUSED = vec_load(At13[index]); CCTK_REAL_VEC At13rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At13rhs[index]); CCTK_REAL_VEC At22L CCTK_ATTRIBUTE_UNUSED = vec_load(At22[index]); CCTK_REAL_VEC At22rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At22rhs[index]); CCTK_REAL_VEC At23L CCTK_ATTRIBUTE_UNUSED = vec_load(At23[index]); CCTK_REAL_VEC At23rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At23rhs[index]); CCTK_REAL_VEC At33L CCTK_ATTRIBUTE_UNUSED = vec_load(At33[index]); CCTK_REAL_VEC At33rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At33rhs[index]); CCTK_REAL_VEC B1L CCTK_ATTRIBUTE_UNUSED = vec_load(B1[index]); CCTK_REAL_VEC B1rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(B1rhs[index]); CCTK_REAL_VEC B2L CCTK_ATTRIBUTE_UNUSED = vec_load(B2[index]); CCTK_REAL_VEC B2rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(B2rhs[index]); CCTK_REAL_VEC B3L CCTK_ATTRIBUTE_UNUSED = vec_load(B3[index]); CCTK_REAL_VEC B3rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(B3rhs[index]); CCTK_REAL_VEC beta1L CCTK_ATTRIBUTE_UNUSED = vec_load(beta1[index]); CCTK_REAL_VEC beta1rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(beta1rhs[index]); CCTK_REAL_VEC beta2L CCTK_ATTRIBUTE_UNUSED = vec_load(beta2[index]); CCTK_REAL_VEC beta2rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(beta2rhs[index]); CCTK_REAL_VEC beta3L CCTK_ATTRIBUTE_UNUSED = vec_load(beta3[index]); CCTK_REAL_VEC beta3rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(beta3rhs[index]); CCTK_REAL_VEC gt11L CCTK_ATTRIBUTE_UNUSED = vec_load(gt11[index]); CCTK_REAL_VEC gt11rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt11rhs[index]); CCTK_REAL_VEC gt12L CCTK_ATTRIBUTE_UNUSED = vec_load(gt12[index]); CCTK_REAL_VEC gt12rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt12rhs[index]); CCTK_REAL_VEC gt13L CCTK_ATTRIBUTE_UNUSED = vec_load(gt13[index]); CCTK_REAL_VEC gt13rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt13rhs[index]); CCTK_REAL_VEC gt22L CCTK_ATTRIBUTE_UNUSED = vec_load(gt22[index]); CCTK_REAL_VEC gt22rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt22rhs[index]); CCTK_REAL_VEC gt23L CCTK_ATTRIBUTE_UNUSED = vec_load(gt23[index]); CCTK_REAL_VEC gt23rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt23rhs[index]); CCTK_REAL_VEC gt33L CCTK_ATTRIBUTE_UNUSED = vec_load(gt33[index]); CCTK_REAL_VEC gt33rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(gt33rhs[index]); CCTK_REAL_VEC phiL CCTK_ATTRIBUTE_UNUSED = vec_load(phi[index]); CCTK_REAL_VEC phirhsL CCTK_ATTRIBUTE_UNUSED = vec_load(phirhs[index]); CCTK_REAL_VEC ThetaL CCTK_ATTRIBUTE_UNUSED = vec_load(Theta[index]); CCTK_REAL_VEC ThetarhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Thetarhs[index]); CCTK_REAL_VEC trKL CCTK_ATTRIBUTE_UNUSED = vec_load(trK[index]); CCTK_REAL_VEC trKrhsL CCTK_ATTRIBUTE_UNUSED = vec_load(trKrhs[index]); CCTK_REAL_VEC Xt1L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt1[index]); CCTK_REAL_VEC Xt1rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Xt1rhs[index]); CCTK_REAL_VEC Xt2L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt2[index]); CCTK_REAL_VEC Xt2rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Xt2rhs[index]); CCTK_REAL_VEC Xt3L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt3[index]); CCTK_REAL_VEC Xt3rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Xt3rhs[index]); CCTK_REAL_VEC J11L, J12L, J13L, J21L, J22L, J23L, J31L, J32L, J33L CCTK_ATTRIBUTE_UNUSED ; if (use_jacobian) { J11L = vec_load(J11[index]); J12L = vec_load(J12[index]); J13L = vec_load(J13[index]); J21L = vec_load(J21[index]); J22L = vec_load(J22[index]); J23L = vec_load(J23[index]); J31L = vec_load(J31[index]); J32L = vec_load(J32[index]); J33L = vec_load(J33[index]); } /* Include user supplied include files */ /* Precompute derivatives */ CCTK_REAL_VEC PDdissipationNth1A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1Xt3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2Xt3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3Xt3 CCTK_ATTRIBUTE_UNUSED; switch (fdOrder) { case 2: { PDdissipationNth1A = PDdissipationNthfdOrder21(&A[index]); PDdissipationNth2A = PDdissipationNthfdOrder22(&A[index]); PDdissipationNth3A = PDdissipationNthfdOrder23(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder21(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder22(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder23(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder21(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder22(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder23(&At11[index]); PDdissipationNth1At12 = PDdissipationNthfdOrder21(&At12[index]); PDdissipationNth2At12 = PDdissipationNthfdOrder22(&At12[index]); PDdissipationNth3At12 = PDdissipationNthfdOrder23(&At12[index]); PDdissipationNth1At13 = PDdissipationNthfdOrder21(&At13[index]); PDdissipationNth2At13 = PDdissipationNthfdOrder22(&At13[index]); PDdissipationNth3At13 = PDdissipationNthfdOrder23(&At13[index]); PDdissipationNth1At22 = PDdissipationNthfdOrder21(&At22[index]); PDdissipationNth2At22 = PDdissipationNthfdOrder22(&At22[index]); PDdissipationNth3At22 = PDdissipationNthfdOrder23(&At22[index]); PDdissipationNth1At23 = PDdissipationNthfdOrder21(&At23[index]); PDdissipationNth2At23 = PDdissipationNthfdOrder22(&At23[index]); PDdissipationNth3At23 = PDdissipationNthfdOrder23(&At23[index]); PDdissipationNth1At33 = PDdissipationNthfdOrder21(&At33[index]); PDdissipationNth2At33 = PDdissipationNthfdOrder22(&At33[index]); PDdissipationNth3At33 = PDdissipationNthfdOrder23(&At33[index]); PDdissipationNth1B1 = PDdissipationNthfdOrder21(&B1[index]); PDdissipationNth2B1 = PDdissipationNthfdOrder22(&B1[index]); PDdissipationNth3B1 = PDdissipationNthfdOrder23(&B1[index]); PDdissipationNth1B2 = PDdissipationNthfdOrder21(&B2[index]); PDdissipationNth2B2 = PDdissipationNthfdOrder22(&B2[index]); PDdissipationNth3B2 = PDdissipationNthfdOrder23(&B2[index]); PDdissipationNth1B3 = PDdissipationNthfdOrder21(&B3[index]); PDdissipationNth2B3 = PDdissipationNthfdOrder22(&B3[index]); PDdissipationNth3B3 = PDdissipationNthfdOrder23(&B3[index]); PDdissipationNth1beta1 = PDdissipationNthfdOrder21(&beta1[index]); PDdissipationNth2beta1 = PDdissipationNthfdOrder22(&beta1[index]); PDdissipationNth3beta1 = PDdissipationNthfdOrder23(&beta1[index]); PDdissipationNth1beta2 = PDdissipationNthfdOrder21(&beta2[index]); PDdissipationNth2beta2 = PDdissipationNthfdOrder22(&beta2[index]); PDdissipationNth3beta2 = PDdissipationNthfdOrder23(&beta2[index]); PDdissipationNth1beta3 = PDdissipationNthfdOrder21(&beta3[index]); PDdissipationNth2beta3 = PDdissipationNthfdOrder22(&beta3[index]); PDdissipationNth3beta3 = PDdissipationNthfdOrder23(&beta3[index]); PDdissipationNth1gt11 = PDdissipationNthfdOrder21(>11[index]); PDdissipationNth2gt11 = PDdissipationNthfdOrder22(>11[index]); PDdissipationNth3gt11 = PDdissipationNthfdOrder23(>11[index]); PDdissipationNth1gt12 = PDdissipationNthfdOrder21(>12[index]); PDdissipationNth2gt12 = PDdissipationNthfdOrder22(>12[index]); PDdissipationNth3gt12 = PDdissipationNthfdOrder23(>12[index]); PDdissipationNth1gt13 = PDdissipationNthfdOrder21(>13[index]); PDdissipationNth2gt13 = PDdissipationNthfdOrder22(>13[index]); PDdissipationNth3gt13 = PDdissipationNthfdOrder23(>13[index]); PDdissipationNth1gt22 = PDdissipationNthfdOrder21(>22[index]); PDdissipationNth2gt22 = PDdissipationNthfdOrder22(>22[index]); PDdissipationNth3gt22 = PDdissipationNthfdOrder23(>22[index]); PDdissipationNth1gt23 = PDdissipationNthfdOrder21(>23[index]); PDdissipationNth2gt23 = PDdissipationNthfdOrder22(>23[index]); PDdissipationNth3gt23 = PDdissipationNthfdOrder23(>23[index]); PDdissipationNth1gt33 = PDdissipationNthfdOrder21(>33[index]); PDdissipationNth2gt33 = PDdissipationNthfdOrder22(>33[index]); PDdissipationNth3gt33 = PDdissipationNthfdOrder23(>33[index]); PDdissipationNth1phi = PDdissipationNthfdOrder21(&phi[index]); PDdissipationNth2phi = PDdissipationNthfdOrder22(&phi[index]); PDdissipationNth3phi = PDdissipationNthfdOrder23(&phi[index]); PDdissipationNth1Theta = PDdissipationNthfdOrder21(&Theta[index]); PDdissipationNth2Theta = PDdissipationNthfdOrder22(&Theta[index]); PDdissipationNth3Theta = PDdissipationNthfdOrder23(&Theta[index]); PDdissipationNth1trK = PDdissipationNthfdOrder21(&trK[index]); PDdissipationNth2trK = PDdissipationNthfdOrder22(&trK[index]); PDdissipationNth3trK = PDdissipationNthfdOrder23(&trK[index]); PDdissipationNth1Xt1 = PDdissipationNthfdOrder21(&Xt1[index]); PDdissipationNth2Xt1 = PDdissipationNthfdOrder22(&Xt1[index]); PDdissipationNth3Xt1 = PDdissipationNthfdOrder23(&Xt1[index]); PDdissipationNth1Xt2 = PDdissipationNthfdOrder21(&Xt2[index]); PDdissipationNth2Xt2 = PDdissipationNthfdOrder22(&Xt2[index]); PDdissipationNth3Xt2 = PDdissipationNthfdOrder23(&Xt2[index]); PDdissipationNth1Xt3 = PDdissipationNthfdOrder21(&Xt3[index]); PDdissipationNth2Xt3 = PDdissipationNthfdOrder22(&Xt3[index]); PDdissipationNth3Xt3 = PDdissipationNthfdOrder23(&Xt3[index]); break; } case 4: { PDdissipationNth1A = PDdissipationNthfdOrder41(&A[index]); PDdissipationNth2A = PDdissipationNthfdOrder42(&A[index]); PDdissipationNth3A = PDdissipationNthfdOrder43(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder41(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder42(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder43(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder41(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder42(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder43(&At11[index]); PDdissipationNth1At12 = PDdissipationNthfdOrder41(&At12[index]); PDdissipationNth2At12 = PDdissipationNthfdOrder42(&At12[index]); PDdissipationNth3At12 = PDdissipationNthfdOrder43(&At12[index]); PDdissipationNth1At13 = PDdissipationNthfdOrder41(&At13[index]); PDdissipationNth2At13 = PDdissipationNthfdOrder42(&At13[index]); PDdissipationNth3At13 = PDdissipationNthfdOrder43(&At13[index]); PDdissipationNth1At22 = PDdissipationNthfdOrder41(&At22[index]); PDdissipationNth2At22 = PDdissipationNthfdOrder42(&At22[index]); PDdissipationNth3At22 = PDdissipationNthfdOrder43(&At22[index]); PDdissipationNth1At23 = PDdissipationNthfdOrder41(&At23[index]); PDdissipationNth2At23 = PDdissipationNthfdOrder42(&At23[index]); PDdissipationNth3At23 = PDdissipationNthfdOrder43(&At23[index]); PDdissipationNth1At33 = PDdissipationNthfdOrder41(&At33[index]); PDdissipationNth2At33 = PDdissipationNthfdOrder42(&At33[index]); PDdissipationNth3At33 = PDdissipationNthfdOrder43(&At33[index]); PDdissipationNth1B1 = PDdissipationNthfdOrder41(&B1[index]); PDdissipationNth2B1 = PDdissipationNthfdOrder42(&B1[index]); PDdissipationNth3B1 = PDdissipationNthfdOrder43(&B1[index]); PDdissipationNth1B2 = PDdissipationNthfdOrder41(&B2[index]); PDdissipationNth2B2 = PDdissipationNthfdOrder42(&B2[index]); PDdissipationNth3B2 = PDdissipationNthfdOrder43(&B2[index]); PDdissipationNth1B3 = PDdissipationNthfdOrder41(&B3[index]); PDdissipationNth2B3 = PDdissipationNthfdOrder42(&B3[index]); PDdissipationNth3B3 = PDdissipationNthfdOrder43(&B3[index]); PDdissipationNth1beta1 = PDdissipationNthfdOrder41(&beta1[index]); PDdissipationNth2beta1 = PDdissipationNthfdOrder42(&beta1[index]); PDdissipationNth3beta1 = PDdissipationNthfdOrder43(&beta1[index]); PDdissipationNth1beta2 = PDdissipationNthfdOrder41(&beta2[index]); PDdissipationNth2beta2 = PDdissipationNthfdOrder42(&beta2[index]); PDdissipationNth3beta2 = PDdissipationNthfdOrder43(&beta2[index]); PDdissipationNth1beta3 = PDdissipationNthfdOrder41(&beta3[index]); PDdissipationNth2beta3 = PDdissipationNthfdOrder42(&beta3[index]); PDdissipationNth3beta3 = PDdissipationNthfdOrder43(&beta3[index]); PDdissipationNth1gt11 = PDdissipationNthfdOrder41(>11[index]); PDdissipationNth2gt11 = PDdissipationNthfdOrder42(>11[index]); PDdissipationNth3gt11 = PDdissipationNthfdOrder43(>11[index]); PDdissipationNth1gt12 = PDdissipationNthfdOrder41(>12[index]); PDdissipationNth2gt12 = PDdissipationNthfdOrder42(>12[index]); PDdissipationNth3gt12 = PDdissipationNthfdOrder43(>12[index]); PDdissipationNth1gt13 = PDdissipationNthfdOrder41(>13[index]); PDdissipationNth2gt13 = PDdissipationNthfdOrder42(>13[index]); PDdissipationNth3gt13 = PDdissipationNthfdOrder43(>13[index]); PDdissipationNth1gt22 = PDdissipationNthfdOrder41(>22[index]); PDdissipationNth2gt22 = PDdissipationNthfdOrder42(>22[index]); PDdissipationNth3gt22 = PDdissipationNthfdOrder43(>22[index]); PDdissipationNth1gt23 = PDdissipationNthfdOrder41(>23[index]); PDdissipationNth2gt23 = PDdissipationNthfdOrder42(>23[index]); PDdissipationNth3gt23 = PDdissipationNthfdOrder43(>23[index]); PDdissipationNth1gt33 = PDdissipationNthfdOrder41(>33[index]); PDdissipationNth2gt33 = PDdissipationNthfdOrder42(>33[index]); PDdissipationNth3gt33 = PDdissipationNthfdOrder43(>33[index]); PDdissipationNth1phi = PDdissipationNthfdOrder41(&phi[index]); PDdissipationNth2phi = PDdissipationNthfdOrder42(&phi[index]); PDdissipationNth3phi = PDdissipationNthfdOrder43(&phi[index]); PDdissipationNth1Theta = PDdissipationNthfdOrder41(&Theta[index]); PDdissipationNth2Theta = PDdissipationNthfdOrder42(&Theta[index]); PDdissipationNth3Theta = PDdissipationNthfdOrder43(&Theta[index]); PDdissipationNth1trK = PDdissipationNthfdOrder41(&trK[index]); PDdissipationNth2trK = PDdissipationNthfdOrder42(&trK[index]); PDdissipationNth3trK = PDdissipationNthfdOrder43(&trK[index]); PDdissipationNth1Xt1 = PDdissipationNthfdOrder41(&Xt1[index]); PDdissipationNth2Xt1 = PDdissipationNthfdOrder42(&Xt1[index]); PDdissipationNth3Xt1 = PDdissipationNthfdOrder43(&Xt1[index]); PDdissipationNth1Xt2 = PDdissipationNthfdOrder41(&Xt2[index]); PDdissipationNth2Xt2 = PDdissipationNthfdOrder42(&Xt2[index]); PDdissipationNth3Xt2 = PDdissipationNthfdOrder43(&Xt2[index]); PDdissipationNth1Xt3 = PDdissipationNthfdOrder41(&Xt3[index]); PDdissipationNth2Xt3 = PDdissipationNthfdOrder42(&Xt3[index]); PDdissipationNth3Xt3 = PDdissipationNthfdOrder43(&Xt3[index]); break; } case 6: { PDdissipationNth1A = PDdissipationNthfdOrder61(&A[index]); PDdissipationNth2A = PDdissipationNthfdOrder62(&A[index]); PDdissipationNth3A = PDdissipationNthfdOrder63(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder61(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder62(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder63(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder61(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder62(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder63(&At11[index]); PDdissipationNth1At12 = PDdissipationNthfdOrder61(&At12[index]); PDdissipationNth2At12 = PDdissipationNthfdOrder62(&At12[index]); PDdissipationNth3At12 = PDdissipationNthfdOrder63(&At12[index]); PDdissipationNth1At13 = PDdissipationNthfdOrder61(&At13[index]); PDdissipationNth2At13 = PDdissipationNthfdOrder62(&At13[index]); PDdissipationNth3At13 = PDdissipationNthfdOrder63(&At13[index]); PDdissipationNth1At22 = PDdissipationNthfdOrder61(&At22[index]); PDdissipationNth2At22 = PDdissipationNthfdOrder62(&At22[index]); PDdissipationNth3At22 = PDdissipationNthfdOrder63(&At22[index]); PDdissipationNth1At23 = PDdissipationNthfdOrder61(&At23[index]); PDdissipationNth2At23 = PDdissipationNthfdOrder62(&At23[index]); PDdissipationNth3At23 = PDdissipationNthfdOrder63(&At23[index]); PDdissipationNth1At33 = PDdissipationNthfdOrder61(&At33[index]); PDdissipationNth2At33 = PDdissipationNthfdOrder62(&At33[index]); PDdissipationNth3At33 = PDdissipationNthfdOrder63(&At33[index]); PDdissipationNth1B1 = PDdissipationNthfdOrder61(&B1[index]); PDdissipationNth2B1 = PDdissipationNthfdOrder62(&B1[index]); PDdissipationNth3B1 = PDdissipationNthfdOrder63(&B1[index]); PDdissipationNth1B2 = PDdissipationNthfdOrder61(&B2[index]); PDdissipationNth2B2 = PDdissipationNthfdOrder62(&B2[index]); PDdissipationNth3B2 = PDdissipationNthfdOrder63(&B2[index]); PDdissipationNth1B3 = PDdissipationNthfdOrder61(&B3[index]); PDdissipationNth2B3 = PDdissipationNthfdOrder62(&B3[index]); PDdissipationNth3B3 = PDdissipationNthfdOrder63(&B3[index]); PDdissipationNth1beta1 = PDdissipationNthfdOrder61(&beta1[index]); PDdissipationNth2beta1 = PDdissipationNthfdOrder62(&beta1[index]); PDdissipationNth3beta1 = PDdissipationNthfdOrder63(&beta1[index]); PDdissipationNth1beta2 = PDdissipationNthfdOrder61(&beta2[index]); PDdissipationNth2beta2 = PDdissipationNthfdOrder62(&beta2[index]); PDdissipationNth3beta2 = PDdissipationNthfdOrder63(&beta2[index]); PDdissipationNth1beta3 = PDdissipationNthfdOrder61(&beta3[index]); PDdissipationNth2beta3 = PDdissipationNthfdOrder62(&beta3[index]); PDdissipationNth3beta3 = PDdissipationNthfdOrder63(&beta3[index]); PDdissipationNth1gt11 = PDdissipationNthfdOrder61(>11[index]); PDdissipationNth2gt11 = PDdissipationNthfdOrder62(>11[index]); PDdissipationNth3gt11 = PDdissipationNthfdOrder63(>11[index]); PDdissipationNth1gt12 = PDdissipationNthfdOrder61(>12[index]); PDdissipationNth2gt12 = PDdissipationNthfdOrder62(>12[index]); PDdissipationNth3gt12 = PDdissipationNthfdOrder63(>12[index]); PDdissipationNth1gt13 = PDdissipationNthfdOrder61(>13[index]); PDdissipationNth2gt13 = PDdissipationNthfdOrder62(>13[index]); PDdissipationNth3gt13 = PDdissipationNthfdOrder63(>13[index]); PDdissipationNth1gt22 = PDdissipationNthfdOrder61(>22[index]); PDdissipationNth2gt22 = PDdissipationNthfdOrder62(>22[index]); PDdissipationNth3gt22 = PDdissipationNthfdOrder63(>22[index]); PDdissipationNth1gt23 = PDdissipationNthfdOrder61(>23[index]); PDdissipationNth2gt23 = PDdissipationNthfdOrder62(>23[index]); PDdissipationNth3gt23 = PDdissipationNthfdOrder63(>23[index]); PDdissipationNth1gt33 = PDdissipationNthfdOrder61(>33[index]); PDdissipationNth2gt33 = PDdissipationNthfdOrder62(>33[index]); PDdissipationNth3gt33 = PDdissipationNthfdOrder63(>33[index]); PDdissipationNth1phi = PDdissipationNthfdOrder61(&phi[index]); PDdissipationNth2phi = PDdissipationNthfdOrder62(&phi[index]); PDdissipationNth3phi = PDdissipationNthfdOrder63(&phi[index]); PDdissipationNth1Theta = PDdissipationNthfdOrder61(&Theta[index]); PDdissipationNth2Theta = PDdissipationNthfdOrder62(&Theta[index]); PDdissipationNth3Theta = PDdissipationNthfdOrder63(&Theta[index]); PDdissipationNth1trK = PDdissipationNthfdOrder61(&trK[index]); PDdissipationNth2trK = PDdissipationNthfdOrder62(&trK[index]); PDdissipationNth3trK = PDdissipationNthfdOrder63(&trK[index]); PDdissipationNth1Xt1 = PDdissipationNthfdOrder61(&Xt1[index]); PDdissipationNth2Xt1 = PDdissipationNthfdOrder62(&Xt1[index]); PDdissipationNth3Xt1 = PDdissipationNthfdOrder63(&Xt1[index]); PDdissipationNth1Xt2 = PDdissipationNthfdOrder61(&Xt2[index]); PDdissipationNth2Xt2 = PDdissipationNthfdOrder62(&Xt2[index]); PDdissipationNth3Xt2 = PDdissipationNthfdOrder63(&Xt2[index]); PDdissipationNth1Xt3 = PDdissipationNthfdOrder61(&Xt3[index]); PDdissipationNth2Xt3 = PDdissipationNthfdOrder62(&Xt3[index]); PDdissipationNth3Xt3 = PDdissipationNthfdOrder63(&Xt3[index]); break; } case 8: { PDdissipationNth1A = PDdissipationNthfdOrder81(&A[index]); PDdissipationNth2A = PDdissipationNthfdOrder82(&A[index]); PDdissipationNth3A = PDdissipationNthfdOrder83(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder81(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder82(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder83(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder81(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder82(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder83(&At11[index]); PDdissipationNth1At12 = PDdissipationNthfdOrder81(&At12[index]); PDdissipationNth2At12 = PDdissipationNthfdOrder82(&At12[index]); PDdissipationNth3At12 = PDdissipationNthfdOrder83(&At12[index]); PDdissipationNth1At13 = PDdissipationNthfdOrder81(&At13[index]); PDdissipationNth2At13 = PDdissipationNthfdOrder82(&At13[index]); PDdissipationNth3At13 = PDdissipationNthfdOrder83(&At13[index]); PDdissipationNth1At22 = PDdissipationNthfdOrder81(&At22[index]); PDdissipationNth2At22 = PDdissipationNthfdOrder82(&At22[index]); PDdissipationNth3At22 = PDdissipationNthfdOrder83(&At22[index]); PDdissipationNth1At23 = PDdissipationNthfdOrder81(&At23[index]); PDdissipationNth2At23 = PDdissipationNthfdOrder82(&At23[index]); PDdissipationNth3At23 = PDdissipationNthfdOrder83(&At23[index]); PDdissipationNth1At33 = PDdissipationNthfdOrder81(&At33[index]); PDdissipationNth2At33 = PDdissipationNthfdOrder82(&At33[index]); PDdissipationNth3At33 = PDdissipationNthfdOrder83(&At33[index]); PDdissipationNth1B1 = PDdissipationNthfdOrder81(&B1[index]); PDdissipationNth2B1 = PDdissipationNthfdOrder82(&B1[index]); PDdissipationNth3B1 = PDdissipationNthfdOrder83(&B1[index]); PDdissipationNth1B2 = PDdissipationNthfdOrder81(&B2[index]); PDdissipationNth2B2 = PDdissipationNthfdOrder82(&B2[index]); PDdissipationNth3B2 = PDdissipationNthfdOrder83(&B2[index]); PDdissipationNth1B3 = PDdissipationNthfdOrder81(&B3[index]); PDdissipationNth2B3 = PDdissipationNthfdOrder82(&B3[index]); PDdissipationNth3B3 = PDdissipationNthfdOrder83(&B3[index]); PDdissipationNth1beta1 = PDdissipationNthfdOrder81(&beta1[index]); PDdissipationNth2beta1 = PDdissipationNthfdOrder82(&beta1[index]); PDdissipationNth3beta1 = PDdissipationNthfdOrder83(&beta1[index]); PDdissipationNth1beta2 = PDdissipationNthfdOrder81(&beta2[index]); PDdissipationNth2beta2 = PDdissipationNthfdOrder82(&beta2[index]); PDdissipationNth3beta2 = PDdissipationNthfdOrder83(&beta2[index]); PDdissipationNth1beta3 = PDdissipationNthfdOrder81(&beta3[index]); PDdissipationNth2beta3 = PDdissipationNthfdOrder82(&beta3[index]); PDdissipationNth3beta3 = PDdissipationNthfdOrder83(&beta3[index]); PDdissipationNth1gt11 = PDdissipationNthfdOrder81(>11[index]); PDdissipationNth2gt11 = PDdissipationNthfdOrder82(>11[index]); PDdissipationNth3gt11 = PDdissipationNthfdOrder83(>11[index]); PDdissipationNth1gt12 = PDdissipationNthfdOrder81(>12[index]); PDdissipationNth2gt12 = PDdissipationNthfdOrder82(>12[index]); PDdissipationNth3gt12 = PDdissipationNthfdOrder83(>12[index]); PDdissipationNth1gt13 = PDdissipationNthfdOrder81(>13[index]); PDdissipationNth2gt13 = PDdissipationNthfdOrder82(>13[index]); PDdissipationNth3gt13 = PDdissipationNthfdOrder83(>13[index]); PDdissipationNth1gt22 = PDdissipationNthfdOrder81(>22[index]); PDdissipationNth2gt22 = PDdissipationNthfdOrder82(>22[index]); PDdissipationNth3gt22 = PDdissipationNthfdOrder83(>22[index]); PDdissipationNth1gt23 = PDdissipationNthfdOrder81(>23[index]); PDdissipationNth2gt23 = PDdissipationNthfdOrder82(>23[index]); PDdissipationNth3gt23 = PDdissipationNthfdOrder83(>23[index]); PDdissipationNth1gt33 = PDdissipationNthfdOrder81(>33[index]); PDdissipationNth2gt33 = PDdissipationNthfdOrder82(>33[index]); PDdissipationNth3gt33 = PDdissipationNthfdOrder83(>33[index]); PDdissipationNth1phi = PDdissipationNthfdOrder81(&phi[index]); PDdissipationNth2phi = PDdissipationNthfdOrder82(&phi[index]); PDdissipationNth3phi = PDdissipationNthfdOrder83(&phi[index]); PDdissipationNth1Theta = PDdissipationNthfdOrder81(&Theta[index]); PDdissipationNth2Theta = PDdissipationNthfdOrder82(&Theta[index]); PDdissipationNth3Theta = PDdissipationNthfdOrder83(&Theta[index]); PDdissipationNth1trK = PDdissipationNthfdOrder81(&trK[index]); PDdissipationNth2trK = PDdissipationNthfdOrder82(&trK[index]); PDdissipationNth3trK = PDdissipationNthfdOrder83(&trK[index]); PDdissipationNth1Xt1 = PDdissipationNthfdOrder81(&Xt1[index]); PDdissipationNth2Xt1 = PDdissipationNthfdOrder82(&Xt1[index]); PDdissipationNth3Xt1 = PDdissipationNthfdOrder83(&Xt1[index]); PDdissipationNth1Xt2 = PDdissipationNthfdOrder81(&Xt2[index]); PDdissipationNth2Xt2 = PDdissipationNthfdOrder82(&Xt2[index]); PDdissipationNth3Xt2 = PDdissipationNthfdOrder83(&Xt2[index]); PDdissipationNth1Xt3 = PDdissipationNthfdOrder81(&Xt3[index]); PDdissipationNth2Xt3 = PDdissipationNthfdOrder82(&Xt3[index]); PDdissipationNth3Xt3 = PDdissipationNthfdOrder83(&Xt3[index]); break; } default: CCTK_BUILTIN_UNREACHABLE(); } /* Calculate temporaries and grid functions */ CCTK_REAL_VEC JacPDdissipationNth1A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Xt3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Xt3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3A CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3B1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3B2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3B3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3beta1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3beta2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3beta3 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt13 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt22 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt23 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3gt33 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3phi CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3Theta CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3trK CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3Xt3 CCTK_ATTRIBUTE_UNUSED; if (use_jacobian) { JacPDdissipationNth1A = kmadd(J11L,PDdissipationNth1A,kmadd(J21L,PDdissipationNth2A,kmul(J31L,PDdissipationNth3A))); JacPDdissipationNth1alpha = kmadd(J11L,PDdissipationNth1alpha,kmadd(J21L,PDdissipationNth2alpha,kmul(J31L,PDdissipationNth3alpha))); JacPDdissipationNth1At11 = kmadd(J11L,PDdissipationNth1At11,kmadd(J21L,PDdissipationNth2At11,kmul(J31L,PDdissipationNth3At11))); JacPDdissipationNth1At12 = kmadd(J11L,PDdissipationNth1At12,kmadd(J21L,PDdissipationNth2At12,kmul(J31L,PDdissipationNth3At12))); JacPDdissipationNth1At13 = kmadd(J11L,PDdissipationNth1At13,kmadd(J21L,PDdissipationNth2At13,kmul(J31L,PDdissipationNth3At13))); JacPDdissipationNth1At22 = kmadd(J11L,PDdissipationNth1At22,kmadd(J21L,PDdissipationNth2At22,kmul(J31L,PDdissipationNth3At22))); JacPDdissipationNth1At23 = kmadd(J11L,PDdissipationNth1At23,kmadd(J21L,PDdissipationNth2At23,kmul(J31L,PDdissipationNth3At23))); JacPDdissipationNth1At33 = kmadd(J11L,PDdissipationNth1At33,kmadd(J21L,PDdissipationNth2At33,kmul(J31L,PDdissipationNth3At33))); JacPDdissipationNth1B1 = kmadd(J11L,PDdissipationNth1B1,kmadd(J21L,PDdissipationNth2B1,kmul(J31L,PDdissipationNth3B1))); JacPDdissipationNth1B2 = kmadd(J11L,PDdissipationNth1B2,kmadd(J21L,PDdissipationNth2B2,kmul(J31L,PDdissipationNth3B2))); JacPDdissipationNth1B3 = kmadd(J11L,PDdissipationNth1B3,kmadd(J21L,PDdissipationNth2B3,kmul(J31L,PDdissipationNth3B3))); JacPDdissipationNth1beta1 = kmadd(J11L,PDdissipationNth1beta1,kmadd(J21L,PDdissipationNth2beta1,kmul(J31L,PDdissipationNth3beta1))); JacPDdissipationNth1beta2 = kmadd(J11L,PDdissipationNth1beta2,kmadd(J21L,PDdissipationNth2beta2,kmul(J31L,PDdissipationNth3beta2))); JacPDdissipationNth1beta3 = kmadd(J11L,PDdissipationNth1beta3,kmadd(J21L,PDdissipationNth2beta3,kmul(J31L,PDdissipationNth3beta3))); JacPDdissipationNth1gt11 = kmadd(J11L,PDdissipationNth1gt11,kmadd(J21L,PDdissipationNth2gt11,kmul(J31L,PDdissipationNth3gt11))); JacPDdissipationNth1gt12 = kmadd(J11L,PDdissipationNth1gt12,kmadd(J21L,PDdissipationNth2gt12,kmul(J31L,PDdissipationNth3gt12))); JacPDdissipationNth1gt13 = kmadd(J11L,PDdissipationNth1gt13,kmadd(J21L,PDdissipationNth2gt13,kmul(J31L,PDdissipationNth3gt13))); JacPDdissipationNth1gt22 = kmadd(J11L,PDdissipationNth1gt22,kmadd(J21L,PDdissipationNth2gt22,kmul(J31L,PDdissipationNth3gt22))); JacPDdissipationNth1gt23 = kmadd(J11L,PDdissipationNth1gt23,kmadd(J21L,PDdissipationNth2gt23,kmul(J31L,PDdissipationNth3gt23))); JacPDdissipationNth1gt33 = kmadd(J11L,PDdissipationNth1gt33,kmadd(J21L,PDdissipationNth2gt33,kmul(J31L,PDdissipationNth3gt33))); JacPDdissipationNth1phi = kmadd(J11L,PDdissipationNth1phi,kmadd(J21L,PDdissipationNth2phi,kmul(J31L,PDdissipationNth3phi))); JacPDdissipationNth1Theta = kmadd(J11L,PDdissipationNth1Theta,kmadd(J21L,PDdissipationNth2Theta,kmul(J31L,PDdissipationNth3Theta))); JacPDdissipationNth1trK = kmadd(J11L,PDdissipationNth1trK,kmadd(J21L,PDdissipationNth2trK,kmul(J31L,PDdissipationNth3trK))); JacPDdissipationNth1Xt1 = kmadd(J11L,PDdissipationNth1Xt1,kmadd(J21L,PDdissipationNth2Xt1,kmul(J31L,PDdissipationNth3Xt1))); JacPDdissipationNth1Xt2 = kmadd(J11L,PDdissipationNth1Xt2,kmadd(J21L,PDdissipationNth2Xt2,kmul(J31L,PDdissipationNth3Xt2))); JacPDdissipationNth1Xt3 = kmadd(J11L,PDdissipationNth1Xt3,kmadd(J21L,PDdissipationNth2Xt3,kmul(J31L,PDdissipationNth3Xt3))); JacPDdissipationNth2A = kmadd(J12L,PDdissipationNth1A,kmadd(J22L,PDdissipationNth2A,kmul(J32L,PDdissipationNth3A))); JacPDdissipationNth2alpha = kmadd(J12L,PDdissipationNth1alpha,kmadd(J22L,PDdissipationNth2alpha,kmul(J32L,PDdissipationNth3alpha))); JacPDdissipationNth2At11 = kmadd(J12L,PDdissipationNth1At11,kmadd(J22L,PDdissipationNth2At11,kmul(J32L,PDdissipationNth3At11))); JacPDdissipationNth2At12 = kmadd(J12L,PDdissipationNth1At12,kmadd(J22L,PDdissipationNth2At12,kmul(J32L,PDdissipationNth3At12))); JacPDdissipationNth2At13 = kmadd(J12L,PDdissipationNth1At13,kmadd(J22L,PDdissipationNth2At13,kmul(J32L,PDdissipationNth3At13))); JacPDdissipationNth2At22 = kmadd(J12L,PDdissipationNth1At22,kmadd(J22L,PDdissipationNth2At22,kmul(J32L,PDdissipationNth3At22))); JacPDdissipationNth2At23 = kmadd(J12L,PDdissipationNth1At23,kmadd(J22L,PDdissipationNth2At23,kmul(J32L,PDdissipationNth3At23))); JacPDdissipationNth2At33 = kmadd(J12L,PDdissipationNth1At33,kmadd(J22L,PDdissipationNth2At33,kmul(J32L,PDdissipationNth3At33))); JacPDdissipationNth2B1 = kmadd(J12L,PDdissipationNth1B1,kmadd(J22L,PDdissipationNth2B1,kmul(J32L,PDdissipationNth3B1))); JacPDdissipationNth2B2 = kmadd(J12L,PDdissipationNth1B2,kmadd(J22L,PDdissipationNth2B2,kmul(J32L,PDdissipationNth3B2))); JacPDdissipationNth2B3 = kmadd(J12L,PDdissipationNth1B3,kmadd(J22L,PDdissipationNth2B3,kmul(J32L,PDdissipationNth3B3))); JacPDdissipationNth2beta1 = kmadd(J12L,PDdissipationNth1beta1,kmadd(J22L,PDdissipationNth2beta1,kmul(J32L,PDdissipationNth3beta1))); JacPDdissipationNth2beta2 = kmadd(J12L,PDdissipationNth1beta2,kmadd(J22L,PDdissipationNth2beta2,kmul(J32L,PDdissipationNth3beta2))); JacPDdissipationNth2beta3 = kmadd(J12L,PDdissipationNth1beta3,kmadd(J22L,PDdissipationNth2beta3,kmul(J32L,PDdissipationNth3beta3))); JacPDdissipationNth2gt11 = kmadd(J12L,PDdissipationNth1gt11,kmadd(J22L,PDdissipationNth2gt11,kmul(J32L,PDdissipationNth3gt11))); JacPDdissipationNth2gt12 = kmadd(J12L,PDdissipationNth1gt12,kmadd(J22L,PDdissipationNth2gt12,kmul(J32L,PDdissipationNth3gt12))); JacPDdissipationNth2gt13 = kmadd(J12L,PDdissipationNth1gt13,kmadd(J22L,PDdissipationNth2gt13,kmul(J32L,PDdissipationNth3gt13))); JacPDdissipationNth2gt22 = kmadd(J12L,PDdissipationNth1gt22,kmadd(J22L,PDdissipationNth2gt22,kmul(J32L,PDdissipationNth3gt22))); JacPDdissipationNth2gt23 = kmadd(J12L,PDdissipationNth1gt23,kmadd(J22L,PDdissipationNth2gt23,kmul(J32L,PDdissipationNth3gt23))); JacPDdissipationNth2gt33 = kmadd(J12L,PDdissipationNth1gt33,kmadd(J22L,PDdissipationNth2gt33,kmul(J32L,PDdissipationNth3gt33))); JacPDdissipationNth2phi = kmadd(J12L,PDdissipationNth1phi,kmadd(J22L,PDdissipationNth2phi,kmul(J32L,PDdissipationNth3phi))); JacPDdissipationNth2Theta = kmadd(J12L,PDdissipationNth1Theta,kmadd(J22L,PDdissipationNth2Theta,kmul(J32L,PDdissipationNth3Theta))); JacPDdissipationNth2trK = kmadd(J12L,PDdissipationNth1trK,kmadd(J22L,PDdissipationNth2trK,kmul(J32L,PDdissipationNth3trK))); JacPDdissipationNth2Xt1 = kmadd(J12L,PDdissipationNth1Xt1,kmadd(J22L,PDdissipationNth2Xt1,kmul(J32L,PDdissipationNth3Xt1))); JacPDdissipationNth2Xt2 = kmadd(J12L,PDdissipationNth1Xt2,kmadd(J22L,PDdissipationNth2Xt2,kmul(J32L,PDdissipationNth3Xt2))); JacPDdissipationNth2Xt3 = kmadd(J12L,PDdissipationNth1Xt3,kmadd(J22L,PDdissipationNth2Xt3,kmul(J32L,PDdissipationNth3Xt3))); JacPDdissipationNth3A = kmadd(J13L,PDdissipationNth1A,kmadd(J23L,PDdissipationNth2A,kmul(J33L,PDdissipationNth3A))); JacPDdissipationNth3alpha = kmadd(J13L,PDdissipationNth1alpha,kmadd(J23L,PDdissipationNth2alpha,kmul(J33L,PDdissipationNth3alpha))); JacPDdissipationNth3At11 = kmadd(J13L,PDdissipationNth1At11,kmadd(J23L,PDdissipationNth2At11,kmul(J33L,PDdissipationNth3At11))); JacPDdissipationNth3At12 = kmadd(J13L,PDdissipationNth1At12,kmadd(J23L,PDdissipationNth2At12,kmul(J33L,PDdissipationNth3At12))); JacPDdissipationNth3At13 = kmadd(J13L,PDdissipationNth1At13,kmadd(J23L,PDdissipationNth2At13,kmul(J33L,PDdissipationNth3At13))); JacPDdissipationNth3At22 = kmadd(J13L,PDdissipationNth1At22,kmadd(J23L,PDdissipationNth2At22,kmul(J33L,PDdissipationNth3At22))); JacPDdissipationNth3At23 = kmadd(J13L,PDdissipationNth1At23,kmadd(J23L,PDdissipationNth2At23,kmul(J33L,PDdissipationNth3At23))); JacPDdissipationNth3At33 = kmadd(J13L,PDdissipationNth1At33,kmadd(J23L,PDdissipationNth2At33,kmul(J33L,PDdissipationNth3At33))); JacPDdissipationNth3B1 = kmadd(J13L,PDdissipationNth1B1,kmadd(J23L,PDdissipationNth2B1,kmul(J33L,PDdissipationNth3B1))); JacPDdissipationNth3B2 = kmadd(J13L,PDdissipationNth1B2,kmadd(J23L,PDdissipationNth2B2,kmul(J33L,PDdissipationNth3B2))); JacPDdissipationNth3B3 = kmadd(J13L,PDdissipationNth1B3,kmadd(J23L,PDdissipationNth2B3,kmul(J33L,PDdissipationNth3B3))); JacPDdissipationNth3beta1 = kmadd(J13L,PDdissipationNth1beta1,kmadd(J23L,PDdissipationNth2beta1,kmul(J33L,PDdissipationNth3beta1))); JacPDdissipationNth3beta2 = kmadd(J13L,PDdissipationNth1beta2,kmadd(J23L,PDdissipationNth2beta2,kmul(J33L,PDdissipationNth3beta2))); JacPDdissipationNth3beta3 = kmadd(J13L,PDdissipationNth1beta3,kmadd(J23L,PDdissipationNth2beta3,kmul(J33L,PDdissipationNth3beta3))); JacPDdissipationNth3gt11 = kmadd(J13L,PDdissipationNth1gt11,kmadd(J23L,PDdissipationNth2gt11,kmul(J33L,PDdissipationNth3gt11))); JacPDdissipationNth3gt12 = kmadd(J13L,PDdissipationNth1gt12,kmadd(J23L,PDdissipationNth2gt12,kmul(J33L,PDdissipationNth3gt12))); JacPDdissipationNth3gt13 = kmadd(J13L,PDdissipationNth1gt13,kmadd(J23L,PDdissipationNth2gt13,kmul(J33L,PDdissipationNth3gt13))); JacPDdissipationNth3gt22 = kmadd(J13L,PDdissipationNth1gt22,kmadd(J23L,PDdissipationNth2gt22,kmul(J33L,PDdissipationNth3gt22))); JacPDdissipationNth3gt23 = kmadd(J13L,PDdissipationNth1gt23,kmadd(J23L,PDdissipationNth2gt23,kmul(J33L,PDdissipationNth3gt23))); JacPDdissipationNth3gt33 = kmadd(J13L,PDdissipationNth1gt33,kmadd(J23L,PDdissipationNth2gt33,kmul(J33L,PDdissipationNth3gt33))); JacPDdissipationNth3phi = kmadd(J13L,PDdissipationNth1phi,kmadd(J23L,PDdissipationNth2phi,kmul(J33L,PDdissipationNth3phi))); JacPDdissipationNth3Theta = kmadd(J13L,PDdissipationNth1Theta,kmadd(J23L,PDdissipationNth2Theta,kmul(J33L,PDdissipationNth3Theta))); JacPDdissipationNth3trK = kmadd(J13L,PDdissipationNth1trK,kmadd(J23L,PDdissipationNth2trK,kmul(J33L,PDdissipationNth3trK))); JacPDdissipationNth3Xt1 = kmadd(J13L,PDdissipationNth1Xt1,kmadd(J23L,PDdissipationNth2Xt1,kmul(J33L,PDdissipationNth3Xt1))); JacPDdissipationNth3Xt2 = kmadd(J13L,PDdissipationNth1Xt2,kmadd(J23L,PDdissipationNth2Xt2,kmul(J33L,PDdissipationNth3Xt2))); JacPDdissipationNth3Xt3 = kmadd(J13L,PDdissipationNth1Xt3,kmadd(J23L,PDdissipationNth2Xt3,kmul(J33L,PDdissipationNth3Xt3))); } else { JacPDdissipationNth1A = PDdissipationNth1A; JacPDdissipationNth1alpha = PDdissipationNth1alpha; JacPDdissipationNth1At11 = PDdissipationNth1At11; JacPDdissipationNth1At12 = PDdissipationNth1At12; JacPDdissipationNth1At13 = PDdissipationNth1At13; JacPDdissipationNth1At22 = PDdissipationNth1At22; JacPDdissipationNth1At23 = PDdissipationNth1At23; JacPDdissipationNth1At33 = PDdissipationNth1At33; JacPDdissipationNth1B1 = PDdissipationNth1B1; JacPDdissipationNth1B2 = PDdissipationNth1B2; JacPDdissipationNth1B3 = PDdissipationNth1B3; JacPDdissipationNth1beta1 = PDdissipationNth1beta1; JacPDdissipationNth1beta2 = PDdissipationNth1beta2; JacPDdissipationNth1beta3 = PDdissipationNth1beta3; JacPDdissipationNth1gt11 = PDdissipationNth1gt11; JacPDdissipationNth1gt12 = PDdissipationNth1gt12; JacPDdissipationNth1gt13 = PDdissipationNth1gt13; JacPDdissipationNth1gt22 = PDdissipationNth1gt22; JacPDdissipationNth1gt23 = PDdissipationNth1gt23; JacPDdissipationNth1gt33 = PDdissipationNth1gt33; JacPDdissipationNth1phi = PDdissipationNth1phi; JacPDdissipationNth1Theta = PDdissipationNth1Theta; JacPDdissipationNth1trK = PDdissipationNth1trK; JacPDdissipationNth1Xt1 = PDdissipationNth1Xt1; JacPDdissipationNth1Xt2 = PDdissipationNth1Xt2; JacPDdissipationNth1Xt3 = PDdissipationNth1Xt3; JacPDdissipationNth2A = PDdissipationNth2A; JacPDdissipationNth2alpha = PDdissipationNth2alpha; JacPDdissipationNth2At11 = PDdissipationNth2At11; JacPDdissipationNth2At12 = PDdissipationNth2At12; JacPDdissipationNth2At13 = PDdissipationNth2At13; JacPDdissipationNth2At22 = PDdissipationNth2At22; JacPDdissipationNth2At23 = PDdissipationNth2At23; JacPDdissipationNth2At33 = PDdissipationNth2At33; JacPDdissipationNth2B1 = PDdissipationNth2B1; JacPDdissipationNth2B2 = PDdissipationNth2B2; JacPDdissipationNth2B3 = PDdissipationNth2B3; JacPDdissipationNth2beta1 = PDdissipationNth2beta1; JacPDdissipationNth2beta2 = PDdissipationNth2beta2; JacPDdissipationNth2beta3 = PDdissipationNth2beta3; JacPDdissipationNth2gt11 = PDdissipationNth2gt11; JacPDdissipationNth2gt12 = PDdissipationNth2gt12; JacPDdissipationNth2gt13 = PDdissipationNth2gt13; JacPDdissipationNth2gt22 = PDdissipationNth2gt22; JacPDdissipationNth2gt23 = PDdissipationNth2gt23; JacPDdissipationNth2gt33 = PDdissipationNth2gt33; JacPDdissipationNth2phi = PDdissipationNth2phi; JacPDdissipationNth2Theta = PDdissipationNth2Theta; JacPDdissipationNth2trK = PDdissipationNth2trK; JacPDdissipationNth2Xt1 = PDdissipationNth2Xt1; JacPDdissipationNth2Xt2 = PDdissipationNth2Xt2; JacPDdissipationNth2Xt3 = PDdissipationNth2Xt3; JacPDdissipationNth3A = PDdissipationNth3A; JacPDdissipationNth3alpha = PDdissipationNth3alpha; JacPDdissipationNth3At11 = PDdissipationNth3At11; JacPDdissipationNth3At12 = PDdissipationNth3At12; JacPDdissipationNth3At13 = PDdissipationNth3At13; JacPDdissipationNth3At22 = PDdissipationNth3At22; JacPDdissipationNth3At23 = PDdissipationNth3At23; JacPDdissipationNth3At33 = PDdissipationNth3At33; JacPDdissipationNth3B1 = PDdissipationNth3B1; JacPDdissipationNth3B2 = PDdissipationNth3B2; JacPDdissipationNth3B3 = PDdissipationNth3B3; JacPDdissipationNth3beta1 = PDdissipationNth3beta1; JacPDdissipationNth3beta2 = PDdissipationNth3beta2; JacPDdissipationNth3beta3 = PDdissipationNth3beta3; JacPDdissipationNth3gt11 = PDdissipationNth3gt11; JacPDdissipationNth3gt12 = PDdissipationNth3gt12; JacPDdissipationNth3gt13 = PDdissipationNth3gt13; JacPDdissipationNth3gt22 = PDdissipationNth3gt22; JacPDdissipationNth3gt23 = PDdissipationNth3gt23; JacPDdissipationNth3gt33 = PDdissipationNth3gt33; JacPDdissipationNth3phi = PDdissipationNth3phi; JacPDdissipationNth3Theta = PDdissipationNth3Theta; JacPDdissipationNth3trK = PDdissipationNth3trK; JacPDdissipationNth3Xt1 = PDdissipationNth3Xt1; JacPDdissipationNth3Xt2 = PDdissipationNth3Xt2; JacPDdissipationNth3Xt3 = PDdissipationNth3Xt3; } CCTK_REAL_VEC epsdiss1 CCTK_ATTRIBUTE_UNUSED = ToReal(EpsDiss); CCTK_REAL_VEC epsdiss2 CCTK_ATTRIBUTE_UNUSED = ToReal(EpsDiss); CCTK_REAL_VEC epsdiss3 CCTK_ATTRIBUTE_UNUSED = ToReal(EpsDiss); phirhsL = kadd(phirhsL,kmadd(epsdiss1,JacPDdissipationNth1phi,kmadd(epsdiss2,JacPDdissipationNth2phi,kmul(epsdiss3,JacPDdissipationNth3phi)))); gt11rhsL = kadd(gt11rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt11,kmadd(epsdiss2,JacPDdissipationNth2gt11,kmul(epsdiss3,JacPDdissipationNth3gt11)))); gt12rhsL = kadd(gt12rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt12,kmadd(epsdiss2,JacPDdissipationNth2gt12,kmul(epsdiss3,JacPDdissipationNth3gt12)))); gt13rhsL = kadd(gt13rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt13,kmadd(epsdiss2,JacPDdissipationNth2gt13,kmul(epsdiss3,JacPDdissipationNth3gt13)))); gt22rhsL = kadd(gt22rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt22,kmadd(epsdiss2,JacPDdissipationNth2gt22,kmul(epsdiss3,JacPDdissipationNth3gt22)))); gt23rhsL = kadd(gt23rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt23,kmadd(epsdiss2,JacPDdissipationNth2gt23,kmul(epsdiss3,JacPDdissipationNth3gt23)))); gt33rhsL = kadd(gt33rhsL,kmadd(epsdiss1,JacPDdissipationNth1gt33,kmadd(epsdiss2,JacPDdissipationNth2gt33,kmul(epsdiss3,JacPDdissipationNth3gt33)))); Xt1rhsL = kadd(Xt1rhsL,kmadd(epsdiss1,JacPDdissipationNth1Xt1,kmadd(epsdiss2,JacPDdissipationNth2Xt1,kmul(epsdiss3,JacPDdissipationNth3Xt1)))); Xt2rhsL = kadd(Xt2rhsL,kmadd(epsdiss1,JacPDdissipationNth1Xt2,kmadd(epsdiss2,JacPDdissipationNth2Xt2,kmul(epsdiss3,JacPDdissipationNth3Xt2)))); Xt3rhsL = kadd(Xt3rhsL,kmadd(epsdiss1,JacPDdissipationNth1Xt3,kmadd(epsdiss2,JacPDdissipationNth2Xt3,kmul(epsdiss3,JacPDdissipationNth3Xt3)))); ThetarhsL = kadd(ThetarhsL,kmadd(epsdiss1,JacPDdissipationNth1Theta,kmadd(epsdiss2,JacPDdissipationNth2Theta,kmul(epsdiss3,JacPDdissipationNth3Theta)))); trKrhsL = kadd(trKrhsL,kmadd(epsdiss1,JacPDdissipationNth1trK,kmadd(epsdiss2,JacPDdissipationNth2trK,kmul(epsdiss3,JacPDdissipationNth3trK)))); At11rhsL = kadd(At11rhsL,kmadd(epsdiss1,JacPDdissipationNth1At11,kmadd(epsdiss2,JacPDdissipationNth2At11,kmul(epsdiss3,JacPDdissipationNth3At11)))); At12rhsL = kadd(At12rhsL,kmadd(epsdiss1,JacPDdissipationNth1At12,kmadd(epsdiss2,JacPDdissipationNth2At12,kmul(epsdiss3,JacPDdissipationNth3At12)))); At13rhsL = kadd(At13rhsL,kmadd(epsdiss1,JacPDdissipationNth1At13,kmadd(epsdiss2,JacPDdissipationNth2At13,kmul(epsdiss3,JacPDdissipationNth3At13)))); At22rhsL = kadd(At22rhsL,kmadd(epsdiss1,JacPDdissipationNth1At22,kmadd(epsdiss2,JacPDdissipationNth2At22,kmul(epsdiss3,JacPDdissipationNth3At22)))); At23rhsL = kadd(At23rhsL,kmadd(epsdiss1,JacPDdissipationNth1At23,kmadd(epsdiss2,JacPDdissipationNth2At23,kmul(epsdiss3,JacPDdissipationNth3At23)))); At33rhsL = kadd(At33rhsL,kmadd(epsdiss1,JacPDdissipationNth1At33,kmadd(epsdiss2,JacPDdissipationNth2At33,kmul(epsdiss3,JacPDdissipationNth3At33)))); alpharhsL = kadd(alpharhsL,kmadd(epsdiss1,JacPDdissipationNth1alpha,kmadd(epsdiss2,JacPDdissipationNth2alpha,kmul(epsdiss3,JacPDdissipationNth3alpha)))); ArhsL = kadd(ArhsL,kmadd(epsdiss1,JacPDdissipationNth1A,kmadd(epsdiss2,JacPDdissipationNth2A,kmul(epsdiss3,JacPDdissipationNth3A)))); beta1rhsL = kadd(beta1rhsL,kmadd(epsdiss1,JacPDdissipationNth1beta1,kmadd(epsdiss2,JacPDdissipationNth2beta1,kmul(epsdiss3,JacPDdissipationNth3beta1)))); beta2rhsL = kadd(beta2rhsL,kmadd(epsdiss1,JacPDdissipationNth1beta2,kmadd(epsdiss2,JacPDdissipationNth2beta2,kmul(epsdiss3,JacPDdissipationNth3beta2)))); beta3rhsL = kadd(beta3rhsL,kmadd(epsdiss1,JacPDdissipationNth1beta3,kmadd(epsdiss2,JacPDdissipationNth2beta3,kmul(epsdiss3,JacPDdissipationNth3beta3)))); B1rhsL = kadd(B1rhsL,kmadd(epsdiss1,JacPDdissipationNth1B1,kmadd(epsdiss2,JacPDdissipationNth2B1,kmul(epsdiss3,JacPDdissipationNth3B1)))); B2rhsL = kadd(B2rhsL,kmadd(epsdiss1,JacPDdissipationNth1B2,kmadd(epsdiss2,JacPDdissipationNth2B2,kmul(epsdiss3,JacPDdissipationNth3B2)))); B3rhsL = kadd(B3rhsL,kmadd(epsdiss1,JacPDdissipationNth1B3,kmadd(epsdiss2,JacPDdissipationNth2B3,kmul(epsdiss3,JacPDdissipationNth3B3)))); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); vec_store_nta_partial(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); vec_store_nta_partial(At12rhs[index],At12rhsL); vec_store_nta_partial(At13rhs[index],At13rhsL); vec_store_nta_partial(At22rhs[index],At22rhsL); vec_store_nta_partial(At23rhs[index],At23rhsL); vec_store_nta_partial(At33rhs[index],At33rhsL); vec_store_nta_partial(B1rhs[index],B1rhsL); vec_store_nta_partial(B2rhs[index],B2rhsL); vec_store_nta_partial(B3rhs[index],B3rhsL); vec_store_nta_partial(beta1rhs[index],beta1rhsL); vec_store_nta_partial(beta2rhs[index],beta2rhsL); vec_store_nta_partial(beta3rhs[index],beta3rhsL); vec_store_nta_partial(gt11rhs[index],gt11rhsL); vec_store_nta_partial(gt12rhs[index],gt12rhsL); vec_store_nta_partial(gt13rhs[index],gt13rhsL); vec_store_nta_partial(gt22rhs[index],gt22rhsL); vec_store_nta_partial(gt23rhs[index],gt23rhsL); vec_store_nta_partial(gt33rhs[index],gt33rhsL); vec_store_nta_partial(phirhs[index],phirhsL); vec_store_nta_partial(Thetarhs[index],ThetarhsL); vec_store_nta_partial(trKrhs[index],trKrhsL); vec_store_nta_partial(Xt1rhs[index],Xt1rhsL); vec_store_nta_partial(Xt2rhs[index],Xt2rhsL); vec_store_nta_partial(Xt3rhs[index],Xt3rhsL); } CCTK_ENDLOOP3STR(ML_CCZ4_Dissipation); } extern "C" void ML_CCZ4_Dissipation(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; if (verbose > 1) { CCTK_VInfo(CCTK_THORNSTRING,"Entering ML_CCZ4_Dissipation_Body"); } if (cctk_iteration % ML_CCZ4_Dissipation_calc_every != ML_CCZ4_Dissipation_calc_offset) { return; } const char* const groups[] = { "ML_CCZ4::ML_curv", "ML_CCZ4::ML_curvrhs", "ML_CCZ4::ML_dtlapse", "ML_CCZ4::ML_dtlapserhs", "ML_CCZ4::ML_dtshift", "ML_CCZ4::ML_dtshiftrhs", "ML_CCZ4::ML_Gamma", "ML_CCZ4::ML_Gammarhs", "ML_CCZ4::ML_lapse", "ML_CCZ4::ML_lapserhs", "ML_CCZ4::ML_log_confac", "ML_CCZ4::ML_log_confacrhs", "ML_CCZ4::ML_metric", "ML_CCZ4::ML_metricrhs", "ML_CCZ4::ML_shift", "ML_CCZ4::ML_shiftrhs", "ML_CCZ4::ML_Theta", "ML_CCZ4::ML_Thetarhs", "ML_CCZ4::ML_trace_curv", "ML_CCZ4::ML_trace_curvrhs"}; GenericFD_AssertGroupStorage(cctkGH, "ML_CCZ4_Dissipation", 20, groups); switch (fdOrder) { case 2: { GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_Dissipation", 2, 2, 2); break; } case 4: { GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_Dissipation", 3, 3, 3); break; } case 6: { GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_Dissipation", 4, 4, 4); break; } case 8: { GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_Dissipation", 5, 5, 5); break; } default: CCTK_BUILTIN_UNREACHABLE(); } GenericFD_LoopOverInterior(cctkGH, ML_CCZ4_Dissipation_Body); if (verbose > 1) { CCTK_VInfo(CCTK_THORNSTRING,"Leaving ML_CCZ4_Dissipation_Body"); } }