Skip to content

Commit

Permalink
ADMconstraints: treat DexK separatedly to remove warnings
Browse files Browse the repository at this point in the history
  • Loading branch information
lwJi committed Nov 19, 2024
1 parent a880d20 commit cb828b1
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 20 deletions.
13 changes: 0 additions & 13 deletions ADMconstraints/wolfram/ADM_set_constraint.hxx
Original file line number Diff line number Diff line change
Expand Up @@ -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
=
Expand Down Expand Up @@ -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 -
Expand Down Expand Up @@ -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
=
Expand Down
12 changes: 5 additions & 7 deletions ADMconstraints/wolfram/ADM_set_constraint.wl
Original file line number Diff line number Diff line change
Expand Up @@ -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];

Expand Down Expand Up @@ -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];
Expand Down

0 comments on commit cb828b1

Please sign in to comment.