aboutsummaryrefslogtreecommitdiff
path: root/ML_ADMConstraints
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@gmail.com>2013-01-17 15:58:33 -0500
committerErik Schnetter <schnetter@gmail.com>2013-01-17 15:58:33 -0500
commitb457ddc6338d15f4a9c8109619f39a88698ebbf7 (patch)
treece3d973bbaeb6fcbfb675df473fad2494dcf191e /ML_ADMConstraints
parent76d692d840e042006df4b2355ecdf506f6b3d1d6 (diff)
Regenerate code after Kranc update
Diffstat (limited to 'ML_ADMConstraints')
-rw-r--r--ML_ADMConstraints/src/Differencing.h72
-rw-r--r--ML_ADMConstraints/src/ML_ADMConstraints.cc683
2 files changed, 375 insertions, 380 deletions
diff --git a/ML_ADMConstraints/src/Differencing.h b/ML_ADMConstraints/src/Differencing.h
index b27ce50..35539bd 100644
--- a/ML_ADMConstraints/src/Differencing.h
+++ b/ML_ADMConstraints/src/Differencing.h
@@ -1,144 +1,144 @@
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth1(u) ((-8*KRANC_GFOFFSET3D(u,-1,0,0) + 8*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx)
+# define PDstandardNth1(u) ((-8.*KRANC_GFOFFSET3D(u,-1,0,0) + 8.*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - 1.*KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx)
#else
# define PDstandardNth1(u) (PDstandardNth1_impl(u,p1o12dx,cdj,cdk))
static CCTK_REAL PDstandardNth1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dx, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-8*KRANC_GFOFFSET3D(u,-1,0,0) + 8*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-8.*KRANC_GFOFFSET3D(u,-1,0,0) + 8.*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - 1.*KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth2(u) ((-8*KRANC_GFOFFSET3D(u,0,-1,0) + 8*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy)
+# define PDstandardNth2(u) ((-8.*KRANC_GFOFFSET3D(u,0,-1,0) + 8.*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - 1.*KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy)
#else
# define PDstandardNth2(u) (PDstandardNth2_impl(u,p1o12dy,cdj,cdk))
static CCTK_REAL PDstandardNth2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dy, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-8*KRANC_GFOFFSET3D(u,0,-1,0) + 8*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-8.*KRANC_GFOFFSET3D(u,0,-1,0) + 8.*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - 1.*KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth3(u) ((-8*KRANC_GFOFFSET3D(u,0,0,-1) + 8*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz)
+# define PDstandardNth3(u) ((-8.*KRANC_GFOFFSET3D(u,0,0,-1) + 8.*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - 1.*KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz)
#else
# define PDstandardNth3(u) (PDstandardNth3_impl(u,p1o12dz,cdj,cdk))
static CCTK_REAL PDstandardNth3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-8*KRANC_GFOFFSET3D(u,0,0,-1) + 8*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-8.*KRANC_GFOFFSET3D(u,0,0,-1) + 8.*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - 1.*KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth11(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2)
+# define PDstandardNth11(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2)
#else
# define PDstandardNth11(u) (PDstandardNth11_impl(u,pm1o12dx2,cdj,cdk))
static CCTK_REAL PDstandardNth11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dx2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dx2, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth22(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2)
+# define PDstandardNth22(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2)
#else
# define PDstandardNth22(u) (PDstandardNth22_impl(u,pm1o12dy2,cdj,cdk))
static CCTK_REAL PDstandardNth22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dy2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dy2, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth33(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2)
+# define PDstandardNth33(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2)
#else
# define PDstandardNth33(u) (PDstandardNth33_impl(u,pm1o12dz2,cdj,cdk))
static CCTK_REAL PDstandardNth33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dz2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dz2, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth12(u) ((-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy)
+# define PDstandardNth12(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy)
#else
# define PDstandardNth12(u) (PDstandardNth12_impl(u,p1o144dxdy,cdj,cdk))
static CCTK_REAL PDstandardNth12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth13(u) ((-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz)
+# define PDstandardNth13(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz)
#else
# define PDstandardNth13(u) (PDstandardNth13_impl(u,p1o144dxdz,cdj,cdk))
static CCTK_REAL PDstandardNth13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth21(u) ((-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy)
+# define PDstandardNth21(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy)
#else
# define PDstandardNth21(u) (PDstandardNth21_impl(u,p1o144dxdy,cdj,cdk))
static CCTK_REAL PDstandardNth21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth23(u) ((-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz)
+# define PDstandardNth23(u) ((-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz)
#else
# define PDstandardNth23(u) (PDstandardNth23_impl(u,p1o144dydz,cdj,cdk))
static CCTK_REAL PDstandardNth23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth31(u) ((-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz)
+# define PDstandardNth31(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz)
#else
# define PDstandardNth31(u) (PDstandardNth31_impl(u,p1o144dxdz,cdj,cdk))
static CCTK_REAL PDstandardNth31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDstandardNth32(u) ((-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz)
+# define PDstandardNth32(u) ((-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz)
#else
# define PDstandardNth32(u) (PDstandardNth32_impl(u,p1o144dydz,cdj,cdk))
static CCTK_REAL PDstandardNth32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDstandardNth32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz;
}
#endif
diff --git a/ML_ADMConstraints/src/ML_ADMConstraints.cc b/ML_ADMConstraints/src/ML_ADMConstraints.cc
index a9f6db1..0b32fb1 100644
--- a/ML_ADMConstraints/src/ML_ADMConstraints.cc
+++ b/ML_ADMConstraints/src/ML_ADMConstraints.cc
@@ -46,42 +46,42 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
/* 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 const dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
- CCTK_REAL const dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
- CCTK_REAL const dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
- CCTK_REAL const dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
- CCTK_REAL const khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
- CCTK_REAL const kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
- CCTK_REAL const ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
- CCTK_REAL const kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
- CCTK_REAL const keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
- CCTK_REAL const hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
- CCTK_REAL const hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
- CCTK_REAL const hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx);
- CCTK_REAL const p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy);
- CCTK_REAL const p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz);
- CCTK_REAL const p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy);
- CCTK_REAL const p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz);
- CCTK_REAL const p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz);
- CCTK_REAL const pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx));
- CCTK_REAL const pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy));
- CCTK_REAL const pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz));
+ CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx);
+ CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy);
+ CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz);
+ CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy);
+ CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz);
+ CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz);
+ CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx));
+ CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy));
+ CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz));
/* Jacobian variable pointers */
- bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map)
+ bool const use_jacobian CCTK_ATTRIBUTE_UNUSED = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map)
&& strlen(jacobian_group) > 0;
if (use_jacobian && strlen(jacobian_derivative_group) == 0)
{
@@ -139,7 +139,7 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
@@ -225,78 +225,78 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
/* Include user supplied include files */
/* Precompute derivatives */
- CCTK_REAL const PDstandardNth1gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxx[index]);
- CCTK_REAL const PDstandardNth2gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxx[index]);
- CCTK_REAL const PDstandardNth3gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxx[index]);
- CCTK_REAL const PDstandardNth11gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxx[index]);
- CCTK_REAL const PDstandardNth22gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxx[index]);
- CCTK_REAL const PDstandardNth33gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxx[index]);
- CCTK_REAL const PDstandardNth12gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxx[index]);
- CCTK_REAL const PDstandardNth13gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxx[index]);
- CCTK_REAL const PDstandardNth23gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxx[index]);
- CCTK_REAL const PDstandardNth1gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxy[index]);
- CCTK_REAL const PDstandardNth2gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxy[index]);
- CCTK_REAL const PDstandardNth3gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxy[index]);
- CCTK_REAL const PDstandardNth11gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxy[index]);
- CCTK_REAL const PDstandardNth22gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxy[index]);
- CCTK_REAL const PDstandardNth33gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxy[index]);
- CCTK_REAL const PDstandardNth12gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxy[index]);
- CCTK_REAL const PDstandardNth13gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxy[index]);
- CCTK_REAL const PDstandardNth23gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxy[index]);
- CCTK_REAL const PDstandardNth1gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxz[index]);
- CCTK_REAL const PDstandardNth2gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxz[index]);
- CCTK_REAL const PDstandardNth3gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxz[index]);
- CCTK_REAL const PDstandardNth11gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxz[index]);
- CCTK_REAL const PDstandardNth22gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxz[index]);
- CCTK_REAL const PDstandardNth33gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxz[index]);
- CCTK_REAL const PDstandardNth12gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxz[index]);
- CCTK_REAL const PDstandardNth13gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxz[index]);
- CCTK_REAL const PDstandardNth23gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxz[index]);
- CCTK_REAL const PDstandardNth1gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gyy[index]);
- CCTK_REAL const PDstandardNth2gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gyy[index]);
- CCTK_REAL const PDstandardNth3gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gyy[index]);
- CCTK_REAL const PDstandardNth11gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gyy[index]);
- CCTK_REAL const PDstandardNth22gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gyy[index]);
- CCTK_REAL const PDstandardNth33gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gyy[index]);
- CCTK_REAL const PDstandardNth12gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gyy[index]);
- CCTK_REAL const PDstandardNth13gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gyy[index]);
- CCTK_REAL const PDstandardNth23gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gyy[index]);
- CCTK_REAL const PDstandardNth1gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gyz[index]);
- CCTK_REAL const PDstandardNth2gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gyz[index]);
- CCTK_REAL const PDstandardNth3gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gyz[index]);
- CCTK_REAL const PDstandardNth11gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gyz[index]);
- CCTK_REAL const PDstandardNth22gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gyz[index]);
- CCTK_REAL const PDstandardNth33gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gyz[index]);
- CCTK_REAL const PDstandardNth12gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gyz[index]);
- CCTK_REAL const PDstandardNth13gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gyz[index]);
- CCTK_REAL const PDstandardNth23gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gyz[index]);
- CCTK_REAL const PDstandardNth1gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gzz[index]);
- CCTK_REAL const PDstandardNth2gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gzz[index]);
- CCTK_REAL const PDstandardNth3gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gzz[index]);
- CCTK_REAL const PDstandardNth11gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gzz[index]);
- CCTK_REAL const PDstandardNth22gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gzz[index]);
- CCTK_REAL const PDstandardNth33gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gzz[index]);
- CCTK_REAL const PDstandardNth12gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gzz[index]);
- CCTK_REAL const PDstandardNth13gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gzz[index]);
- CCTK_REAL const PDstandardNth23gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gzz[index]);
- CCTK_REAL const PDstandardNth1kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxx[index]);
- CCTK_REAL const PDstandardNth2kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxx[index]);
- CCTK_REAL const PDstandardNth3kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxx[index]);
- CCTK_REAL const PDstandardNth1kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxy[index]);
- CCTK_REAL const PDstandardNth2kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxy[index]);
- CCTK_REAL const PDstandardNth3kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxy[index]);
- CCTK_REAL const PDstandardNth1kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxz[index]);
- CCTK_REAL const PDstandardNth2kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxz[index]);
- CCTK_REAL const PDstandardNth3kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxz[index]);
- CCTK_REAL const PDstandardNth1kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kyy[index]);
- CCTK_REAL const PDstandardNth2kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kyy[index]);
- CCTK_REAL const PDstandardNth3kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kyy[index]);
- CCTK_REAL const PDstandardNth1kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kyz[index]);
- CCTK_REAL const PDstandardNth2kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kyz[index]);
- CCTK_REAL const PDstandardNth3kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kyz[index]);
- CCTK_REAL const PDstandardNth1kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kzz[index]);
- CCTK_REAL const PDstandardNth2kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kzz[index]);
- CCTK_REAL const PDstandardNth3kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth11gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth22gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth33gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth12gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth13gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth23gzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(&gzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kxx CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxx[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kxy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxy[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kxz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kxz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kyy CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kyy[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kyz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kyz[index]);
+ CCTK_REAL /*const*/ PDstandardNth1kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&kzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth2kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&kzz[index]);
+ CCTK_REAL /*const*/ PDstandardNth3kzz CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&kzz[index]);
/* Calculate temporaries and grid functions */
CCTK_REAL JacPDstandardNth11gyy CCTK_ATTRIBUTE_UNUSED ;
@@ -480,55 +480,55 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
J23L*PDstandardNth2kyz + J33L*PDstandardNth3kyz;
JacPDstandardNth11gyy = dJ111L*PDstandardNth1gyy +
- 2*(J11L*(J21L*PDstandardNth12gyy + J31L*PDstandardNth13gyy) +
+ 2.*(J11L*(J21L*PDstandardNth12gyy + J31L*PDstandardNth13gyy) +
J21L*J31L*PDstandardNth23gyy) + dJ211L*PDstandardNth2gyy +
dJ311L*PDstandardNth3gyy + PDstandardNth11gyy*SQR(J11L) +
PDstandardNth22gyy*SQR(J21L) + PDstandardNth33gyy*SQR(J31L);
JacPDstandardNth11gyz = dJ111L*PDstandardNth1gyz +
- 2*(J11L*(J21L*PDstandardNth12gyz + J31L*PDstandardNth13gyz) +
+ 2.*(J11L*(J21L*PDstandardNth12gyz + J31L*PDstandardNth13gyz) +
J21L*J31L*PDstandardNth23gyz) + dJ211L*PDstandardNth2gyz +
dJ311L*PDstandardNth3gyz + PDstandardNth11gyz*SQR(J11L) +
PDstandardNth22gyz*SQR(J21L) + PDstandardNth33gyz*SQR(J31L);
JacPDstandardNth11gzz = dJ111L*PDstandardNth1gzz +
- 2*(J11L*(J21L*PDstandardNth12gzz + J31L*PDstandardNth13gzz) +
+ 2.*(J11L*(J21L*PDstandardNth12gzz + J31L*PDstandardNth13gzz) +
J21L*J31L*PDstandardNth23gzz) + dJ211L*PDstandardNth2gzz +
dJ311L*PDstandardNth3gzz + PDstandardNth11gzz*SQR(J11L) +
PDstandardNth22gzz*SQR(J21L) + PDstandardNth33gzz*SQR(J31L);
JacPDstandardNth22gxx = dJ122L*PDstandardNth1gxx +
- 2*(J12L*(J22L*PDstandardNth12gxx + J32L*PDstandardNth13gxx) +
+ 2.*(J12L*(J22L*PDstandardNth12gxx + J32L*PDstandardNth13gxx) +
J22L*J32L*PDstandardNth23gxx) + dJ222L*PDstandardNth2gxx +
dJ322L*PDstandardNth3gxx + PDstandardNth11gxx*SQR(J12L) +
PDstandardNth22gxx*SQR(J22L) + PDstandardNth33gxx*SQR(J32L);
JacPDstandardNth22gxz = dJ122L*PDstandardNth1gxz +
- 2*(J12L*(J22L*PDstandardNth12gxz + J32L*PDstandardNth13gxz) +
+ 2.*(J12L*(J22L*PDstandardNth12gxz + J32L*PDstandardNth13gxz) +
J22L*J32L*PDstandardNth23gxz) + dJ222L*PDstandardNth2gxz +
dJ322L*PDstandardNth3gxz + PDstandardNth11gxz*SQR(J12L) +
PDstandardNth22gxz*SQR(J22L) + PDstandardNth33gxz*SQR(J32L);
JacPDstandardNth22gzz = dJ122L*PDstandardNth1gzz +
- 2*(J12L*(J22L*PDstandardNth12gzz + J32L*PDstandardNth13gzz) +
+ 2.*(J12L*(J22L*PDstandardNth12gzz + J32L*PDstandardNth13gzz) +
J22L*J32L*PDstandardNth23gzz) + dJ222L*PDstandardNth2gzz +
dJ322L*PDstandardNth3gzz + PDstandardNth11gzz*SQR(J12L) +
PDstandardNth22gzz*SQR(J22L) + PDstandardNth33gzz*SQR(J32L);
JacPDstandardNth33gxx = dJ133L*PDstandardNth1gxx +
- 2*(J13L*(J23L*PDstandardNth12gxx + J33L*PDstandardNth13gxx) +
+ 2.*(J13L*(J23L*PDstandardNth12gxx + J33L*PDstandardNth13gxx) +
J23L*J33L*PDstandardNth23gxx) + dJ233L*PDstandardNth2gxx +
dJ333L*PDstandardNth3gxx + PDstandardNth11gxx*SQR(J13L) +
PDstandardNth22gxx*SQR(J23L) + PDstandardNth33gxx*SQR(J33L);
JacPDstandardNth33gxy = dJ133L*PDstandardNth1gxy +
- 2*(J13L*(J23L*PDstandardNth12gxy + J33L*PDstandardNth13gxy) +
+ 2.*(J13L*(J23L*PDstandardNth12gxy + J33L*PDstandardNth13gxy) +
J23L*J33L*PDstandardNth23gxy) + dJ233L*PDstandardNth2gxy +
dJ333L*PDstandardNth3gxy + PDstandardNth11gxy*SQR(J13L) +
PDstandardNth22gxy*SQR(J23L) + PDstandardNth33gxy*SQR(J33L);
JacPDstandardNth33gyy = dJ133L*PDstandardNth1gyy +
- 2*(J13L*(J23L*PDstandardNth12gyy + J33L*PDstandardNth13gyy) +
+ 2.*(J13L*(J23L*PDstandardNth12gyy + J33L*PDstandardNth13gyy) +
J23L*J33L*PDstandardNth23gyy) + dJ233L*PDstandardNth2gyy +
dJ333L*PDstandardNth3gyy + PDstandardNth11gyy*SQR(J13L) +
PDstandardNth22gyy*SQR(J23L) + PDstandardNth33gyy*SQR(J33L);
@@ -980,136 +980,135 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
JacPDstandardNth32gzz = PDstandardNth23gzz;
}
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED detg = 2*gxyL*gxzL*gyzL +
- gzzL*(gxxL*gyyL - SQR(gxyL)) - gyyL*SQR(gxzL) -
- gxxL*SQR(gyzL);
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED detg = 2.*gxyL*gxzL*gyzL +
+ gzzL*(gxxL*gyyL - 1.*SQR(gxyL)) - 1.*(gyyL*SQR(gxzL) +
+ gxxL*SQR(gyzL));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu11 = INV(detg)*(gyyL*gzzL -
- SQR(gyzL));
+ 1.*SQR(gyzL));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu21 = (gxzL*gyzL -
- gxyL*gzzL)*INV(detg);
+ 1.*gxyL*gzzL)*INV(detg);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu31 = (-(gxzL*gyyL) +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu31 = (-1.*gxzL*gyyL +
gxyL*gyzL)*INV(detg);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu22 = INV(detg)*(gxxL*gzzL -
- SQR(gxzL));
+ 1.*SQR(gxzL));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu32 = (gxyL*gxzL -
- gxxL*gyzL)*INV(detg);
+ 1.*gxxL*gyzL)*INV(detg);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu33 = INV(detg)*(gxxL*gyyL -
- SQR(gxyL));
+ 1.*SQR(gxyL));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G111 = 0.5*(gu11*JacPDstandardNth1gxx
- + 2*(gu21*JacPDstandardNth1gxy + gu31*JacPDstandardNth1gxz) -
- gu21*JacPDstandardNth2gxx - gu31*JacPDstandardNth3gxx);
+ + 2.*(gu21*JacPDstandardNth1gxy + gu31*JacPDstandardNth1gxz) -
+ 1.*(gu21*JacPDstandardNth2gxx + gu31*JacPDstandardNth3gxx));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G211 = 0.5*(gu21*JacPDstandardNth1gxx
- + 2*(gu22*JacPDstandardNth1gxy + gu32*JacPDstandardNth1gxz) -
- gu22*JacPDstandardNth2gxx - gu32*JacPDstandardNth3gxx);
+ + 2.*(gu22*JacPDstandardNth1gxy + gu32*JacPDstandardNth1gxz) -
+ 1.*(gu22*JacPDstandardNth2gxx + gu32*JacPDstandardNth3gxx));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G311 = 0.5*(gu31*JacPDstandardNth1gxx
- + 2*(gu32*JacPDstandardNth1gxy + gu33*JacPDstandardNth1gxz) -
- gu32*JacPDstandardNth2gxx - gu33*JacPDstandardNth3gxx);
+ + 2.*(gu32*JacPDstandardNth1gxy + gu33*JacPDstandardNth1gxz) -
+ 1.*(gu32*JacPDstandardNth2gxx + gu33*JacPDstandardNth3gxx));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G112 = 0.5*(gu21*JacPDstandardNth1gyy
+ gu11*JacPDstandardNth2gxx + gu31*(JacPDstandardNth1gyz +
- JacPDstandardNth2gxz - JacPDstandardNth3gxy));
+ JacPDstandardNth2gxz - 1.*JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G212 = 0.5*(gu22*JacPDstandardNth1gyy
+ gu21*JacPDstandardNth2gxx + gu32*(JacPDstandardNth1gyz +
- JacPDstandardNth2gxz - JacPDstandardNth3gxy));
+ JacPDstandardNth2gxz - 1.*JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G312 = 0.5*(gu32*JacPDstandardNth1gyy
+ gu31*JacPDstandardNth2gxx + gu33*(JacPDstandardNth1gyz +
- JacPDstandardNth2gxz - JacPDstandardNth3gxy));
+ JacPDstandardNth2gxz - 1.*JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G113 = 0.5*(gu31*JacPDstandardNth1gzz
+ gu11*JacPDstandardNth3gxx + gu21*(JacPDstandardNth1gyz -
- JacPDstandardNth2gxz + JacPDstandardNth3gxy));
+ 1.*JacPDstandardNth2gxz + JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G213 = 0.5*(gu32*JacPDstandardNth1gzz
+ gu21*JacPDstandardNth3gxx + gu22*(JacPDstandardNth1gyz -
- JacPDstandardNth2gxz + JacPDstandardNth3gxy));
+ 1.*JacPDstandardNth2gxz + JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G313 = 0.5*(gu33*JacPDstandardNth1gzz
+ gu31*JacPDstandardNth3gxx + gu32*(JacPDstandardNth1gyz -
- JacPDstandardNth2gxz + JacPDstandardNth3gxy));
+ 1.*JacPDstandardNth2gxz + JacPDstandardNth3gxy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G122 =
- 0.5*(gu11*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) +
- gu21*JacPDstandardNth2gyy + gu31*(2*JacPDstandardNth2gyz -
- JacPDstandardNth3gyy));
+ 0.5*(gu11*(-1.*JacPDstandardNth1gyy + 2.*JacPDstandardNth2gxy) +
+ gu21*JacPDstandardNth2gyy + gu31*(2.*JacPDstandardNth2gyz -
+ 1.*JacPDstandardNth3gyy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G222 =
- 0.5*(gu21*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) +
- gu22*JacPDstandardNth2gyy + gu32*(2*JacPDstandardNth2gyz -
- JacPDstandardNth3gyy));
+ 0.5*(gu21*(-1.*JacPDstandardNth1gyy + 2.*JacPDstandardNth2gxy) +
+ gu22*JacPDstandardNth2gyy + gu32*(2.*JacPDstandardNth2gyz -
+ 1.*JacPDstandardNth3gyy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G322 =
- 0.5*(gu31*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) +
- gu32*JacPDstandardNth2gyy + gu33*(2*JacPDstandardNth2gyz -
- JacPDstandardNth3gyy));
+ 0.5*(gu31*(-1.*JacPDstandardNth1gyy + 2.*JacPDstandardNth2gxy) +
+ gu32*JacPDstandardNth2gyy + gu33*(2.*JacPDstandardNth2gyz -
+ 1.*JacPDstandardNth3gyy));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G123 = 0.5*(gu31*JacPDstandardNth2gzz
- + gu11*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz +
+ + gu11*(-1.*JacPDstandardNth1gyz + JacPDstandardNth2gxz +
JacPDstandardNth3gxy) + gu21*JacPDstandardNth3gyy);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G223 = 0.5*(gu32*JacPDstandardNth2gzz
- + gu21*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz +
+ + gu21*(-1.*JacPDstandardNth1gyz + JacPDstandardNth2gxz +
JacPDstandardNth3gxy) + gu22*JacPDstandardNth3gyy);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G323 = 0.5*(gu33*JacPDstandardNth2gzz
- + gu31*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz +
+ + gu31*(-1.*JacPDstandardNth1gyz + JacPDstandardNth2gxz +
JacPDstandardNth3gxy) + gu32*JacPDstandardNth3gyy);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G133 =
- 0.5*(gu11*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) +
- gu21*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) +
+ 0.5*(-1.*(gu11*JacPDstandardNth1gzz + gu21*JacPDstandardNth2gzz) +
+ 2.*(gu11*JacPDstandardNth3gxz + gu21*JacPDstandardNth3gyz) +
gu31*JacPDstandardNth3gzz);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G233 =
- 0.5*(gu21*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) +
- gu22*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) +
+ 0.5*(-1.*(gu21*JacPDstandardNth1gzz + gu22*JacPDstandardNth2gzz) +
+ 2.*(gu21*JacPDstandardNth3gxz + gu22*JacPDstandardNth3gyz) +
gu32*JacPDstandardNth3gzz);
CCTK_REAL CCTK_ATTRIBUTE_UNUSED G333 =
- 0.5*(gu31*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) +
- gu32*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) +
+ 0.5*(-1.*(gu31*JacPDstandardNth1gzz + gu32*JacPDstandardNth2gzz) +
+ 2.*(gu31*JacPDstandardNth3gxz + gu32*JacPDstandardNth3gyz) +
gu33*JacPDstandardNth3gzz);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R11 = 0.5*(gxyL*(-4*G123*G211*gu32 +
- 4*G113*G212*gu32 + 4*G113*G213*gu33) + gxzL*(-4*G123*G311*gu32 +
- 4*G113*G312*gu32 - 2*G133*G311*gu33) + 4*(gyzL*G213*G312*gu32 +
- gyzL*G212*G313*gu32 + G112*(gxyL*G212*gu22 + gxzL*G312*gu22 +
- gxxL*G113*gu32 + gxyL*G213*gu32 + gxzL*G313*gu32) +
- gxzL*G113*G313*gu33) + gyzL*(4*G212*G312*gu22 - 2*G211*G322*gu22 -
- 4*G223*G311*gu32 - 2*G233*G311*gu33 + 4*G213*G313*gu33 -
- 2*G211*G333*gu33) - 2*(gxyL*G122*G211*gu22 + gyyL*G211*G222*gu22 +
- gxzL*G122*G311*gu22 + gyzL*G222*G311*gu22 + gzzL*G311*G322*gu22 +
- gxyL*G133*G211*gu33 + gzzL*G311*G333*gu33 + G111*(gxxL*G122*gu22
- + gxyL*G222*gu22 + gxzL*G322*gu22 + 2*gxxL*G123*gu32 +
- 2*gxyL*G223*gu32 + 2*gxzL*G323*gu32 + gxxL*G133*gu33 +
- gxyL*G233*gu33 + gxzL*G333*gu33)) + gu21*(-JacPDstandardNth12gxx +
- JacPDstandardNth21gxx) + gu22*(-JacPDstandardNth11gyy +
- 2*JacPDstandardNth21gxy - JacPDstandardNth22gxx) +
- gu32*(4*gzzL*G312*G313 - 4*gyzL*G211*G323 - 4*gzzL*G311*G323 -
- 2*JacPDstandardNth11gyz - JacPDstandardNth23gxx) +
- gu31*(-JacPDstandardNth13gxx + JacPDstandardNth31gxx) +
- 2*gu32*JacPDstandardNth31gxy + gu32*(2*JacPDstandardNth21gxz -
- JacPDstandardNth32gxx) + gu33*(-JacPDstandardNth11gzz +
- 2*JacPDstandardNth31gxz - JacPDstandardNth33gxx) +
- 2*gxxL*gu22*SQR(G112) + 2*gxxL*gu33*SQR(G113) +
- gyyL*(4*G212*G213*gu32 - 4*G211*G223*gu32 - 2*G211*G233*gu33 +
- 2*gu22*SQR(G212)) + 2*gyyL*gu33*SQR(G213) + 2*gzzL*gu22*SQR(G312) +
- 2*gzzL*gu33*SQR(G313));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R11 = 0.5*(gxyL*((-4.*G123*G211 +
+ 4.*G113*G212)*gu32 + 4.*G113*G213*gu33) + gxzL*((-4.*G123*G311 +
+ 4.*G113*G312)*gu32 - 2.*G133*G311*gu33) + 4.*(gyzL*(G213*G312 +
+ G212*G313)*gu32 + G112*((gxyL*G212 + gxzL*G312)*gu22 + (gxxL*G113
+ + gxyL*G213 + gxzL*G313)*gu32) + gxzL*G113*G313*gu33) +
+ gyzL*((4.*G212*G312 - 2.*G211*G322)*gu22 + (4.*G213*G313 -
+ 2.*G211*G333)*gu33 + G311*(-4.*G223*gu32 - 2.*G233*gu33)) -
+ 2.*(G211*((gxyL*G122 + gyyL*G222)*gu22 + gxyL*G133*gu33) +
+ G311*((gxzL*G122 + gyzL*G222 + gzzL*G322)*gu22 +
+ gzzL*G333*gu33) + G111*((gxxL*G122 + gxyL*G222 +
+ gxzL*G322)*gu22 + 2.*(gxxL*G123 + gxyL*G223 + gxzL*G323)*gu32 +
+ (gxxL*G133 + gxyL*G233 + gxzL*G333)*gu33)) +
+ gu21*JacPDstandardNth21gxx + gu32*(-4.*gyzL*G211*G323 +
+ gzzL*(4.*G312*G313 - 4.*G311*G323) - 2.*JacPDstandardNth11gyz +
+ 2.*JacPDstandardNth21gxz - 1.*JacPDstandardNth23gxx) +
+ gu31*JacPDstandardNth31gxx - 1.*(gu21*JacPDstandardNth12gxx +
+ gu31*JacPDstandardNth13gxx + gu22*(JacPDstandardNth11gyy +
+ JacPDstandardNth22gxx) + gu32*JacPDstandardNth32gxx +
+ gu33*(JacPDstandardNth11gzz + JacPDstandardNth33gxx)) +
+ gyyL*((4.*G212*G213 - 4.*G211*G223)*gu32 - 2.*G211*G233*gu33 +
+ 2.*gu22*SQR(G212)) + 2.*(gu32*JacPDstandardNth31gxy +
+ gu22*(JacPDstandardNth21gxy + gxxL*SQR(G112) + gzzL*SQR(G312)) +
+ gu33*(JacPDstandardNth31gxz + gxxL*SQR(G113) + gyyL*SQR(G213) +
+ gzzL*SQR(G313))));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R12 = 0.5*(gxzL*((-2*G123*G312 +
- 2*G122*G313)*gu32 - 2*G133*G312*gu33) + gyzL*((-4*G212*G312 +
- 2*G211*G322)*gu21 + G223*(2*G311*gu31 - 2*G312*gu32) -
- 2*G233*G312*gu33) + 2*(G111*((gxxL*G122 + gxyL*G222 +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R12 = 0.5*(gxzL*((-2.*G123*G312 +
+ 2.*G122*G313)*gu32 - 2.*G133*G312*gu33) + gyzL*((-4.*G212*G312 +
+ 2.*G211*G322)*gu21 + G223*(2.*G311*gu31 - 2.*G312*gu32) -
+ 2.*G233*G312*gu33) + 2.*(G111*((gxxL*G122 + gxyL*G222 +
gxzL*G322)*gu21 + (gxxL*G123 + gxyL*G223 + gxzL*G323)*gu31) +
G211*((gxyL*G122 + gyyL*G222)*gu21 + (gxyL*G123 + gyyL*G223 +
gyzL*G323)*gu31) + G311*((gxzL*G122 + gyzL*G222 +
@@ -1117,142 +1116,137 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
(G122*(gxxL*G113 + gxyL*G213) + G222*(gxyL*G113 + gyzL*G313) +
(gxzL*G113 + gyzL*G213 + gzzL*G313)*G322)*gu32 +
(G123*(gxxL*G113 + gxyL*G213) + gxyL*G113*G223 + (gyzL*G213 +
- gzzL*G313)*G323)*gu33) + gu22*(-JacPDstandardNth12gyy +
+ gzzL*G313)*G323)*gu33) + gu22*(-1.*JacPDstandardNth12gyy +
JacPDstandardNth21gyy) + gu32*(JacPDstandardNth21gyz +
- JacPDstandardNth22gxz - JacPDstandardNth23gxy + JacPDstandardNth31gyy)
- + gu31*(JacPDstandardNth11gyz - JacPDstandardNth12gxz -
- JacPDstandardNth13gxy + JacPDstandardNth32gxx) +
- gu33*(2*(gxzL*G123*G313 + gyzL*G223*G313 + gxzL*G113*G323) -
- 2*(gyyL*G212*G233 + gyzL*G212*G333 + gzzL*G312*G333) -
- JacPDstandardNth12gzz + JacPDstandardNth31gyz + JacPDstandardNth32gxz -
- JacPDstandardNth33gxy) - 2*((gyzL*G212 + gzzL*G312)*G313*gu31 +
- G312*((gxzL*G113 + gyzL*G213)*gu31 + gzzL*G323*gu32) +
- G212*((gxyL*G113 + gyyL*G213)*gu31 + (gxyL*G123 +
- gyzL*G323)*gu32 + gxyL*G133*gu33) + G112*(2*(gxyL*G212 +
- gxzL*G312)*gu21 + (gxxL*G113 + gxyL*G213 + gxzL*G313)*gu31 +
- (gxxL*G123 + gxyL*G223 + gxzL*G323)*gu32 + (gxxL*G133 +
- gxyL*G233 + gxzL*G333)*gu33) + gu32*JacPDstandardNth12gyz +
- gu21*(JacPDstandardNth12gxy + gxxL*SQR(G112))) + gyyL*((2*G213*G222
- - 2*G212*G223)*gu32 + 2*G213*G223*gu33 - 2*gu21*SQR(G212)) +
- gu21*(JacPDstandardNth11gyy + JacPDstandardNth22gxx -
- 2*gzzL*SQR(G312)));
+ JacPDstandardNth22gxz - 1.*JacPDstandardNth23gxy +
+ JacPDstandardNth31gyy) + gu31*(JacPDstandardNth11gyz -
+ 1.*(JacPDstandardNth12gxz + JacPDstandardNth13gxy) +
+ JacPDstandardNth32gxx) + gu33*(2.*((gxzL*G123 + gyzL*G223)*G313 +
+ gxzL*G113*G323) - 2.*(gzzL*G312*G333 + G212*(gyyL*G233 +
+ gyzL*G333)) + JacPDstandardNth31gyz + JacPDstandardNth32gxz -
+ 1.*(JacPDstandardNth12gzz + JacPDstandardNth33gxy)) - 2.*((gyzL*G212
+ + gzzL*G312)*G313*gu31 + G312*((gxzL*G113 + gyzL*G213)*gu31 +
+ gzzL*G323*gu32) + G212*((gxyL*G113 + gyyL*G213)*gu31 +
+ (gxyL*G123 + gyzL*G323)*gu32 + gxyL*G133*gu33) +
+ G112*(2.*(gxyL*G212 + gxzL*G312)*gu21 + (gxxL*G113 + gxyL*G213
+ + gxzL*G313)*gu31 + (gxxL*G123 + gxyL*G223 + gxzL*G323)*gu32 +
+ (gxxL*G133 + gxyL*G233 + gxzL*G333)*gu33) +
+ gu32*JacPDstandardNth12gyz + gu21*(JacPDstandardNth12gxy +
+ gxxL*SQR(G112))) + gyyL*((2.*G213*G222 - 2.*G212*G223)*gu32 +
+ 2.*G213*G223*gu33 - 2.*gu21*SQR(G212)) + gu21*(JacPDstandardNth11gyy +
+ JacPDstandardNth22gxx - 2.*gzzL*SQR(G312)));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R13 = 0.5*(gxyL*(2*G123*G211*gu21 -
- 2*G113*G212*gu21 + 2*G123*G212*gu22 + 2*G133*G211*gu31 +
- 2*G133*G212*gu32) + 2*(gxzL*G123*G311*gu21 + gyzL*G223*G311*gu21 +
- gyzL*G211*G323*gu21 + gzzL*G311*G323*gu21 + gyyL*G212*G223*gu22 +
- gxzL*G123*G312*gu22 + gyzL*G223*G312*gu22 + gxzL*G133*G311*gu31 +
- G111*(gxxL*G123*gu21 + gxyL*G223*gu21 + gxzL*G323*gu21 +
- gxxL*G133*gu31 + gxyL*G233*gu31 + gxzL*G333*gu31) +
- gyyL*G212*G233*gu32 + gxzL*G133*G312*gu32 + gyzL*G233*G312*gu32 +
- G112*((-(gxxL*G113) - gxyL*G213 - gxzL*G313)*gu21 +
- gxxL*G123*gu22 + gxyL*G223*gu22 + gxzL*G323*gu22 +
- gxxL*G133*gu32 + gxyL*G233*gu32 + gxzL*G333*gu32)) +
- gu22*(-2*gzzL*G313*G322 + 2*gyzL*G212*G323 + 2*gzzL*G312*G323 -
- JacPDstandardNth13gyy + JacPDstandardNth21gyz - JacPDstandardNth22gxz)
- + gu21*(JacPDstandardNth11gyz - JacPDstandardNth12gxz -
- JacPDstandardNth13gxy + JacPDstandardNth23gxx) +
- gu22*JacPDstandardNth23gxy + gu32*JacPDstandardNth31gyz +
- gu33*(-JacPDstandardNth13gzz + JacPDstandardNth31gzz) +
- gu32*(-2*gzzL*G313*G323 + 2*gyzL*G212*G333 + 2*gzzL*G312*G333 -
- 2*JacPDstandardNth13gyz + JacPDstandardNth21gzz -
- JacPDstandardNth32gxz) + gu31*(-4*gxyL*G113*G213 + 2*gyzL*G233*G311
- - 4*gxzL*G113*G313 - 4*gyzL*G213*G313 + 2*gyzL*G211*G333 +
- 2*gzzL*G311*G333 + JacPDstandardNth11gzz + JacPDstandardNth33gxx) +
- gu32*JacPDstandardNth33gxy - 2*(gxzL*G113*G312*gu21 +
- gyzL*G213*G312*gu21 + gyzL*G212*G313*gu21 + gzzL*G312*G313*gu21 +
- gxxL*G113*G122*gu22 + gxyL*G122*G213*gu22 + gxyL*G113*G222*gu22 +
- gxzL*G122*G313*gu22 + gyzL*G222*G313*gu22 + gxzL*G113*G322*gu22 +
- gyzL*G213*G322*gu22 + gxxL*G113*G123*gu32 + gxyL*G123*G213*gu32 +
- gxyL*G113*G223*gu32 + gxzL*G123*G313*gu32 + gyzL*G223*G313*gu32 +
- gxzL*G113*G323*gu32 + gyzL*G213*G323*gu32 +
- gu31*JacPDstandardNth13gxz + gxxL*gu31*SQR(G113)) +
- gyyL*(-2*G212*G213*gu21 + 2*G211*G223*gu21 - 2*G213*G222*gu22 +
- 2*G211*G233*gu31 - 2*G213*G223*gu32 - 2*gu31*SQR(G213)) -
- 2*gzzL*gu31*SQR(G313));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R13 = 0.5*(2.*((gyzL*G211 +
+ gzzL*G311)*G323*gu21 + (gxzL*G123*G312 + G223*(gyyL*G212 +
+ gyzL*G312))*gu22 + G311*((gxzL*G123 + gyzL*G223)*gu21 +
+ gxzL*G133*gu31) + G111*((gxxL*G123 + gxyL*G223 +
+ gxzL*G323)*gu21 + (gxxL*G133 + gxyL*G233 + gxzL*G333)*gu31) +
+ (gxzL*G133*G312 + G233*(gyyL*G212 + gyzL*G312))*gu32 +
+ G112*(-1.*(gxxL*G113 + gxyL*G213 + gxzL*G313)*gu21 + (gxxL*G123
+ + gxyL*G223 + gxzL*G323)*gu22 + (gxxL*G133 + gxyL*G233 +
+ gxzL*G333)*gu32)) + gxyL*((2.*G123*G211 - 2.*G113*G212)*gu21 +
+ 2.*(G133*G211*gu31 + G212*(G123*gu22 + G133*gu32))) +
+ gu21*(JacPDstandardNth11gyz - 1.*(JacPDstandardNth12gxz +
+ JacPDstandardNth13gxy) + JacPDstandardNth23gxx) +
+ gu22*(2.*gyzL*G212*G323 + gzzL*(-2.*G313*G322 + 2.*G312*G323) -
+ 1.*JacPDstandardNth13gyy + JacPDstandardNth21gyz +
+ JacPDstandardNth23gxy) + gu33*JacPDstandardNth31gzz -
+ 1.*(gu33*JacPDstandardNth13gzz + gu22*JacPDstandardNth22gxz +
+ gu32*JacPDstandardNth32gxz) + gu32*(2.*gyzL*G212*G333 +
+ gzzL*(-2.*G313*G323 + 2.*G312*G333) - 2.*JacPDstandardNth13gyz +
+ JacPDstandardNth21gzz + JacPDstandardNth31gyz + JacPDstandardNth33gxy)
+ - 2.*((G122*(gxxL*G113 + gxyL*G213) + G222*(gxyL*G113 +
+ gyzL*G313))*gu22 + (G123*(gxxL*G113 + gxyL*G213) +
+ G223*(gxyL*G113 + gyzL*G313))*gu32 + (gxzL*G113 +
+ gyzL*G213)*(G312*gu21 + G322*gu22 + G323*gu32) + G313*((gyzL*G212 +
+ gzzL*G312)*gu21 + gxzL*(G122*gu22 + G123*gu32)) +
+ gu31*(JacPDstandardNth13gxz + gxxL*SQR(G113))) +
+ gyyL*((-2.*G212*G213 + 2.*G211*G223)*gu21 - 2.*G213*(G222*gu22 +
+ G223*gu32) + gu31*(2.*G211*G233 - 2.*SQR(G213))) +
+ gu31*(-4.*G113*(gxyL*G213 + gxzL*G313) + gyzL*(-4.*G213*G313 +
+ 2.*(G233*G311 + G211*G333)) + JacPDstandardNth11gzz +
+ JacPDstandardNth33gxx + gzzL*(2.*G311*G333 - 2.*SQR(G313))));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R22 = 0.5*(-4*(gyzL*G213*G322*gu31 +
- gzzL*G313*G322*gu31) + gyzL*(4*G212*G312*gu11 - 2*G211*G322*gu11 +
- 4*G223*G312*gu31 - 2*G233*G322*gu33) + 4*(G112*(gxyL*G212*gu11 +
- gxzL*G312*gu11 + gxxL*G123*gu31 + gxyL*G223*gu31 +
- gxzL*G323*gu31) + gxyL*G123*G223*gu33 + gxzL*G123*G323*gu33 +
- gyzL*G223*G323*gu33) + gu21*(JacPDstandardNth12gyy -
- JacPDstandardNth21gyy) - 2*(gxyL*G122*G211*gu11 +
- gyyL*G211*G222*gu11 + gxzL*G122*G311*gu11 + gyzL*G222*G311*gu11 +
- gzzL*G311*G322*gu11 + G111*(gxxL*G122 + gxyL*G222 +
- gxzL*G322)*gu11 + gxxL*G122*G133*gu33 + gxyL*G133*G222*gu33 +
- gxyL*G122*G233*gu33 + gyyL*G222*G233*gu33 + gxzL*G133*G322*gu33 +
- gxzL*G122*G333*gu33 + gyzL*G222*G333*gu33 + gzzL*G322*G333*gu33 +
- gu31*JacPDstandardNth22gxz) + gu31*(-4*gxxL*G113*G122 +
- 4*gxyL*G123*G212 - 4*gxyL*G122*G213 - 4*gxyL*G113*G222 -
- 4*gyyL*G213*G222 + 4*gyyL*G212*G223 + 4*gxzL*G123*G312 -
- 4*gxzL*G122*G313 - 4*gyzL*G222*G313 - 4*gxzL*G113*G322 +
- 4*gyzL*G212*G323 + 4*gzzL*G312*G323 + 2*JacPDstandardNth12gyz -
- JacPDstandardNth13gyy - JacPDstandardNth31gyy) +
- 2*gu31*JacPDstandardNth32gxy + gu32*(-JacPDstandardNth23gyy +
- JacPDstandardNth32gyy) + gu33*(-JacPDstandardNth22gzz +
- 2*JacPDstandardNth32gyz - JacPDstandardNth33gyy) +
- gu11*(-JacPDstandardNth11gyy + 2*JacPDstandardNth12gxy -
- JacPDstandardNth22gxx + 2*gxxL*SQR(G112)) + 2*gxxL*gu33*SQR(G123) +
- 2*gyyL*gu11*SQR(G212) + 2*gyyL*gu33*SQR(G223) +
- 2*gzzL*gu11*SQR(G312) + 2*gzzL*gu33*SQR(G323));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R22 = 0.5*(gyzL*((4.*G212*G312 -
+ 2.*G211*G322)*gu11 + 4.*G223*G312*gu31 - 2.*G233*G322*gu33) +
+ 4.*(G112*((gxyL*G212 + gxzL*G312)*gu11 + (gxxL*G123 + gxyL*G223
+ + gxzL*G323)*gu31) + (gyzL*G223*G323 + G123*(gxyL*G223 +
+ gxzL*G323))*gu33) + gu21*(JacPDstandardNth12gyy -
+ 1.*JacPDstandardNth21gyy) - 2.*((G211*(gxyL*G122 + gyyL*G222) +
+ G111*(gxxL*G122 + gxyL*G222 + gxzL*G322) + G311*(gxzL*G122 +
+ gyzL*G222 + gzzL*G322))*gu11 + (G133*(gxxL*G122 + gxyL*G222) +
+ (gxyL*G122 + gyyL*G222)*G233 + (gyzL*G222 + gzzL*G322)*G333 +
+ gxzL*(G133*G322 + G122*G333))*gu33 + gu31*JacPDstandardNth22gxz) +
+ gu31*(gxyL*(4.*G123*G212 - 4.*G113*G222) - 4.*(G122*(gxxL*G113 +
+ gxyL*G213 + gxzL*G313) + G222*(gyyL*G213 + gyzL*G313) +
+ (gxzL*G113 + gyzL*G213 + gzzL*G313)*G322) + 4.*(G212*(gyyL*G223
+ + gyzL*G323) + G312*(gxzL*G123 + gzzL*G323)) -
+ 1.*(JacPDstandardNth13gyy + JacPDstandardNth31gyy) +
+ 2.*(JacPDstandardNth12gyz + JacPDstandardNth32gxy)) +
+ gu32*JacPDstandardNth32gyy - 1.*(gu32*JacPDstandardNth23gyy +
+ gu33*(JacPDstandardNth22gzz + JacPDstandardNth33gyy)) +
+ gu11*(-1.*(JacPDstandardNth11gyy + JacPDstandardNth22gxx) +
+ 2.*(JacPDstandardNth12gxy + gxxL*SQR(G112) + gyyL*SQR(G212) +
+ gzzL*SQR(G312))) + 2.*gu33*(JacPDstandardNth32gyz + gxxL*SQR(G123)
+ + gyyL*SQR(G223) + gzzL*SQR(G323)));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R23 = 0.5*(gxyL*((-2*G123*G211 +
- 2*G113*G212)*gu11 + G212*(-2*G123*gu21 + 2*G133*gu31) + (-4*G123*G223 +
- 2*(G133*G222 + G122*G233))*gu32) + 2*(G313*((gyzL*G212 +
- gzzL*G312)*gu11 + gxzL*G122*gu21) + G312*((gxzL*G113 +
- gyzL*G213)*gu11 + gxzL*G133*gu31) + G112*((gxxL*G113 +
- gxyL*G213 + gxzL*G313)*gu11 + (-(gxxL*G123) - gxyL*G223 -
- gxzL*G323)*gu21 + gxxL*G133*gu31 + gxyL*G233*gu31 +
- gxzL*G333*gu31) + (gxzL*G122 + gzzL*G322)*G333*gu32 +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R23 = 0.5*(gxyL*((-2.*G123*G211 +
+ 2.*G113*G212)*gu11 + G212*(-2.*G123*gu21 + 2.*G133*gu31) +
+ (-4.*G123*G223 + 2.*(G133*G222 + G122*G233))*gu32) +
+ 2.*(G313*((gyzL*G212 + gzzL*G312)*gu11 + gxzL*G122*gu21) +
+ G312*((gxzL*G113 + gyzL*G213)*gu11 + gxzL*G133*gu31) +
+ G112*((gxxL*G113 + gxyL*G213 + gxzL*G313)*gu11 - 1.*(gxxL*G123
+ + gxyL*G223 + gxzL*G323)*gu21 + (gxxL*G133 + gxyL*G233 +
+ gxzL*G333)*gu31) + (gxzL*G122 + gzzL*G322)*G333*gu32 +
G122*((gxxL*G113 + gxyL*G213)*gu21 + gxxL*G133*gu32) +
G322*((gxzL*G113 + gyzL*G213)*gu21 + gxzL*G133*gu32) +
G233*((gyyL*G212 + gyzL*G312)*gu31 + gyzL*G322*gu32) +
G222*((gxyL*G113 + gyzL*G313)*gu21 + gyzL*G333*gu32)) +
- gu11*(-JacPDstandardNth11gyz + JacPDstandardNth12gxz +
- JacPDstandardNth13gxy - JacPDstandardNth23gxx) +
- gu21*(-2*gyzL*G212*G323 + gzzL*(2*G313*G322 - 2*G312*G323) +
- JacPDstandardNth13gyy - JacPDstandardNth21gyz + JacPDstandardNth22gxz -
- JacPDstandardNth23gxy) + gu33*(-JacPDstandardNth23gzz +
- JacPDstandardNth32gzz) + gu31*(2*gyzL*G212*G333 +
- gzzL*(-2*G313*G323 + 2*G312*G333) + JacPDstandardNth12gzz -
- JacPDstandardNth31gyz + JacPDstandardNth32gxz + JacPDstandardNth33gxy)
- - 2*(((gxzL*G123 + gyzL*G223)*G311 + G111*(gxxL*G123 +
- gxyL*G223 + gxzL*G323))*gu11 + (gxzL*G123*G312 +
- G223*(gyyL*G212 + gyzL*G312))*gu21 + G323*((gyzL*G211 +
- gzzL*G311)*gu11 + gxzL*G113*gu31) + gu31*(G123*(gxxL*G113 +
- gxyL*G213 + gxzL*G313) + G223*(gxyL*G113 + gyzL*G313) +
- gyzL*G213*G323 + JacPDstandardNth23gxz) + gxxL*gu32*SQR(G123)) +
- gyyL*((2*G212*G213 - 2*G211*G223)*gu11 + G213*(2*G222*gu21 -
- 2*G223*gu31) + gu32*(2*G222*G233 - 2*SQR(G223))) +
- gu32*(-4*(gxzL*G123 + gyzL*G223)*G323 + JacPDstandardNth22gzz +
- JacPDstandardNth33gyy - 2*(JacPDstandardNth23gyz + gzzL*SQR(G323))));
+ gu11*(JacPDstandardNth12gxz + JacPDstandardNth13gxy -
+ 1.*(JacPDstandardNth11gyz + JacPDstandardNth23gxx)) +
+ gu21*(-2.*gyzL*G212*G323 + gzzL*(2.*G313*G322 - 2.*G312*G323) +
+ JacPDstandardNth13gyy + JacPDstandardNth22gxz -
+ 1.*(JacPDstandardNth21gyz + JacPDstandardNth23gxy)) +
+ gu33*(-1.*JacPDstandardNth23gzz + JacPDstandardNth32gzz) +
+ gu31*(2.*gyzL*G212*G333 + gzzL*(-2.*G313*G323 + 2.*G312*G333) +
+ JacPDstandardNth12gzz - 1.*JacPDstandardNth31gyz +
+ JacPDstandardNth32gxz + JacPDstandardNth33gxy) - 2.*(((gxzL*G123 +
+ gyzL*G223)*G311 + G111*(gxxL*G123 + gxyL*G223 +
+ gxzL*G323))*gu11 + (gxzL*G123*G312 + G223*(gyyL*G212 +
+ gyzL*G312))*gu21 + G323*((gyzL*G211 + gzzL*G311)*gu11 +
+ gxzL*G113*gu31) + gu31*(G123*(gxxL*G113 + gxyL*G213 +
+ gxzL*G313) + G223*(gxyL*G113 + gyzL*G313) + gyzL*G213*G323 +
+ JacPDstandardNth23gxz) + gxxL*gu32*SQR(G123)) + gyyL*((2.*G212*G213
+ - 2.*G211*G223)*gu11 + G213*(2.*G222*gu21 - 2.*G223*gu31) +
+ gu32*(2.*G222*G233 - 2.*SQR(G223))) + gu32*(-4.*(gxzL*G123 +
+ gyzL*G223)*G323 + JacPDstandardNth22gzz + JacPDstandardNth33gyy -
+ 2.*(JacPDstandardNth23gyz + gzzL*SQR(G323))));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED R33 = 0.5*(4*((G213*(gxyL*G123 +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED R33 = 0.5*(4.*((G213*(gxyL*G123 +
gyyL*G223) + gyzL*(G223*G313 + G213*G323))*gu21 +
G113*((gxyL*G213 + gxzL*G313)*gu11 + (gxxL*G123 + gxyL*G223 +
- gxzL*G323)*gu21) + gxyL*G123*G223*gu22) + gxzL*((-4*G133*G312 +
- 4*G123*G313)*gu21 - 2*G133*G322*gu22) + gyzL*((4*G213*G313 -
- 2*G211*G333)*gu11 + G233*(-4*G312*gu21 - 2*G322*gu22)) +
- gu21*(4*gzzL*G313*G323 - 4*(G133*(gxxL*G112 + gxyL*G212) +
- (gxyL*G112 + gyyL*G212)*G233 + gxzL*G112*G333 + (gyzL*G212 +
- gzzL*G312)*G333) - JacPDstandardNth12gzz + 2*JacPDstandardNth13gyz -
- JacPDstandardNth21gzz + 2*JacPDstandardNth23gxz) +
- gu31*(JacPDstandardNth13gzz - JacPDstandardNth31gzz) +
- gu32*(JacPDstandardNth23gzz - JacPDstandardNth32gzz) -
- 2*((G211*(gxyL*G133 + gyyL*G233) + G111*(gxxL*G133 + gxyL*G233
+ gxzL*G323)*gu21) + gxyL*G123*G223*gu22) + gxzL*((-4.*G133*G312 +
+ 4.*G123*G313)*gu21 - 2.*G133*G322*gu22) + gyzL*((4.*G213*G313 -
+ 2.*G211*G333)*gu11 + G233*(-4.*G312*gu21 - 2.*G322*gu22)) +
+ gu21*(4.*gzzL*G313*G323 - 4.*(G133*(gxxL*G112 + gxyL*G212) +
+ (gxyL*G112 + gyyL*G212)*G233 + (gxzL*G112 + gyzL*G212 +
+ gzzL*G312)*G333) - 1.*(JacPDstandardNth12gzz + JacPDstandardNth21gzz)
+ + 2.*(JacPDstandardNth13gyz + JacPDstandardNth23gxz)) +
+ gu31*(JacPDstandardNth13gzz - 1.*JacPDstandardNth31gzz) +
+ gu32*(JacPDstandardNth23gzz - 1.*JacPDstandardNth32gzz) -
+ 2.*((G211*(gxyL*G133 + gyyL*G233) + G111*(gxxL*G133 + gxyL*G233
+ gxzL*G333) + G311*(gxzL*G133 + gyzL*G233 + gzzL*G333))*gu11 +
(G133*(gxxL*G122 + gxyL*G222) + gxyL*G122*G233 + (gyzL*G222 +
gzzL*G322)*G333)*gu22 + gu21*JacPDstandardNth33gxy) +
- gu11*(-JacPDstandardNth11gzz + 2*JacPDstandardNth13gxz -
- JacPDstandardNth33gxx + 2*gxxL*SQR(G113) + 2*gyyL*SQR(G213) +
- 2*gzzL*SQR(G313)) + gu22*(4*(gxzL*G123*G323 + gyzL*G223*G323) -
- 2*(gyyL*G222*G233 + gxzL*G122*G333) - JacPDstandardNth22gzz +
- 2*JacPDstandardNth23gyz - JacPDstandardNth33gyy + 2*gxxL*SQR(G123) +
- 2*gyyL*SQR(G223) + 2*gzzL*SQR(G323)));
+ gu11*(-1.*(JacPDstandardNth11gzz + JacPDstandardNth33gxx) +
+ 2.*(JacPDstandardNth13gxz + gxxL*SQR(G113) + gyyL*SQR(G213) +
+ gzzL*SQR(G313))) + gu22*(4.*(gxzL*G123 + gyzL*G223)*G323 -
+ 2.*(gyyL*G222*G233 + gxzL*G122*G333) - 1.*(JacPDstandardNth22gzz +
+ JacPDstandardNth33gyy) + 2.*(JacPDstandardNth23gyz + gxxL*SQR(G123) +
+ gyyL*SQR(G223) + gzzL*SQR(G323))));
CCTK_REAL CCTK_ATTRIBUTE_UNUSED trR = gu11*R11 + gu22*R22 +
- 2*(gu21*R12 + gu31*R13 + gu32*R23) + gu33*R33;
+ 2.*(gu21*R12 + gu31*R13 + gu32*R23) + gu33*R33;
CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km11 = kxxL*gu11 + kxyL*gu21 +
kxzL*gu31;
@@ -1284,61 +1278,62 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const
CCTK_REAL CCTK_ATTRIBUTE_UNUSED trK = Km11 + Km22 + Km33;
CCTK_REAL CCTK_ATTRIBUTE_UNUSED rho = INV(SQR(alpL))*(eTttL -
- 2*(betayL*eTtyL + betazL*eTtzL) + 2*(betaxL*(-eTtxL +
+ 2.*(betayL*eTtyL + betazL*eTtzL) + 2.*(betaxL*(-1.*eTtxL +
betayL*eTxyL + betazL*eTxzL) + betayL*betazL*eTyzL) +
eTxxL*SQR(betaxL) + eTyyL*SQR(betayL) + eTzzL*SQR(betazL));
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1 = (-eTtxL + betaxL*eTxxL +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1 = (-1.*eTtxL + betaxL*eTxxL +
betayL*eTxyL + betazL*eTxzL)*INV(alpL);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-eTtyL + betaxL*eTxyL +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-1.*eTtyL + betaxL*eTxyL +
betayL*eTyyL + betazL*eTyzL)*INV(alpL);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-eTtzL + betaxL*eTxzL +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-1.*eTtzL + betaxL*eTxzL +
betayL*eTyzL + betazL*eTzzL)*INV(alpL);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED HL = -2*(Km12*Km21 + Km13*Km31 +
- Km23*Km32) - 16*Pi*rho + trR - SQR(Km11) - SQR(Km22) - SQR(Km33) +
- SQR(trK);
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED HL = -2.*(Km12*Km21 + Km13*Km31 +
+ Km23*Km32) - 50.26548245743669*rho + trR - 1.*(SQR(Km11) + SQR(Km22) +
+ SQR(Km33)) + SQR(trK);
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED M1L = gu21*(-(kxxL*G112) +
- kyyL*G211 + kxyL*(G111 - G212) + kyzL*G311 - kxzL*G312 -
- JacPDstandardNth1kxy + JacPDstandardNth2kxx) + gu22*(-(kxxL*G122) +
- kyyL*G212 + kxyL*(G112 - G222) + kyzL*G312 - kxzL*G322 -
- JacPDstandardNth1kyy + JacPDstandardNth2kxy) + gu31*(-(kxxL*G113) +
- kyzL*G211 - kxyL*G213 + kzzL*G311 + kxzL*(G111 - G313) -
- JacPDstandardNth1kxz + JacPDstandardNth3kxx) + gu32*(kyyL*G213 +
- kxyL*(G113 - 2*G223) + kzzL*G312 + kyzL*(G212 + G313) +
- kxzL*(G112 - 2*G323) - 2*(kxxL*G123 + JacPDstandardNth1kyz) +
- JacPDstandardNth2kxz + JacPDstandardNth3kxy) + gu33*(-(kxxL*G133) +
- kyzL*G213 - kxyL*G233 + kzzL*G313 + kxzL*(G113 - G333) -
- JacPDstandardNth1kzz + JacPDstandardNth3kxz) - 8*Pi*S1;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED M1L = gu21*(kyyL*G211 + kxyL*(G111
+ - 1.*G212) + kyzL*G311 - 1.*(kxxL*G112 + kxzL*G312 +
+ JacPDstandardNth1kxy) + JacPDstandardNth2kxx) + gu22*(kyyL*G212 +
+ kxyL*(G112 - 1.*G222) + kyzL*G312 - 1.*(kxxL*G122 + kxzL*G322 +
+ JacPDstandardNth1kyy) + JacPDstandardNth2kxy) + gu31*(kyzL*G211 +
+ kzzL*G311 + kxzL*(G111 - 1.*G313) - 1.*(kxxL*G113 + kxyL*G213 +
+ JacPDstandardNth1kxz) + JacPDstandardNth3kxx) + gu32*(kyyL*G213 +
+ kxyL*(G113 - 2.*G223) + kzzL*G312 + kyzL*(G212 + G313) +
+ kxzL*(G112 - 2.*G323) - 2.*(kxxL*G123 + JacPDstandardNth1kyz) +
+ JacPDstandardNth2kxz + JacPDstandardNth3kxy) + gu33*(kyzL*G213 +
+ kzzL*G313 + kxzL*(G113 - 1.*G333) - 1.*(kxxL*G133 + kxyL*G233 +
+ JacPDstandardNth1kzz) + JacPDstandardNth3kxz) - 25.132741228718345*S1;
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED M2L = gu11*(kxxL*G112 - kyyL*G211
- + kxyL*(-G111 + G212) - kyzL*G311 + kxzL*G312 +
- JacPDstandardNth1kxy - JacPDstandardNth2kxx) + gu21*(kxxL*G122 -
- kyyL*G212 + kxyL*(-G112 + G222) - kyzL*G312 + kxzL*G322 +
- JacPDstandardNth1kyy - JacPDstandardNth2kxy) + gu31*(kxxL*G123 +
- kxyL*G223 + kzzL*G312 + kyzL*(G212 - 2*G313) + kxzL*(G112 +
- G323) + JacPDstandardNth1kyz - 2*(kxyL*G113 + kyyL*G213 +
- JacPDstandardNth2kxz) + JacPDstandardNth3kxy) + gu32*(kxzL*G122 -
- kxyL*G123 - kyyL*G223 + kzzL*G322 + kyzL*(G222 - G323) -
- JacPDstandardNth2kyz + JacPDstandardNth3kyy) + gu33*(kxzL*G123 -
- kxyL*G133 - kyyL*G233 + kzzL*G323 + kyzL*(G223 - G333) -
- JacPDstandardNth2kzz + JacPDstandardNth3kyz) - 8*Pi*S2;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED M2L = gu11*(kxxL*G112 + kxyL*G212
+ + kxzL*G312 + JacPDstandardNth1kxy - 1.*(kxyL*G111 + kyyL*G211 +
+ kyzL*G311 + JacPDstandardNth2kxx)) + gu21*(kxxL*G122 + kxyL*G222
+ + kxzL*G322 + JacPDstandardNth1kyy - 1.*(kxyL*G112 + kyyL*G212 +
+ kyzL*G312 + JacPDstandardNth2kxy)) + gu31*(kxxL*G123 + kxyL*G223
+ + kzzL*G312 + kyzL*(G212 - 2.*G313) + kxzL*(G112 + G323) +
+ JacPDstandardNth1kyz - 2.*(kxyL*G113 + kyyL*G213 +
+ JacPDstandardNth2kxz) + JacPDstandardNth3kxy) + gu32*(kxzL*G122 +
+ kzzL*G322 + kyzL*(G222 - 1.*G323) - 1.*(kxyL*G123 + kyyL*G223 +
+ JacPDstandardNth2kyz) + JacPDstandardNth3kyy) + gu33*(kxzL*G123 +
+ kzzL*G323 + kyzL*(G223 - 1.*G333) - 1.*(kxyL*G133 + kyyL*G233 +
+ JacPDstandardNth2kzz) + JacPDstandardNth3kyz) - 25.132741228718345*S2;
- CCTK_REAL CCTK_ATTRIBUTE_UNUSED M3L = gu11*(kxxL*G113 - kyzL*G211
- + kxyL*G213 - kzzL*G311 + kxzL*(-G111 + G313) +
- JacPDstandardNth1kxz - JacPDstandardNth3kxx) + gu21*(kxxL*G123 +
- kyyL*G213 + kxyL*(G113 + G223) + kyzL*G313 + kxzL*G323 +
- JacPDstandardNth1kyz + JacPDstandardNth2kxz - 2*(kxzL*G112 +
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED M3L = gu11*(kxxL*G113 + kxyL*G213
+ + kxzL*G313 + JacPDstandardNth1kxz - 1.*(kxzL*G111 + kyzL*G211 +
+ kzzL*G311 + JacPDstandardNth3kxx)) + gu21*(kxxL*G123 + kyyL*G213
+ + kxyL*(G113 + G223) + kyzL*G313 + kxzL*G323 +
+ JacPDstandardNth1kyz + JacPDstandardNth2kxz - 2.*(kxzL*G112 +
kyzL*G212 + kzzL*G312 + JacPDstandardNth3kxy)) + gu31*(kxxL*G133
- - kyzL*G213 + kxyL*G233 - kzzL*G313 + kxzL*(-G113 + G333) +
- JacPDstandardNth1kzz - JacPDstandardNth3kxz) + gu22*(-(kxzL*G122) +
- kxyL*G123 + kyyL*G223 - kzzL*G322 + kyzL*(-G222 + G323) +
- JacPDstandardNth2kyz - JacPDstandardNth3kyy) + gu32*(-(kxzL*G123) +
- kxyL*G133 + kyyL*G233 - kzzL*G323 + kyzL*(-G223 + G333) +
- JacPDstandardNth2kzz - JacPDstandardNth3kyz) - 8*Pi*S3;
+ + kxyL*G233 + kxzL*G333 + JacPDstandardNth1kzz - 1.*(kxzL*G113 +
+ kyzL*G213 + kzzL*G313 + JacPDstandardNth3kxz)) + gu22*(kxyL*G123
+ + kyyL*G223 + kyzL*G323 + JacPDstandardNth2kyz - 1.*(kxzL*G122 +
+ kyzL*G222 + kzzL*G322 + JacPDstandardNth3kyy)) + gu32*(kxyL*G133
+ + kyyL*G233 + kyzL*G333 + JacPDstandardNth2kzz - 1.*(kxzL*G123 +
+ kyzL*G223 + kzzL*G323 + JacPDstandardNth3kyz)) -
+ 25.132741228718345*S3;
/* Copy local copies back to grid functions */
H[index] = HL;
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-16 18:19:35 +0000