Z4co: organize wl files

Z4co/wolfram/Z4c_set_rhs.hxx

@@ -1,14 +1,14 @@
 /* Z4c_set_rhs.hxx */
 /* Produced with Mathematica */
 
-vreal &chi = gf_chi(mask, index5).elts[0];
+vreal &chi = gf_chi(mask, index5);
 vreal &gamt11 = gf_gamt(mask, index5).elts[0];
 vreal &gamt12 = gf_gamt(mask, index5).elts[1];
 vreal &gamt13 = gf_gamt(mask, index5).elts[2];
 vreal &gamt22 = gf_gamt(mask, index5).elts[3];
 vreal &gamt23 = gf_gamt(mask, index5).elts[4];
 vreal &gamt33 = gf_gamt(mask, index5).elts[5];
-vreal &exKh = gf_exKh(mask, index5).elts[0];
+vreal &exKh = gf_exKh(mask, index5);
 vreal &exAt11 = gf_exAt(mask, index5).elts[0];
 vreal &exAt12 = gf_exAt(mask, index5).elts[1];
 vreal &exAt13 = gf_exAt(mask, index5).elts[2];
@@ -18,8 +18,8 @@ vreal &exAt33 = gf_exAt(mask, index5).elts[5];
 vreal &trGt1 = gf_trGt(mask, index5).elts[0];
 vreal &trGt2 = gf_trGt(mask, index5).elts[1];
 vreal &trGt3 = gf_trGt(mask, index5).elts[2];
-vreal &Theta = gf_Theta(mask, index5).elts[0];
-vreal &alpha = gf_alpha(mask, index5).elts[0];
+vreal &Theta = gf_Theta(mask, index5);
+vreal &alpha = gf_alpha(mask, index5);
 vreal &beta1 = gf_beta(mask, index5).elts[0];
 vreal &beta2 = gf_beta(mask, index5).elts[1];
 vreal &beta3 = gf_beta(mask, index5).elts[2];
@@ -89,5 +89,71 @@ vreal &dbeta23 = gf_dbeta(mask, index5).elts[5];
 vreal &dbeta31 = gf_dbeta(mask, index5).elts[6];
 vreal &dbeta32 = gf_dbeta(mask, index5).elts[7];
 vreal &dbeta33 = gf_dbeta(mask, index5).elts[8];
+vreal &ddchi11 = gf_ddchi(mask, index5).elts[0];
+vreal &ddchi12 = gf_ddchi(mask, index5).elts[1];
+vreal &ddchi13 = gf_ddchi(mask, index5).elts[2];
+vreal &ddchi22 = gf_ddchi(mask, index5).elts[3];
+vreal &ddchi23 = gf_ddchi(mask, index5).elts[4];
+vreal &ddchi33 = gf_ddchi(mask, index5).elts[5];
+vreal &ddgamt1111 = gf_ddgamt(mask, index5).elts[0];
+vreal &ddgamt1112 = gf_ddgamt(mask, index5).elts[1];
+vreal &ddgamt1113 = gf_ddgamt(mask, index5).elts[2];
+vreal &ddgamt1122 = gf_ddgamt(mask, index5).elts[3];
+vreal &ddgamt1123 = gf_ddgamt(mask, index5).elts[4];
+vreal &ddgamt1133 = gf_ddgamt(mask, index5).elts[5];
+vreal &ddgamt1211 = gf_ddgamt(mask, index5).elts[6];
+vreal &ddgamt1212 = gf_ddgamt(mask, index5).elts[7];
+vreal &ddgamt1213 = gf_ddgamt(mask, index5).elts[8];
+vreal &ddgamt1222 = gf_ddgamt(mask, index5).elts[9];
+vreal &ddgamt1223 = gf_ddgamt(mask, index5).elts[10];
+vreal &ddgamt1233 = gf_ddgamt(mask, index5).elts[11];
+vreal &ddgamt1311 = gf_ddgamt(mask, index5).elts[12];
+vreal &ddgamt1312 = gf_ddgamt(mask, index5).elts[13];
+vreal &ddgamt1313 = gf_ddgamt(mask, index5).elts[14];
+vreal &ddgamt1322 = gf_ddgamt(mask, index5).elts[15];
+vreal &ddgamt1323 = gf_ddgamt(mask, index5).elts[16];
+vreal &ddgamt1333 = gf_ddgamt(mask, index5).elts[17];
+vreal &ddgamt2211 = gf_ddgamt(mask, index5).elts[18];
+vreal &ddgamt2212 = gf_ddgamt(mask, index5).elts[19];
+vreal &ddgamt2213 = gf_ddgamt(mask, index5).elts[20];
+vreal &ddgamt2222 = gf_ddgamt(mask, index5).elts[21];
+vreal &ddgamt2223 = gf_ddgamt(mask, index5).elts[22];
+vreal &ddgamt2233 = gf_ddgamt(mask, index5).elts[23];
+vreal &ddgamt2311 = gf_ddgamt(mask, index5).elts[24];
+vreal &ddgamt2312 = gf_ddgamt(mask, index5).elts[25];
+vreal &ddgamt2313 = gf_ddgamt(mask, index5).elts[26];
+vreal &ddgamt2322 = gf_ddgamt(mask, index5).elts[27];
+vreal &ddgamt2323 = gf_ddgamt(mask, index5).elts[28];
+vreal &ddgamt2333 = gf_ddgamt(mask, index5).elts[29];
+vreal &ddgamt3311 = gf_ddgamt(mask, index5).elts[30];
+vreal &ddgamt3312 = gf_ddgamt(mask, index5).elts[31];
+vreal &ddgamt3313 = gf_ddgamt(mask, index5).elts[32];
+vreal &ddgamt3322 = gf_ddgamt(mask, index5).elts[33];
+vreal &ddgamt3323 = gf_ddgamt(mask, index5).elts[34];
+vreal &ddgamt3333 = gf_ddgamt(mask, index5).elts[35];
+vreal &ddalpha11 = gf_ddalpha(mask, index5).elts[0];
+vreal &ddalpha12 = gf_ddalpha(mask, index5).elts[1];
+vreal &ddalpha13 = gf_ddalpha(mask, index5).elts[2];
+vreal &ddalpha22 = gf_ddalpha(mask, index5).elts[3];
+vreal &ddalpha23 = gf_ddalpha(mask, index5).elts[4];
+vreal &ddalpha33 = gf_ddalpha(mask, index5).elts[5];
+vreal &ddbeta111 = gf_ddbeta(mask, index5).elts[0];
+vreal &ddbeta112 = gf_ddbeta(mask, index5).elts[1];
+vreal &ddbeta113 = gf_ddbeta(mask, index5).elts[2];
+vreal &ddbeta121 = gf_ddbeta(mask, index5).elts[3];
+vreal &ddbeta122 = gf_ddbeta(mask, index5).elts[4];
+vreal &ddbeta123 = gf_ddbeta(mask, index5).elts[5];
+vreal &ddbeta131 = gf_ddbeta(mask, index5).elts[6];
+vreal &ddbeta132 = gf_ddbeta(mask, index5).elts[7];
+vreal &ddbeta133 = gf_ddbeta(mask, index5).elts[8];
+vreal &ddbeta221 = gf_ddbeta(mask, index5).elts[9];
+vreal &ddbeta222 = gf_ddbeta(mask, index5).elts[10];
+vreal &ddbeta223 = gf_ddbeta(mask, index5).elts[11];
+vreal &ddbeta231 = gf_ddbeta(mask, index5).elts[12];
+vreal &ddbeta232 = gf_ddbeta(mask, index5).elts[13];
+vreal &ddbeta233 = gf_ddbeta(mask, index5).elts[14];
+vreal &ddbeta331 = gf_ddbeta(mask, index5).elts[15];
+vreal &ddbeta332 = gf_ddbeta(mask, index5).elts[16];
+vreal &ddbeta333 = gf_ddbeta(mask, index5).elts[17];
 
 /* Z4c_set_rhs.hxx */

Z4co/wolfram/Z4c_set_rhs.wl

@@ -21,14 +21,15 @@ DefChart[cart, M3, {1, 2, 3}, {X[], Y[], Z[]}, ChartColor -> Blue];
 
 (* Define Variables *)
 
-<<Z4c_vars.wl
+<<wl/Z4c_vars.wl
 
 SetOutputFile[FileNameJoin[{Directory[], "Z4c_set_rhs.hxx"}]];
 
 $MainPrint[] :=
   Module[{},
     PrintInitializations[{Mode -> "MainIn"}, EvolVarlist];
     PrintInitializations[{Mode -> "MainIn"}, dEvolVarlist];
+    PrintInitializations[{Mode -> "MainIn"}, ddEvolVarlist];
   ];
 
 Import[FileNameJoin[{Environment["GENERATO"], "codes/CarpetX.wl"}]];