From cb828b1b9b2e41e545154b34a8526e55acd458cf Mon Sep 17 00:00:00 2001
From: Liwei Ji <jiliwei.phys@gmail.com>
Date: Tue, 19 Nov 2024 10:43:36 -0500
Subject: [PATCH] ADMconstraints: treat DexK separatedly to remove warnings

---
 ADMconstraints/wolfram/ADM_set_constraint.hxx | 13 -------------
 ADMconstraints/wolfram/ADM_set_constraint.wl  | 12 +++++-------
 2 files changed, 5 insertions(+), 20 deletions(-)

diff --git a/ADMconstraints/wolfram/ADM_set_constraint.hxx b/ADMconstraints/wolfram/ADM_set_constraint.hxx
index f244272b..b49f76d1 100644
--- a/ADMconstraints/wolfram/ADM_set_constraint.hxx
+++ b/ADMconstraints/wolfram/ADM_set_constraint.hxx
@@ -1084,10 +1084,6 @@ exK11*invgam11 + 2*exK12*invgam12 + 2*exK13*invgam13 + exK22*invgam22 +
   2*exK23*invgam23 + exK33*invgam33
 ;
 
-vreal DexK111
-=
-dexK111 - 2*(exK11*Gam111 + exK12*Gam211 + exK13*Gam311)
-;
 
 vreal DexK112
 =
@@ -1134,11 +1130,6 @@ dexK213 - exK11*Gam123 - exK23*Gam212 - exK12*Gam223 - exK33*Gam312 -
   exK13*(Gam112 + Gam323)
 ;
 
-vreal DexK222
-=
-dexK222 - 2*(exK12*Gam122 + exK22*Gam222 + exK23*Gam322)
-;
-
 vreal DexK223
 =
 dexK223 - exK13*Gam122 - exK12*Gam123 - exK23*Gam222 - exK22*Gam223 -
@@ -1178,10 +1169,6 @@ dexK323 - exK13*Gam123 - exK12*Gam133 - exK23*Gam223 - exK22*Gam233 -
   exK33*Gam323 - exK23*Gam333
 ;
 
-vreal DexK333
-=
-dexK333 - 2*(exK13*Gam133 + exK23*Gam233 + exK33*Gam333)
-;
 
 vreal rho
 =
diff --git a/ADMconstraints/wolfram/ADM_set_constraint.wl b/ADMconstraints/wolfram/ADM_set_constraint.wl
index f51ef720..61881bd2 100644
--- a/ADMconstraints/wolfram/ADM_set_constraint.wl
+++ b/ADMconstraints/wolfram/ADM_set_constraint.wl
@@ -4,8 +4,7 @@
 
 (* (c) Liwei Ji, 07/2024 *)
 
-Needs["xAct`xCoba`", FileNameJoin[{Environment["GENERATO"], "src/Generato.wl"
-  }]]
+Needs["xAct`xCoba`", FileNameJoin[{Environment["GENERATO"], "src/Generato.wl"}]]
 
 SetPVerbose[False];
 
@@ -62,15 +61,14 @@ $MainPrint[] :=
                          ddADMVarlist];
     pr[];
     PrintEquations[{Mode -> "Temp"}, Drop[IntermediateVarlist, {-1}]];
-
+    pr[];
+    (* treat DexK separately *)
     Module[{printDK},
-      printDK[kk_,ii_,jj_] := PrintComponentEquations[GetDefaultChart[],
-                                                      DexK[{kk,-GetDefaultChart[]},{ii,-GetDefaultChart[]},{jj,-GetDefaultChart[]}]];
+      printDK[kk_, ii_, jj_] := PrintComponentEquation[GetDefaultChart[], DexK[{kk, -GetDefaultChart[]}, {ii, -GetDefaultChart[]}, {jj, -GetDefaultChart[]}]];
       SetParsePrintCompEQNMode[{NewVar -> True}];
-      Do[If[ii!=jj || ii!=kk || jj!=kk, printDK[kk,ii,jj]], {kk,1,3}, {ii,1,3}, {jj,ii,3}]
+      Do[If[ii != jj || ii != kk || jj != kk, printDK[kk, ii, jj]], {kk, 1, 3}, {ii, 1, 3}, {jj, ii, 3}]
     ];
     pr[];
-
     PrintEquations[{Mode -> "Temp"}, MatterVarlist];
     pr[];
     PrintEquations[{Mode -> "Main"}, ConstraintVarlist];