Z4co: initialize masked gamt chi and alp to one in Z4co_set_rhs.hxx

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);