Minus sign of Weyl moved

GRTLCollaboration · Jan 9, 2025 · 3a9213c · 3a9213c
1 parent cb5a069
commit 3a9213c
Show file tree

Hide file tree

Showing 3 changed files with 15 additions and 7 deletions.
diff --git a/Source/ModifiedCCZ4/ModifiedCCZ4RHS.impl.hpp b/Source/ModifiedCCZ4/ModifiedCCZ4RHS.impl.hpp
@@ -242,8 +242,8 @@ ModifiedCCZ4RHS<theory_t, gauge_t, deriv_t>::get_full_kappa_times_Sij_TF(
     // solve linear system for the theory fields that require it (e.g. 4dST)
     my_theory.solve_lhs(theory_rhs, theory_vars, d1, d2, advec, coords);
 
-    Tensor<2, data_t> out = theory_rhs.A;
-    FOR(i, j) out[i][j] += -rhs.A[i][j];
+    Tensor<2, data_t> out = -theory_rhs.A;
+    FOR(i, j) out[i][j] += rhs.A[i][j];
     FOR(i, j) out[i][j] /= chi_regularised;
 
     return out;

diff --git a/Source/ModifiedCCZ4/ModifiedGravityWeyl4.impl.hpp b/Source/ModifiedCCZ4/ModifiedGravityWeyl4.impl.hpp
@@ -61,7 +61,7 @@ void ModifiedGravityWeyl4<theory_t, gauge_t, deriv_t>::add_theory_EB(
         this->get_full_kappa_times_Sij_TF(vars, d1, d2, advec, coords);
     // as we made the vacuum expression of Bij explictly symmetric and Eij
     // explictly trace-free, only Eij has matter terms
-    FOR(i, j) { ebfields.E[i][j] += 0.5 * kappa_times_Sij_TF[i][j]; }
+    FOR(i, j) { ebfields.E[i][j] += -0.5 * kappa_times_Sij_TF[i][j]; }
 }
 
 #endif /* MODIFIEDGRAVITYWEYL4_IMPL_HPP_ */
diff --git a/Tests/FourDerivScalarTensorTest/values1.hpp b/Tests/FourDerivScalarTensorTest/values1.hpp
@@ -304,12 +304,20 @@ advec.shift[0] = 0.461360733159421;
 advec.shift[1] = 1.0268828982221523;
 advec.shift[2] = 0.6401386881951718;
 double dlapsedt_known = 1.5315265817973822;
-double dshiftdt_known[3] = {0.644927622295347, 1.1036799516420281, 1.1942651362731231};
-double dhdt_known[3][3] = {{-0.41248363281772926, 0.09489728268776099, -1.613067404747195}, {0.09489728268776099, -1.1460244712912242, -0.3891497048282506}, {-1.613067404747195, -0.3891497048282506, -0.941011218218068}};
+double dshiftdt_known[3] = {0.644927622295347, 1.1036799516420281,
+                            1.1942651362731231};
+double dhdt_known[3][3] = {
+    {-0.41248363281772926, 0.09489728268776099, -1.613067404747195},
+    {0.09489728268776099, -1.1460244712912242, -0.3891497048282506},
+    {-1.613067404747195, -0.3891497048282506, -0.941011218218068}};
 double dchidt_known = -0.20789282532664588;
 double dphidt_known = 1.007294909427872;
 double dThetadt_known = 0.24968744636812856;
-double dGammadt_known[3] = {1.7154099903876132, -3.866470912731934, 0.6107228928800519};
-double dAdt_known[3][3] = {{3.0012653143064405, 5.356715436903114, 10.21380403073045}, {5.356715436903113, 5.066192678570908, 6.031056178402607}, {10.213804030730447, 6.031056178402605, 13.35675402658088}};
+double dGammadt_known[3] = {1.7154099903876132, -3.866470912731934,
+                            0.6107228928800519};
+double dAdt_known[3][3] = {
+    {3.0012653143064405, 5.356715436903114, 10.21380403073045},
+    {5.356715436903113, 5.066192678570908, 6.031056178402607},
+    {10.213804030730447, 6.031056178402605, 13.35675402658088}};
 double dKdt_known = -0.9512787601245327;
 double dPidt_known = 6.7881943530110105;