Z4co: add Z4c_set_rhs.wl

Z4co/wolfram/Z4c_set_rhs.hxx

@@ -0,0 +1,93 @@
+/* Z4c_set_rhs.hxx */
+/* Produced with Mathematica */
+
+vreal &chi = gf_chi(mask, index5).elts[0];
+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 &exAt11 = gf_exAt(mask, index5).elts[0];
+vreal &exAt12 = gf_exAt(mask, index5).elts[1];
+vreal &exAt13 = gf_exAt(mask, index5).elts[2];
+vreal &exAt22 = gf_exAt(mask, index5).elts[3];
+vreal &exAt23 = gf_exAt(mask, index5).elts[4];
+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 &beta1 = gf_beta(mask, index5).elts[0];
+vreal &beta2 = gf_beta(mask, index5).elts[1];
+vreal &beta3 = gf_beta(mask, index5).elts[2];
+vreal &dchi1 = gf_dchi(mask, index5).elts[0];
+vreal &dchi2 = gf_dchi(mask, index5).elts[1];
+vreal &dchi3 = gf_dchi(mask, index5).elts[2];
+vreal &dgamt111 = gf_dgamt(mask, index5).elts[0];
+vreal &dgamt112 = gf_dgamt(mask, index5).elts[1];
+vreal &dgamt113 = gf_dgamt(mask, index5).elts[2];
+vreal &dgamt122 = gf_dgamt(mask, index5).elts[3];
+vreal &dgamt123 = gf_dgamt(mask, index5).elts[4];
+vreal &dgamt133 = gf_dgamt(mask, index5).elts[5];
+vreal &dgamt211 = gf_dgamt(mask, index5).elts[6];
+vreal &dgamt212 = gf_dgamt(mask, index5).elts[7];
+vreal &dgamt213 = gf_dgamt(mask, index5).elts[8];
+vreal &dgamt222 = gf_dgamt(mask, index5).elts[9];
+vreal &dgamt223 = gf_dgamt(mask, index5).elts[10];
+vreal &dgamt233 = gf_dgamt(mask, index5).elts[11];
+vreal &dgamt311 = gf_dgamt(mask, index5).elts[12];
+vreal &dgamt312 = gf_dgamt(mask, index5).elts[13];
+vreal &dgamt313 = gf_dgamt(mask, index5).elts[14];
+vreal &dgamt322 = gf_dgamt(mask, index5).elts[15];
+vreal &dgamt323 = gf_dgamt(mask, index5).elts[16];
+vreal &dgamt333 = gf_dgamt(mask, index5).elts[17];
+vreal &dexKh1 = gf_dexKh(mask, index5).elts[0];
+vreal &dexKh2 = gf_dexKh(mask, index5).elts[1];
+vreal &dexKh3 = gf_dexKh(mask, index5).elts[2];
+vreal &dexAt111 = gf_dexAt(mask, index5).elts[0];
+vreal &dexAt112 = gf_dexAt(mask, index5).elts[1];
+vreal &dexAt113 = gf_dexAt(mask, index5).elts[2];
+vreal &dexAt122 = gf_dexAt(mask, index5).elts[3];
+vreal &dexAt123 = gf_dexAt(mask, index5).elts[4];
+vreal &dexAt133 = gf_dexAt(mask, index5).elts[5];
+vreal &dexAt211 = gf_dexAt(mask, index5).elts[6];
+vreal &dexAt212 = gf_dexAt(mask, index5).elts[7];
+vreal &dexAt213 = gf_dexAt(mask, index5).elts[8];
+vreal &dexAt222 = gf_dexAt(mask, index5).elts[9];
+vreal &dexAt223 = gf_dexAt(mask, index5).elts[10];
+vreal &dexAt233 = gf_dexAt(mask, index5).elts[11];
+vreal &dexAt311 = gf_dexAt(mask, index5).elts[12];
+vreal &dexAt312 = gf_dexAt(mask, index5).elts[13];
+vreal &dexAt313 = gf_dexAt(mask, index5).elts[14];
+vreal &dexAt322 = gf_dexAt(mask, index5).elts[15];
+vreal &dexAt323 = gf_dexAt(mask, index5).elts[16];
+vreal &dexAt333 = gf_dexAt(mask, index5).elts[17];
+vreal &dtrGt11 = gf_dtrGt(mask, index5).elts[0];
+vreal &dtrGt12 = gf_dtrGt(mask, index5).elts[1];
+vreal &dtrGt13 = gf_dtrGt(mask, index5).elts[2];
+vreal &dtrGt21 = gf_dtrGt(mask, index5).elts[3];
+vreal &dtrGt22 = gf_dtrGt(mask, index5).elts[4];
+vreal &dtrGt23 = gf_dtrGt(mask, index5).elts[5];
+vreal &dtrGt31 = gf_dtrGt(mask, index5).elts[6];
+vreal &dtrGt32 = gf_dtrGt(mask, index5).elts[7];
+vreal &dtrGt33 = gf_dtrGt(mask, index5).elts[8];
+vreal &dTheta1 = gf_dTheta(mask, index5).elts[0];
+vreal &dTheta2 = gf_dTheta(mask, index5).elts[1];
+vreal &dTheta3 = gf_dTheta(mask, index5).elts[2];
+vreal &dalpha1 = gf_dalpha(mask, index5).elts[0];
+vreal &dalpha2 = gf_dalpha(mask, index5).elts[1];
+vreal &dalpha3 = gf_dalpha(mask, index5).elts[2];
+vreal &dbeta11 = gf_dbeta(mask, index5).elts[0];
+vreal &dbeta12 = gf_dbeta(mask, index5).elts[1];
+vreal &dbeta13 = gf_dbeta(mask, index5).elts[2];
+vreal &dbeta21 = gf_dbeta(mask, index5).elts[3];
+vreal &dbeta22 = gf_dbeta(mask, index5).elts[4];
+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];
+
+/* Z4c_set_rhs.hxx */

Z4co/wolfram/Z4c_set_rhs.wl

@@ -0,0 +1,34 @@
+(* ::Package:: *)
+
+(* Z4c_set_rhs.wl *)
+
+(* (c) Liwei Ji, 07/2024 *)
+
+Needs["xAct`xCoba`", FileNameJoin[{Environment["GENERATO"], "src/Generato.wl"
+  }]]
+
+SetPVerbose[False];
+
+SetPrintDate[False];
+
+SetGridPointIndex[""];
+
+SetTempVariableType["vreal"];
+
+DefManifold[M3, 3, IndexRange[a, z]];
+
+DefChart[cart, M3, {1, 2, 3}, {X[], Y[], Z[]}, ChartColor -> Blue];
+
+(* Define Variables *)
+
+<<Z4c_vars.wl
+
+SetOutputFile[FileNameJoin[{Directory[], "Z4c_set_rhs.hxx"}]];
+
+$MainPrint[] :=
+  Module[{},
+    PrintInitializations[{Mode -> "MainIn"}, EvolVarlist];
+    PrintInitializations[{Mode -> "MainIn"}, dEvolVarlist];
+  ];
+
+Import[FileNameJoin[{Environment["GENERATO"], "codes/CarpetX.wl"}]];