Regenerate code (Kranc commit 10daf4a83a783a812c7b5a3aba0df115f0eac98f)

1 files changed, 239 insertions, 250 deletions

diff --git a/ML_ADMConstraints/src/ML_ADMConstraints.cc b/ML_ADMConstraints/src/ML_ADMConstraints.cc
index 829a4cf..4c72e7b 100644
--- a/ML_ADMConstraints/src/ML_ADMConstraints.cc
+++ b/ML_ADMConstraints/src/ML_ADMConstraints.cc
@@ -63,13 +63,15 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
   const CCTK_REAL dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
   const CCTK_REAL dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
   const CCTK_REAL khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
-  const CCTK_REAL kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
-  const CCTK_REAL ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
-  const CCTK_REAL kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
-  const CCTK_REAL keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
-  const CCTK_REAL hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
-  const CCTK_REAL hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
-  const CCTK_REAL hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
+  const CCTK_REAL kthird CCTK_ATTRIBUTE_UNUSED = 
+    0.333333333333333333333333333333;
+  const CCTK_REAL ktwothird CCTK_ATTRIBUTE_UNUSED = 
+    0.666666666666666666666666666667;
+  const CCTK_REAL kfourthird CCTK_ATTRIBUTE_UNUSED = 
+    1.33333333333333333333333333333;
+  const CCTK_REAL hdxi CCTK_ATTRIBUTE_UNUSED = 0.5*dxi;
+  const CCTK_REAL hdyi CCTK_ATTRIBUTE_UNUSED = 0.5*dyi;
+  const CCTK_REAL hdzi CCTK_ATTRIBUTE_UNUSED = 0.5*dzi;
   
   /* Initialize predefined quantities */
   const CCTK_REAL p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx);
@@ -411,104 +413,104 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
     
     if (use_jacobian)
     {
-      JacPDstandardNth1gxx = J11L*PDstandardNth1gxx + 
-        J21L*PDstandardNth2gxx + J31L*PDstandardNth3gxx;
+      JacPDstandardNth1gxx = J11L*PDstandardNth1gxx + J21L*PDstandardNth2gxx 
+        + J31L*PDstandardNth3gxx;
       
-      JacPDstandardNth1gxy = J11L*PDstandardNth1gxy + 
-        J21L*PDstandardNth2gxy + J31L*PDstandardNth3gxy;
+      JacPDstandardNth1gxy = J11L*PDstandardNth1gxy + J21L*PDstandardNth2gxy 
+        + J31L*PDstandardNth3gxy;
       
-      JacPDstandardNth1gxz = J11L*PDstandardNth1gxz + 
-        J21L*PDstandardNth2gxz + J31L*PDstandardNth3gxz;
+      JacPDstandardNth1gxz = J11L*PDstandardNth1gxz + J21L*PDstandardNth2gxz 
+        + J31L*PDstandardNth3gxz;
       
-      JacPDstandardNth1gyy = J11L*PDstandardNth1gyy + 
-        J21L*PDstandardNth2gyy + J31L*PDstandardNth3gyy;
+      JacPDstandardNth1gyy = J11L*PDstandardNth1gyy + J21L*PDstandardNth2gyy 
+        + J31L*PDstandardNth3gyy;
       
-      JacPDstandardNth1gyz = J11L*PDstandardNth1gyz + 
-        J21L*PDstandardNth2gyz + J31L*PDstandardNth3gyz;
+      JacPDstandardNth1gyz = J11L*PDstandardNth1gyz + J21L*PDstandardNth2gyz 
+        + J31L*PDstandardNth3gyz;
       
-      JacPDstandardNth1gzz = J11L*PDstandardNth1gzz + 
-        J21L*PDstandardNth2gzz + J31L*PDstandardNth3gzz;
+      JacPDstandardNth1gzz = J11L*PDstandardNth1gzz + J21L*PDstandardNth2gzz 
+        + J31L*PDstandardNth3gzz;
       
-      JacPDstandardNth1kxy = J11L*PDstandardNth1kxy + 
-        J21L*PDstandardNth2kxy + J31L*PDstandardNth3kxy;
+      JacPDstandardNth1kxy = J11L*PDstandardNth1kxy + J21L*PDstandardNth2kxy 
+        + J31L*PDstandardNth3kxy;
       
-      JacPDstandardNth1kxz = J11L*PDstandardNth1kxz + 
-        J21L*PDstandardNth2kxz + J31L*PDstandardNth3kxz;
+      JacPDstandardNth1kxz = J11L*PDstandardNth1kxz + J21L*PDstandardNth2kxz 
+        + J31L*PDstandardNth3kxz;
       
-      JacPDstandardNth1kyy = J11L*PDstandardNth1kyy + 
-        J21L*PDstandardNth2kyy + J31L*PDstandardNth3kyy;
+      JacPDstandardNth1kyy = J11L*PDstandardNth1kyy + J21L*PDstandardNth2kyy 
+        + J31L*PDstandardNth3kyy;
       
-      JacPDstandardNth1kyz = J11L*PDstandardNth1kyz + 
-        J21L*PDstandardNth2kyz + J31L*PDstandardNth3kyz;
+      JacPDstandardNth1kyz = J11L*PDstandardNth1kyz + J21L*PDstandardNth2kyz 
+        + J31L*PDstandardNth3kyz;
       
-      JacPDstandardNth1kzz = J11L*PDstandardNth1kzz + 
-        J21L*PDstandardNth2kzz + J31L*PDstandardNth3kzz;
+      JacPDstandardNth1kzz = J11L*PDstandardNth1kzz + J21L*PDstandardNth2kzz 
+        + J31L*PDstandardNth3kzz;
       
-      JacPDstandardNth2gxx = J12L*PDstandardNth1gxx + 
-        J22L*PDstandardNth2gxx + J32L*PDstandardNth3gxx;
+      JacPDstandardNth2gxx = J12L*PDstandardNth1gxx + J22L*PDstandardNth2gxx 
+        + J32L*PDstandardNth3gxx;
       
-      JacPDstandardNth2gxy = J12L*PDstandardNth1gxy + 
-        J22L*PDstandardNth2gxy + J32L*PDstandardNth3gxy;
+      JacPDstandardNth2gxy = J12L*PDstandardNth1gxy + J22L*PDstandardNth2gxy 
+        + J32L*PDstandardNth3gxy;
       
-      JacPDstandardNth2gxz = J12L*PDstandardNth1gxz + 
-        J22L*PDstandardNth2gxz + J32L*PDstandardNth3gxz;
+      JacPDstandardNth2gxz = J12L*PDstandardNth1gxz + J22L*PDstandardNth2gxz 
+        + J32L*PDstandardNth3gxz;
       
-      JacPDstandardNth2gyy = J12L*PDstandardNth1gyy + 
-        J22L*PDstandardNth2gyy + J32L*PDstandardNth3gyy;
+      JacPDstandardNth2gyy = J12L*PDstandardNth1gyy + J22L*PDstandardNth2gyy 
+        + J32L*PDstandardNth3gyy;
       
-      JacPDstandardNth2gyz = J12L*PDstandardNth1gyz + 
-        J22L*PDstandardNth2gyz + J32L*PDstandardNth3gyz;
+      JacPDstandardNth2gyz = J12L*PDstandardNth1gyz + J22L*PDstandardNth2gyz 
+        + J32L*PDstandardNth3gyz;
       
-      JacPDstandardNth2gzz = J12L*PDstandardNth1gzz + 
-        J22L*PDstandardNth2gzz + J32L*PDstandardNth3gzz;
+      JacPDstandardNth2gzz = J12L*PDstandardNth1gzz + J22L*PDstandardNth2gzz 
+        + J32L*PDstandardNth3gzz;
       
-      JacPDstandardNth2kxx = J12L*PDstandardNth1kxx + 
-        J22L*PDstandardNth2kxx + J32L*PDstandardNth3kxx;
+      JacPDstandardNth2kxx = J12L*PDstandardNth1kxx + J22L*PDstandardNth2kxx 
+        + J32L*PDstandardNth3kxx;
       
-      JacPDstandardNth2kxy = J12L*PDstandardNth1kxy + 
-        J22L*PDstandardNth2kxy + J32L*PDstandardNth3kxy;
+      JacPDstandardNth2kxy = J12L*PDstandardNth1kxy + J22L*PDstandardNth2kxy 
+        + J32L*PDstandardNth3kxy;
       
-      JacPDstandardNth2kxz = J12L*PDstandardNth1kxz + 
-        J22L*PDstandardNth2kxz + J32L*PDstandardNth3kxz;
+      JacPDstandardNth2kxz = J12L*PDstandardNth1kxz + J22L*PDstandardNth2kxz 
+        + J32L*PDstandardNth3kxz;
       
-      JacPDstandardNth2kyz = J12L*PDstandardNth1kyz + 
-        J22L*PDstandardNth2kyz + J32L*PDstandardNth3kyz;
+      JacPDstandardNth2kyz = J12L*PDstandardNth1kyz + J22L*PDstandardNth2kyz 
+        + J32L*PDstandardNth3kyz;
       
-      JacPDstandardNth2kzz = J12L*PDstandardNth1kzz + 
-        J22L*PDstandardNth2kzz + J32L*PDstandardNth3kzz;
+      JacPDstandardNth2kzz = J12L*PDstandardNth1kzz + J22L*PDstandardNth2kzz 
+        + J32L*PDstandardNth3kzz;
       
-      JacPDstandardNth3gxx = J13L*PDstandardNth1gxx + 
-        J23L*PDstandardNth2gxx + J33L*PDstandardNth3gxx;
+      JacPDstandardNth3gxx = J13L*PDstandardNth1gxx + J23L*PDstandardNth2gxx 
+        + J33L*PDstandardNth3gxx;
       
-      JacPDstandardNth3gxy = J13L*PDstandardNth1gxy + 
-        J23L*PDstandardNth2gxy + J33L*PDstandardNth3gxy;
+      JacPDstandardNth3gxy = J13L*PDstandardNth1gxy + J23L*PDstandardNth2gxy 
+        + J33L*PDstandardNth3gxy;
       
-      JacPDstandardNth3gxz = J13L*PDstandardNth1gxz + 
-        J23L*PDstandardNth2gxz + J33L*PDstandardNth3gxz;
+      JacPDstandardNth3gxz = J13L*PDstandardNth1gxz + J23L*PDstandardNth2gxz 
+        + J33L*PDstandardNth3gxz;
       
-      JacPDstandardNth3gyy = J13L*PDstandardNth1gyy + 
-        J23L*PDstandardNth2gyy + J33L*PDstandardNth3gyy;
+      JacPDstandardNth3gyy = J13L*PDstandardNth1gyy + J23L*PDstandardNth2gyy 
+        + J33L*PDstandardNth3gyy;
       
-      JacPDstandardNth3gyz = J13L*PDstandardNth1gyz + 
-        J23L*PDstandardNth2gyz + J33L*PDstandardNth3gyz;
+      JacPDstandardNth3gyz = J13L*PDstandardNth1gyz + J23L*PDstandardNth2gyz 
+        + J33L*PDstandardNth3gyz;
       
-      JacPDstandardNth3gzz = J13L*PDstandardNth1gzz + 
-        J23L*PDstandardNth2gzz + J33L*PDstandardNth3gzz;
+      JacPDstandardNth3gzz = J13L*PDstandardNth1gzz + J23L*PDstandardNth2gzz 
+        + J33L*PDstandardNth3gzz;
       
-      JacPDstandardNth3kxx = J13L*PDstandardNth1kxx + 
-        J23L*PDstandardNth2kxx + J33L*PDstandardNth3kxx;
+      JacPDstandardNth3kxx = J13L*PDstandardNth1kxx + J23L*PDstandardNth2kxx 
+        + J33L*PDstandardNth3kxx;
       
-      JacPDstandardNth3kxy = J13L*PDstandardNth1kxy + 
-        J23L*PDstandardNth2kxy + J33L*PDstandardNth3kxy;
+      JacPDstandardNth3kxy = J13L*PDstandardNth1kxy + J23L*PDstandardNth2kxy 
+        + J33L*PDstandardNth3kxy;
       
-      JacPDstandardNth3kxz = J13L*PDstandardNth1kxz + 
-        J23L*PDstandardNth2kxz + J33L*PDstandardNth3kxz;
+      JacPDstandardNth3kxz = J13L*PDstandardNth1kxz + J23L*PDstandardNth2kxz 
+        + J33L*PDstandardNth3kxz;
       
-      JacPDstandardNth3kyy = J13L*PDstandardNth1kyy + 
-        J23L*PDstandardNth2kyy + J33L*PDstandardNth3kyy;
+      JacPDstandardNth3kyy = J13L*PDstandardNth1kyy + J23L*PDstandardNth2kyy 
+        + J33L*PDstandardNth3kyy;
       
-      JacPDstandardNth3kyz = J13L*PDstandardNth1kyz + 
-        J23L*PDstandardNth2kyz + J33L*PDstandardNth3kyz;
+      JacPDstandardNth3kyz = J13L*PDstandardNth1kyz + J23L*PDstandardNth2kyz 
+        + J33L*PDstandardNth3kyz;
       
       JacPDstandardNth11gyy = dJ111L*PDstandardNth1gyy + 
         2*(J11L*(J21L*PDstandardNth12gyy + J31L*PDstandardNth13gyy) + 
@@ -1011,107 +1013,106 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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 detg CCTK_ATTRIBUTE_UNUSED = 2*gxyL*gxzL*gyzL + 
+      gzzL*(gxxL*gyyL - SQR(gxyL)) - gyyL*SQR(gxzL) - gxxL*SQR(gyzL);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu11 = INV(detg)*(gyyL*gzzL - 
+    CCTK_REAL gu11 CCTK_ATTRIBUTE_UNUSED = INV(detg)*(gyyL*gzzL - 
       SQR(gyzL));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu21 = (gxzL*gyzL - 
+    CCTK_REAL gu21 CCTK_ATTRIBUTE_UNUSED = (gxzL*gyzL - 
       gxyL*gzzL)*INV(detg);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu31 = (-(gxzL*gyyL) + 
+    CCTK_REAL gu31 CCTK_ATTRIBUTE_UNUSED = (-(gxzL*gyyL) + 
       gxyL*gyzL)*INV(detg);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu22 = INV(detg)*(gxxL*gzzL - 
+    CCTK_REAL gu22 CCTK_ATTRIBUTE_UNUSED = INV(detg)*(gxxL*gzzL - 
       SQR(gxzL));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu32 = (gxyL*gxzL - 
+    CCTK_REAL gu32 CCTK_ATTRIBUTE_UNUSED = (gxyL*gxzL - 
       gxxL*gyzL)*INV(detg);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED gu33 = INV(detg)*(gxxL*gyyL - 
+    CCTK_REAL gu33 CCTK_ATTRIBUTE_UNUSED = INV(detg)*(gxxL*gyyL - 
       SQR(gxyL));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G111 = 0.5*(gu11*JacPDstandardNth1gxx 
+    CCTK_REAL G111 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu11*JacPDstandardNth1gxx 
       + 2*(gu21*JacPDstandardNth1gxy + gu31*JacPDstandardNth1gxz) - 
       gu21*JacPDstandardNth2gxx - gu31*JacPDstandardNth3gxx);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G211 = 0.5*(gu21*JacPDstandardNth1gxx 
+    CCTK_REAL G211 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu21*JacPDstandardNth1gxx 
       + 2*(gu22*JacPDstandardNth1gxy + gu32*JacPDstandardNth1gxz) - 
       gu22*JacPDstandardNth2gxx - gu32*JacPDstandardNth3gxx);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G311 = 0.5*(gu31*JacPDstandardNth1gxx 
+    CCTK_REAL G311 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu31*JacPDstandardNth1gxx 
       + 2*(gu32*JacPDstandardNth1gxy + gu33*JacPDstandardNth1gxz) - 
       gu32*JacPDstandardNth2gxx - gu33*JacPDstandardNth3gxx);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G112 = 0.5*(gu21*JacPDstandardNth1gyy 
+    CCTK_REAL G112 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu21*JacPDstandardNth1gyy 
       + gu11*JacPDstandardNth2gxx + gu31*(JacPDstandardNth1gyz + 
       JacPDstandardNth2gxz - JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G212 = 0.5*(gu22*JacPDstandardNth1gyy 
+    CCTK_REAL G212 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu22*JacPDstandardNth1gyy 
       + gu21*JacPDstandardNth2gxx + gu32*(JacPDstandardNth1gyz + 
       JacPDstandardNth2gxz - JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G312 = 0.5*(gu32*JacPDstandardNth1gyy 
+    CCTK_REAL G312 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu32*JacPDstandardNth1gyy 
       + gu31*JacPDstandardNth2gxx + gu33*(JacPDstandardNth1gyz + 
       JacPDstandardNth2gxz - JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G113 = 0.5*(gu31*JacPDstandardNth1gzz 
+    CCTK_REAL G113 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu31*JacPDstandardNth1gzz 
       + gu11*JacPDstandardNth3gxx + gu21*(JacPDstandardNth1gyz - 
       JacPDstandardNth2gxz + JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G213 = 0.5*(gu32*JacPDstandardNth1gzz 
+    CCTK_REAL G213 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu32*JacPDstandardNth1gzz 
       + gu21*JacPDstandardNth3gxx + gu22*(JacPDstandardNth1gyz - 
       JacPDstandardNth2gxz + JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G313 = 0.5*(gu33*JacPDstandardNth1gzz 
+    CCTK_REAL G313 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu33*JacPDstandardNth1gzz 
       + gu31*JacPDstandardNth3gxx + gu32*(JacPDstandardNth1gyz - 
       JacPDstandardNth2gxz + JacPDstandardNth3gxy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G122 = 
+    CCTK_REAL G122 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu11*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) + 
       gu21*JacPDstandardNth2gyy + gu31*(2*JacPDstandardNth2gyz - 
       JacPDstandardNth3gyy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G222 = 
+    CCTK_REAL G222 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu21*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) + 
       gu22*JacPDstandardNth2gyy + gu32*(2*JacPDstandardNth2gyz - 
       JacPDstandardNth3gyy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G322 = 
+    CCTK_REAL G322 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu31*(-JacPDstandardNth1gyy + 2*JacPDstandardNth2gxy) + 
       gu32*JacPDstandardNth2gyy + gu33*(2*JacPDstandardNth2gyz - 
       JacPDstandardNth3gyy));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G123 = 0.5*(gu31*JacPDstandardNth2gzz 
+    CCTK_REAL G123 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu31*JacPDstandardNth2gzz 
       + gu11*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz + 
       JacPDstandardNth3gxy) + gu21*JacPDstandardNth3gyy);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G223 = 0.5*(gu32*JacPDstandardNth2gzz 
+    CCTK_REAL G223 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu32*JacPDstandardNth2gzz 
       + gu21*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz + 
       JacPDstandardNth3gxy) + gu22*JacPDstandardNth3gyy);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G323 = 0.5*(gu33*JacPDstandardNth2gzz 
+    CCTK_REAL G323 CCTK_ATTRIBUTE_UNUSED = 0.5*(gu33*JacPDstandardNth2gzz 
       + gu31*(-JacPDstandardNth1gyz + JacPDstandardNth2gxz + 
       JacPDstandardNth3gxy) + gu32*JacPDstandardNth3gyy);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G133 = 
+    CCTK_REAL G133 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu11*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) + 
       gu21*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) + 
       gu31*JacPDstandardNth3gzz);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G233 = 
+    CCTK_REAL G233 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu21*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) + 
       gu22*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) + 
       gu32*JacPDstandardNth3gzz);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED G333 = 
+    CCTK_REAL G333 CCTK_ATTRIBUTE_UNUSED = 
       0.5*(gu31*(-JacPDstandardNth1gzz + 2*JacPDstandardNth3gxz) + 
       gu32*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) + 
       gu33*JacPDstandardNth3gzz);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED R11 = 0.5*(gxyL*(-4*G123*G211*gu32 + 
+    CCTK_REAL R11 CCTK_ATTRIBUTE_UNUSED = 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 + 
@@ -1120,33 +1121,30 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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) + 
+      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));
+      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 R12 = 0.5*(gxzL*((-2*G123*G312 + 
+    CCTK_REAL R12 CCTK_ATTRIBUTE_UNUSED = 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 + 
-      gzzL*G322)*gu21 + (gxzL*G123 + gzzL*G323)*gu31) + 
-      (G122*(gxxL*G113 + gxyL*G213) + G222*(gxyL*G113 + gyzL*G313) + 
-      (gxzL*G113 + gyzL*G213 + gzzL*G313)*G322)*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 + gzzL*G322)*gu21 + (gxzL*G123 + 
+      gzzL*G323)*gu31) + (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 + 
       JacPDstandardNth21gyy) + gu32*(JacPDstandardNth21gyz + 
@@ -1157,28 +1155,26 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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)) + 
+      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 - 
+    CCTK_REAL R13 CCTK_ATTRIBUTE_UNUSED = 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)) + 
+      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 - 
@@ -1187,8 +1183,8 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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 + 
+      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 + 
@@ -1196,22 +1192,20 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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));
+      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 R22 = 0.5*(-4*(gyzL*G213*G322*gu31 + 
+    CCTK_REAL R22 CCTK_ATTRIBUTE_UNUSED = 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 + 
+      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 + 
@@ -1225,56 +1219,52 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       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));
+      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 R23 = 0.5*(gxyL*((-2*G123*G211 + 
+    CCTK_REAL R23 CCTK_ATTRIBUTE_UNUSED = 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 + 
-      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) + 
+      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 + 
+      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))));
+      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))));
     
-    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*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) + 
+    CCTK_REAL R33 CCTK_ATTRIBUTE_UNUSED = 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*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) - 
@@ -1282,94 +1272,93 @@ static void ML_ADMConstraints_Body(const cGH* restrict const cctkGH, const int d
       2*JacPDstandardNth23gyz - JacPDstandardNth33gyy + 2*gxxL*SQR(G123) + 
       2*gyyL*SQR(G223) + 2*gzzL*SQR(G323)));
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED trR = gu11*R11 + gu22*R22 + 
+    CCTK_REAL trR CCTK_ATTRIBUTE_UNUSED = gu11*R11 + gu22*R22 + 
       2*(gu21*R12 + gu31*R13 + gu32*R23) + gu33*R33;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km11 = kxxL*gu11 + kxyL*gu21 + 
+    CCTK_REAL Km11 CCTK_ATTRIBUTE_UNUSED = kxxL*gu11 + kxyL*gu21 + 
       kxzL*gu31;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km21 = kxxL*gu21 + kxyL*gu22 + 
+    CCTK_REAL Km21 CCTK_ATTRIBUTE_UNUSED = kxxL*gu21 + kxyL*gu22 + 
       kxzL*gu32;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km31 = kxxL*gu31 + kxyL*gu32 + 
+    CCTK_REAL Km31 CCTK_ATTRIBUTE_UNUSED = kxxL*gu31 + kxyL*gu32 + 
       kxzL*gu33;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km12 = kxyL*gu11 + kyyL*gu21 + 
+    CCTK_REAL Km12 CCTK_ATTRIBUTE_UNUSED = kxyL*gu11 + kyyL*gu21 + 
       kyzL*gu31;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km22 = kxyL*gu21 + kyyL*gu22 + 
+    CCTK_REAL Km22 CCTK_ATTRIBUTE_UNUSED = kxyL*gu21 + kyyL*gu22 + 
       kyzL*gu32;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km32 = kxyL*gu31 + kyyL*gu32 + 
+    CCTK_REAL Km32 CCTK_ATTRIBUTE_UNUSED = kxyL*gu31 + kyyL*gu32 + 
       kyzL*gu33;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km13 = kxzL*gu11 + kyzL*gu21 + 
+    CCTK_REAL Km13 CCTK_ATTRIBUTE_UNUSED = kxzL*gu11 + kyzL*gu21 + 
       kzzL*gu31;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km23 = kxzL*gu21 + kyzL*gu22 + 
+    CCTK_REAL Km23 CCTK_ATTRIBUTE_UNUSED = kxzL*gu21 + kyzL*gu22 + 
       kzzL*gu32;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Km33 = kxzL*gu31 + kyzL*gu32 + 
+    CCTK_REAL Km33 CCTK_ATTRIBUTE_UNUSED = kxzL*gu31 + kyzL*gu32 + 
       kzzL*gu33;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED trK = Km11 + Km22 + Km33;
+    CCTK_REAL trK CCTK_ATTRIBUTE_UNUSED = Km11 + Km22 + Km33;
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rho = INV(SQR(alpL))*(eTttL - 
-      2*(betayL*eTtyL + betazL*eTtzL) + 2*(betaxL*(-eTtxL + 
-      betayL*eTxyL + betazL*eTxzL) + betayL*betazL*eTyzL) + 
-      eTxxL*SQR(betaxL) + eTyyL*SQR(betayL) + eTzzL*SQR(betazL));
+    CCTK_REAL rho CCTK_ATTRIBUTE_UNUSED = INV(SQR(alpL))*(eTttL - 
+      2*(betayL*eTtyL + betazL*eTtzL) + 2*(betaxL*(-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 S1 CCTK_ATTRIBUTE_UNUSED = (-eTtxL + betaxL*eTxxL + 
       betayL*eTxyL + betazL*eTxzL)*INV(alpL);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-eTtyL + betaxL*eTxyL + 
+    CCTK_REAL S2 CCTK_ATTRIBUTE_UNUSED = (-eTtyL + betaxL*eTxyL + 
       betayL*eTyyL + betazL*eTyzL)*INV(alpL);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-eTtzL + betaxL*eTxzL + 
+    CCTK_REAL S3 CCTK_ATTRIBUTE_UNUSED = (-eTtzL + betaxL*eTxzL + 
       betayL*eTyzL + betazL*eTzzL)*INV(alpL);
     
-    CCTK_REAL CCTK_ATTRIBUTE_UNUSED HL = -2*(Km12*Km21 + Km13*Km31 + 
+    CCTK_REAL HL CCTK_ATTRIBUTE_UNUSED = -2*(Km12*Km21 + Km13*Km31 + 
       Km23*Km32) - 16*Pi*rho + trR - 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 M1L CCTK_ATTRIBUTE_UNUSED = 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 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 M2L CCTK_ATTRIBUTE_UNUSED = 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 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 + 
-      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;
+    CCTK_REAL M3L CCTK_ATTRIBUTE_UNUSED = 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 + 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;
     
     /* Copy local copies back to grid functions */
     H[index] = HL;