From 66d5ec743a8448aa46f51c7c161d21922fe90f32 Mon Sep 17 00:00:00 2001
From: Liwei Ji <jiliwei.phys@gmail.com>
Date: Fri, 5 Jul 2024 11:27:23 -0400
Subject: [PATCH] Z4co: initialize masked gamt chi and alp to one in
 Z4co_set_rhs.hxx

---
 Z4co/wolfram/Z4co_set_rhs.hxx | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/Z4co/wolfram/Z4co_set_rhs.hxx b/Z4co/wolfram/Z4co_set_rhs.hxx
index a224d9c7..f1a27967 100644
--- a/Z4co/wolfram/Z4co_set_rhs.hxx
+++ b/Z4co/wolfram/Z4co_set_rhs.hxx
@@ -43,13 +43,13 @@ const vreal &eT23 = gf_eT(mask, index2)(1,2);
 const vreal &eT31 = gf_eT(mask, index2)(2,0);
 const vreal &eT32 = gf_eT(mask, index2)(2,1);
 const vreal &eT33 = gf_eT(mask, index2)(2,2);
-const vreal chi = tl_chi(mask, index5);
-const vreal gamt11 = tl_gamt(mask, index5)(0,0);
-const vreal gamt12 = tl_gamt(mask, index5)(0,1);
-const vreal gamt13 = tl_gamt(mask, index5)(0,2);
-const vreal gamt22 = tl_gamt(mask, index5)(1,1);
-const vreal gamt23 = tl_gamt(mask, index5)(1,2);
-const vreal gamt33 = tl_gamt(mask, index5)(2,2);
+const vreal chi = tl_chi(mask, index5, 1);
+const vreal gamt11 = tl_gamt(mask, index5, 1)(0,0);
+const vreal gamt12 = tl_gamt(mask, index5, 0)(0,1);
+const vreal gamt13 = tl_gamt(mask, index5, 0)(0,2);
+const vreal gamt22 = tl_gamt(mask, index5, 1)(1,1);
+const vreal gamt23 = tl_gamt(mask, index5, 0)(1,2);
+const vreal gamt33 = tl_gamt(mask, index5, 1)(2,2);
 const vreal exKh = tl_exKh(mask, index5);
 const vreal exAt11 = tl_exAt(mask, index5)(0,0);
 const vreal exAt12 = tl_exAt(mask, index5)(0,1);
@@ -61,7 +61,7 @@ const vreal trGt1 = tl_trGt(mask, index5)(0);
 const vreal trGt2 = tl_trGt(mask, index5)(1);
 const vreal trGt3 = tl_trGt(mask, index5)(2);
 const vreal Theta = tl_Theta(mask, index5);
-const vreal alpha = tl_alpha(mask, index5);
+const vreal alpha = tl_alpha(mask, index5, 1);
 const vreal beta1 = tl_beta(mask, index5)(0);
 const vreal beta2 = tl_beta(mask, index5)(1);
 const vreal beta3 = tl_beta(mask, index5)(2);