aboutsummaryrefslogtreecommitdiff
path: root/ML_CCZ4/src/ML_CCZ4_Minkowski.cc
diff options
context:
space:
mode:
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_Minkowski.cc')
-rw-r--r--ML_CCZ4/src/ML_CCZ4_Minkowski.cc164
1 files changed, 84 insertions, 80 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_Minkowski.cc b/ML_CCZ4/src/ML_CCZ4_Minkowski.cc
index 77ae232..6c1e349 100644
--- a/ML_CCZ4/src/ML_CCZ4_Minkowski.cc
+++ b/ML_CCZ4/src/ML_CCZ4_Minkowski.cc
@@ -18,10 +18,14 @@
/* Define macros used in calculations */
#define INITVALUE (42)
-#define QAD(x) (SQR(SQR(x)))
+#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)))
static void ML_CCZ4_Minkowski_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[])
{
@@ -56,84 +60,84 @@ static void ML_CCZ4_Minkowski_Body(cGH const * restrict const cctkGH, int const
CCTK_REAL_VEC const hdzi = kmul(ToReal(0.5), dzi);
/* Initialize predefined quantities */
- CCTK_REAL_VEC const p1o1024dx = kmul(INV(dx),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o1024dy = kmul(INV(dy),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o1024dz = kmul(INV(dz),ToReal(0.0009765625));
- CCTK_REAL_VEC const p1o120dx = kmul(INV(dx),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o120dy = kmul(INV(dy),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o120dz = kmul(INV(dz),ToReal(0.00833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dx = kmul(INV(dx),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dy = kmul(INV(dy),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o12dz = kmul(INV(dz),ToReal(0.0833333333333333333333333333333));
- CCTK_REAL_VEC const p1o144dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o144dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o144dydz = kmul(INV(kmul(dy,dz)),ToReal(0.00694444444444444444444444444444));
- CCTK_REAL_VEC const p1o1680dx = kmul(INV(dx),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o1680dy = kmul(INV(dy),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o1680dz = kmul(INV(dz),ToReal(0.000595238095238095238095238095238));
- CCTK_REAL_VEC const p1o16dx = kmul(INV(dx),ToReal(0.0625));
- CCTK_REAL_VEC const p1o16dy = kmul(INV(dy),ToReal(0.0625));
- CCTK_REAL_VEC const p1o16dz = kmul(INV(dz),ToReal(0.0625));
- CCTK_REAL_VEC const p1o180dx2 = kmul(INV(SQR(dx)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o180dy2 = kmul(INV(SQR(dy)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o180dz2 = kmul(INV(SQR(dz)),ToReal(0.00555555555555555555555555555556));
- CCTK_REAL_VEC const p1o24dx = kmul(INV(dx),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o24dy = kmul(INV(dy),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o24dz = kmul(INV(dz),ToReal(0.0416666666666666666666666666667));
- CCTK_REAL_VEC const p1o256dx = kmul(INV(dx),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o256dy = kmul(INV(dy),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o256dz = kmul(INV(dz),ToReal(0.00390625));
- CCTK_REAL_VEC const p1o2dx = kmul(INV(dx),ToReal(0.5));
- CCTK_REAL_VEC const p1o2dy = kmul(INV(dy),ToReal(0.5));
- CCTK_REAL_VEC const p1o2dz = kmul(INV(dz),ToReal(0.5));
- CCTK_REAL_VEC const p1o3600dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o3600dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o3600dydz = kmul(INV(kmul(dy,dz)),ToReal(0.000277777777777777777777777777778));
- CCTK_REAL_VEC const p1o4dx = kmul(INV(dx),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dxdy = kmul(INV(kmul(dx,dy)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dxdz = kmul(INV(kmul(dx,dz)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dy = kmul(INV(dy),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dydz = kmul(INV(kmul(dy,dz)),ToReal(0.25));
- CCTK_REAL_VEC const p1o4dz = kmul(INV(dz),ToReal(0.25));
- CCTK_REAL_VEC const p1o5040dx2 = kmul(INV(SQR(dx)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o5040dy2 = kmul(INV(SQR(dy)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o5040dz2 = kmul(INV(SQR(dz)),ToReal(0.000198412698412698412698412698413));
- CCTK_REAL_VEC const p1o560dx = kmul(INV(dx),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o560dy = kmul(INV(dy),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o560dz = kmul(INV(dz),ToReal(0.00178571428571428571428571428571));
- CCTK_REAL_VEC const p1o60dx = kmul(INV(dx),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o60dy = kmul(INV(dy),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o60dz = kmul(INV(dz),ToReal(0.0166666666666666666666666666667));
- CCTK_REAL_VEC const p1o64dx = kmul(INV(dx),ToReal(0.015625));
- CCTK_REAL_VEC const p1o64dy = kmul(INV(dy),ToReal(0.015625));
- CCTK_REAL_VEC const p1o64dz = kmul(INV(dz),ToReal(0.015625));
- CCTK_REAL_VEC const p1o705600dxdy = kmul(INV(kmul(dx,dy)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o705600dxdz = kmul(INV(kmul(dx,dz)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o705600dydz = kmul(INV(kmul(dy,dz)),ToReal(1.41723356009070294784580498866e-6));
- CCTK_REAL_VEC const p1o840dx = kmul(INV(dx),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1o840dy = kmul(INV(dy),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1o840dz = kmul(INV(dz),ToReal(0.00119047619047619047619047619048));
- CCTK_REAL_VEC const p1odx = INV(dx);
- CCTK_REAL_VEC const p1odx2 = INV(SQR(dx));
- CCTK_REAL_VEC const p1ody = INV(dy);
- CCTK_REAL_VEC const p1ody2 = INV(SQR(dy));
- CCTK_REAL_VEC const p1odz = INV(dz);
- CCTK_REAL_VEC const p1odz2 = INV(SQR(dz));
- CCTK_REAL_VEC const pm1o120dx = kmul(INV(dx),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o120dy = kmul(INV(dy),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o120dz = kmul(INV(dz),ToReal(-0.00833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dx2 = kmul(INV(SQR(dx)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dy2 = kmul(INV(SQR(dy)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o12dz2 = kmul(INV(SQR(dz)),ToReal(-0.0833333333333333333333333333333));
- CCTK_REAL_VEC const pm1o2dx = kmul(INV(dx),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o2dy = kmul(INV(dy),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o2dz = kmul(INV(dz),ToReal(-0.5));
- CCTK_REAL_VEC const pm1o4dx = kmul(INV(dx),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o4dy = kmul(INV(dy),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o4dz = kmul(INV(dz),ToReal(-0.25));
- CCTK_REAL_VEC const pm1o60dx = kmul(INV(dx),ToReal(-0.0166666666666666666666666666667));
- CCTK_REAL_VEC const pm1o60dy = kmul(INV(dy),ToReal(-0.0166666666666666666666666666667));
- CCTK_REAL_VEC const pm1o60dz = kmul(INV(dz),ToReal(-0.0166666666666666666666666666667));
+ CCTK_REAL_VEC const p1o1024dx = kdiv(ToReal(0.0009765625),dx);
+ CCTK_REAL_VEC const p1o1024dy = kdiv(ToReal(0.0009765625),dy);
+ CCTK_REAL_VEC const p1o1024dz = kdiv(ToReal(0.0009765625),dz);
+ CCTK_REAL_VEC const p1o120dx = kdiv(ToReal(0.00833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const p1o120dy = kdiv(ToReal(0.00833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const p1o120dz = kdiv(ToReal(0.00833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const p1o12dx = kdiv(ToReal(0.0833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const p1o12dy = kdiv(ToReal(0.0833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const p1o12dz = kdiv(ToReal(0.0833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const p1o144dxdy = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o144dxdz = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o144dydz = kdiv(ToReal(0.00694444444444444444444444444444),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o1680dx = kdiv(ToReal(0.000595238095238095238095238095238),dx);
+ CCTK_REAL_VEC const p1o1680dy = kdiv(ToReal(0.000595238095238095238095238095238),dy);
+ CCTK_REAL_VEC const p1o1680dz = kdiv(ToReal(0.000595238095238095238095238095238),dz);
+ CCTK_REAL_VEC const p1o16dx = kdiv(ToReal(0.0625),dx);
+ CCTK_REAL_VEC const p1o16dy = kdiv(ToReal(0.0625),dy);
+ CCTK_REAL_VEC const p1o16dz = kdiv(ToReal(0.0625),dz);
+ CCTK_REAL_VEC const p1o180dx2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dx,dx));
+ CCTK_REAL_VEC const p1o180dy2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dy,dy));
+ CCTK_REAL_VEC const p1o180dz2 = kdiv(ToReal(0.00555555555555555555555555555556),kmul(dz,dz));
+ CCTK_REAL_VEC const p1o24dx = kdiv(ToReal(0.0416666666666666666666666666667),dx);
+ CCTK_REAL_VEC const p1o24dy = kdiv(ToReal(0.0416666666666666666666666666667),dy);
+ CCTK_REAL_VEC const p1o24dz = kdiv(ToReal(0.0416666666666666666666666666667),dz);
+ CCTK_REAL_VEC const p1o256dx = kdiv(ToReal(0.00390625),dx);
+ CCTK_REAL_VEC const p1o256dy = kdiv(ToReal(0.00390625),dy);
+ CCTK_REAL_VEC const p1o256dz = kdiv(ToReal(0.00390625),dz);
+ CCTK_REAL_VEC const p1o2dx = kdiv(ToReal(0.5),dx);
+ CCTK_REAL_VEC const p1o2dy = kdiv(ToReal(0.5),dy);
+ CCTK_REAL_VEC const p1o2dz = kdiv(ToReal(0.5),dz);
+ CCTK_REAL_VEC const p1o3600dxdy = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o3600dxdz = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o3600dydz = kdiv(ToReal(0.000277777777777777777777777777778),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o4dx = kdiv(ToReal(0.25),dx);
+ CCTK_REAL_VEC const p1o4dxdy = kdiv(ToReal(0.25),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o4dxdz = kdiv(ToReal(0.25),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o4dy = kdiv(ToReal(0.25),dy);
+ CCTK_REAL_VEC const p1o4dydz = kdiv(ToReal(0.25),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o4dz = kdiv(ToReal(0.25),dz);
+ CCTK_REAL_VEC const p1o5040dx2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dx,dx));
+ CCTK_REAL_VEC const p1o5040dy2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dy,dy));
+ CCTK_REAL_VEC const p1o5040dz2 = kdiv(ToReal(0.000198412698412698412698412698413),kmul(dz,dz));
+ CCTK_REAL_VEC const p1o560dx = kdiv(ToReal(0.00178571428571428571428571428571),dx);
+ CCTK_REAL_VEC const p1o560dy = kdiv(ToReal(0.00178571428571428571428571428571),dy);
+ CCTK_REAL_VEC const p1o560dz = kdiv(ToReal(0.00178571428571428571428571428571),dz);
+ CCTK_REAL_VEC const p1o60dx = kdiv(ToReal(0.0166666666666666666666666666667),dx);
+ CCTK_REAL_VEC const p1o60dy = kdiv(ToReal(0.0166666666666666666666666666667),dy);
+ CCTK_REAL_VEC const p1o60dz = kdiv(ToReal(0.0166666666666666666666666666667),dz);
+ CCTK_REAL_VEC const p1o64dx = kdiv(ToReal(0.015625),dx);
+ CCTK_REAL_VEC const p1o64dy = kdiv(ToReal(0.015625),dy);
+ CCTK_REAL_VEC const p1o64dz = kdiv(ToReal(0.015625),dz);
+ CCTK_REAL_VEC const p1o705600dxdy = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dy,dx));
+ CCTK_REAL_VEC const p1o705600dxdz = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dx));
+ CCTK_REAL_VEC const p1o705600dydz = kdiv(ToReal(1.41723356009070294784580498866e-6),kmul(dz,dy));
+ CCTK_REAL_VEC const p1o840dx = kdiv(ToReal(0.00119047619047619047619047619048),dx);
+ CCTK_REAL_VEC const p1o840dy = kdiv(ToReal(0.00119047619047619047619047619048),dy);
+ CCTK_REAL_VEC const p1o840dz = kdiv(ToReal(0.00119047619047619047619047619048),dz);
+ CCTK_REAL_VEC const p1odx = kdiv(ToReal(1),dx);
+ CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx));
+ CCTK_REAL_VEC const p1ody = kdiv(ToReal(1),dy);
+ CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy));
+ CCTK_REAL_VEC const p1odz = kdiv(ToReal(1),dz);
+ CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz));
+ CCTK_REAL_VEC const pm1o120dx = kdiv(ToReal(-0.00833333333333333333333333333333),dx);
+ CCTK_REAL_VEC const pm1o120dy = kdiv(ToReal(-0.00833333333333333333333333333333),dy);
+ CCTK_REAL_VEC const pm1o120dz = kdiv(ToReal(-0.00833333333333333333333333333333),dz);
+ CCTK_REAL_VEC const pm1o12dx2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dx,dx));
+ CCTK_REAL_VEC const pm1o12dy2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dy,dy));
+ CCTK_REAL_VEC const pm1o12dz2 = kdiv(ToReal(-0.0833333333333333333333333333333),kmul(dz,dz));
+ CCTK_REAL_VEC const pm1o2dx = kdiv(ToReal(-0.5),dx);
+ CCTK_REAL_VEC const pm1o2dy = kdiv(ToReal(-0.5),dy);
+ CCTK_REAL_VEC const pm1o2dz = kdiv(ToReal(-0.5),dz);
+ CCTK_REAL_VEC const pm1o4dx = kdiv(ToReal(-0.25),dx);
+ CCTK_REAL_VEC const pm1o4dy = kdiv(ToReal(-0.25),dy);
+ CCTK_REAL_VEC const pm1o4dz = kdiv(ToReal(-0.25),dz);
+ CCTK_REAL_VEC const pm1o60dx = kdiv(ToReal(-0.0166666666666666666666666666667),dx);
+ CCTK_REAL_VEC const pm1o60dy = kdiv(ToReal(-0.0166666666666666666666666666667),dy);
+ CCTK_REAL_VEC const pm1o60dz = kdiv(ToReal(-0.0166666666666666666666666666667),dz);
/* Jacobian variable pointers */
bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map)
@@ -192,7 +196,7 @@ static void ML_CCZ4_Minkowski_Body(cGH const * restrict const cctkGH, int const
#pragma omp parallel
LC_LOOP3VEC(ML_CCZ4_Minkowski,
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
- cctk_lsh[0],cctk_lsh[1],cctk_lsh[2],
+ cctk_ash[0],cctk_ash[1],cctk_ash[2],
CCTK_REAL_VEC_SIZE)
{
ptrdiff_t const index = di*i + dj*j + dk*k;
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-29 07:17:50 +0000