Skip to content

Commit

Permalink
Z4co: update Z4co_set_rhs.hxx
Browse files Browse the repository at this point in the history
  • Loading branch information
@lwJi
lwJi committed Jul 3, 2024
1 parent 65e4d25 commit ee6ffbf
Showing 1 changed file with 167 additions and 163 deletions.
330 changes: 167 additions & 163 deletions Z4co/wolfram/Z4co_set_rhs.hxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -43,160 +43,160 @@ 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 &exKh = tl_exKh(mask, index5);
const vreal &exAt11 = tl_exAt(mask, index5)(0,0);
const vreal &exAt12 = tl_exAt(mask, index5)(0,1);
const vreal &exAt13 = tl_exAt(mask, index5)(0,2);
const vreal &exAt22 = tl_exAt(mask, index5)(1,1);
const vreal &exAt23 = tl_exAt(mask, index5)(1,2);
const vreal &exAt33 = tl_exAt(mask, index5)(2,2);
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 &beta1 = tl_beta(mask, index5)(0);
const vreal &beta2 = tl_beta(mask, index5)(1);
const vreal &beta3 = tl_beta(mask, index5)(2);
const vreal &dchi1 = tl_dchi(mask, index5)(0);
const vreal &dchi2 = tl_dchi(mask, index5)(1);
const vreal &dchi3 = tl_dchi(mask, index5)(2);
const vreal &dgamt111 = tl_dgamt(mask, index5)(0,0)(0);
const vreal &dgamt112 = tl_dgamt(mask, index5)(0,1)(0);
const vreal &dgamt113 = tl_dgamt(mask, index5)(0,2)(0);
const vreal &dgamt122 = tl_dgamt(mask, index5)(1,1)(0);
const vreal &dgamt123 = tl_dgamt(mask, index5)(1,2)(0);
const vreal &dgamt133 = tl_dgamt(mask, index5)(2,2)(0);
const vreal &dgamt211 = tl_dgamt(mask, index5)(0,0)(1);
const vreal &dgamt212 = tl_dgamt(mask, index5)(0,1)(1);
const vreal &dgamt213 = tl_dgamt(mask, index5)(0,2)(1);
const vreal &dgamt222 = tl_dgamt(mask, index5)(1,1)(1);
const vreal &dgamt223 = tl_dgamt(mask, index5)(1,2)(1);
const vreal &dgamt233 = tl_dgamt(mask, index5)(2,2)(1);
const vreal &dgamt311 = tl_dgamt(mask, index5)(0,0)(2);
const vreal &dgamt312 = tl_dgamt(mask, index5)(0,1)(2);
const vreal &dgamt313 = tl_dgamt(mask, index5)(0,2)(2);
const vreal &dgamt322 = tl_dgamt(mask, index5)(1,1)(2);
const vreal &dgamt323 = tl_dgamt(mask, index5)(1,2)(2);
const vreal &dgamt333 = tl_dgamt(mask, index5)(2,2)(2);
const vreal &dexKh1 = tl_dexKh(mask, index5)(0);
const vreal &dexKh2 = tl_dexKh(mask, index5)(1);
const vreal &dexKh3 = tl_dexKh(mask, index5)(2);
const vreal &dexAt111 = tl_dexAt(mask, index5)(0,0)(0);
const vreal &dexAt112 = tl_dexAt(mask, index5)(0,1)(0);
const vreal &dexAt113 = tl_dexAt(mask, index5)(0,2)(0);
const vreal &dexAt122 = tl_dexAt(mask, index5)(1,1)(0);
const vreal &dexAt123 = tl_dexAt(mask, index5)(1,2)(0);
const vreal &dexAt133 = tl_dexAt(mask, index5)(2,2)(0);
const vreal &dexAt211 = tl_dexAt(mask, index5)(0,0)(1);
const vreal &dexAt212 = tl_dexAt(mask, index5)(0,1)(1);
const vreal &dexAt213 = tl_dexAt(mask, index5)(0,2)(1);
const vreal &dexAt222 = tl_dexAt(mask, index5)(1,1)(1);
const vreal &dexAt223 = tl_dexAt(mask, index5)(1,2)(1);
const vreal &dexAt233 = tl_dexAt(mask, index5)(2,2)(1);
const vreal &dexAt311 = tl_dexAt(mask, index5)(0,0)(2);
const vreal &dexAt312 = tl_dexAt(mask, index5)(0,1)(2);
const vreal &dexAt313 = tl_dexAt(mask, index5)(0,2)(2);
const vreal &dexAt322 = tl_dexAt(mask, index5)(1,1)(2);
const vreal &dexAt323 = tl_dexAt(mask, index5)(1,2)(2);
const vreal &dexAt333 = tl_dexAt(mask, index5)(2,2)(2);
const vreal &dtrGt11 = tl_dtrGt(mask, index5)(0)(0);
const vreal &dtrGt12 = tl_dtrGt(mask, index5)(1)(0);
const vreal &dtrGt13 = tl_dtrGt(mask, index5)(2)(0);
const vreal &dtrGt21 = tl_dtrGt(mask, index5)(0)(1);
const vreal &dtrGt22 = tl_dtrGt(mask, index5)(1)(1);
const vreal &dtrGt23 = tl_dtrGt(mask, index5)(2)(1);
const vreal &dtrGt31 = tl_dtrGt(mask, index5)(0)(2);
const vreal &dtrGt32 = tl_dtrGt(mask, index5)(1)(2);
const vreal &dtrGt33 = tl_dtrGt(mask, index5)(2)(2);
const vreal &dTheta1 = tl_dTheta(mask, index5)(0);
const vreal &dTheta2 = tl_dTheta(mask, index5)(1);
const vreal &dTheta3 = tl_dTheta(mask, index5)(2);
const vreal &dalpha1 = tl_dalpha(mask, index5)(0);
const vreal &dalpha2 = tl_dalpha(mask, index5)(1);
const vreal &dalpha3 = tl_dalpha(mask, index5)(2);
const vreal &dbeta11 = tl_dbeta(mask, index5)(0)(0);
const vreal &dbeta12 = tl_dbeta(mask, index5)(1)(0);
const vreal &dbeta13 = tl_dbeta(mask, index5)(2)(0);
const vreal &dbeta21 = tl_dbeta(mask, index5)(0)(1);
const vreal &dbeta22 = tl_dbeta(mask, index5)(1)(1);
const vreal &dbeta23 = tl_dbeta(mask, index5)(2)(1);
const vreal &dbeta31 = tl_dbeta(mask, index5)(0)(2);
const vreal &dbeta32 = tl_dbeta(mask, index5)(1)(2);
const vreal &dbeta33 = tl_dbeta(mask, index5)(2)(2);
const vreal &ddchi11 = tl_ddchi(mask, index5)(0,0);
const vreal &ddchi12 = tl_ddchi(mask, index5)(0,1);
const vreal &ddchi13 = tl_ddchi(mask, index5)(0,2);
const vreal &ddchi22 = tl_ddchi(mask, index5)(1,1);
const vreal &ddchi23 = tl_ddchi(mask, index5)(1,2);
const vreal &ddchi33 = tl_ddchi(mask, index5)(2,2);
const vreal &ddgamt1111 = tl_ddgamt(mask, index5)(0,0)(0,0);
const vreal &ddgamt1112 = tl_ddgamt(mask, index5)(0,1)(0,0);
const vreal &ddgamt1113 = tl_ddgamt(mask, index5)(0,2)(0,0);
const vreal &ddgamt1122 = tl_ddgamt(mask, index5)(1,1)(0,0);
const vreal &ddgamt1123 = tl_ddgamt(mask, index5)(1,2)(0,0);
const vreal &ddgamt1133 = tl_ddgamt(mask, index5)(2,2)(0,0);
const vreal &ddgamt1211 = tl_ddgamt(mask, index5)(0,0)(0,1);
const vreal &ddgamt1212 = tl_ddgamt(mask, index5)(0,1)(0,1);
const vreal &ddgamt1213 = tl_ddgamt(mask, index5)(0,2)(0,1);
const vreal &ddgamt1222 = tl_ddgamt(mask, index5)(1,1)(0,1);
const vreal &ddgamt1223 = tl_ddgamt(mask, index5)(1,2)(0,1);
const vreal &ddgamt1233 = tl_ddgamt(mask, index5)(2,2)(0,1);
const vreal &ddgamt1311 = tl_ddgamt(mask, index5)(0,0)(0,2);
const vreal &ddgamt1312 = tl_ddgamt(mask, index5)(0,1)(0,2);
const vreal &ddgamt1313 = tl_ddgamt(mask, index5)(0,2)(0,2);
const vreal &ddgamt1322 = tl_ddgamt(mask, index5)(1,1)(0,2);
const vreal &ddgamt1323 = tl_ddgamt(mask, index5)(1,2)(0,2);
const vreal &ddgamt1333 = tl_ddgamt(mask, index5)(2,2)(0,2);
const vreal &ddgamt2211 = tl_ddgamt(mask, index5)(0,0)(1,1);
const vreal &ddgamt2212 = tl_ddgamt(mask, index5)(0,1)(1,1);
const vreal &ddgamt2213 = tl_ddgamt(mask, index5)(0,2)(1,1);
const vreal &ddgamt2222 = tl_ddgamt(mask, index5)(1,1)(1,1);
const vreal &ddgamt2223 = tl_ddgamt(mask, index5)(1,2)(1,1);
const vreal &ddgamt2233 = tl_ddgamt(mask, index5)(2,2)(1,1);
const vreal &ddgamt2311 = tl_ddgamt(mask, index5)(0,0)(1,2);
const vreal &ddgamt2312 = tl_ddgamt(mask, index5)(0,1)(1,2);
const vreal &ddgamt2313 = tl_ddgamt(mask, index5)(0,2)(1,2);
const vreal &ddgamt2322 = tl_ddgamt(mask, index5)(1,1)(1,2);
const vreal &ddgamt2323 = tl_ddgamt(mask, index5)(1,2)(1,2);
const vreal &ddgamt2333 = tl_ddgamt(mask, index5)(2,2)(1,2);
const vreal &ddgamt3311 = tl_ddgamt(mask, index5)(0,0)(2,2);
const vreal &ddgamt3312 = tl_ddgamt(mask, index5)(0,1)(2,2);
const vreal &ddgamt3313 = tl_ddgamt(mask, index5)(0,2)(2,2);
const vreal &ddgamt3322 = tl_ddgamt(mask, index5)(1,1)(2,2);
const vreal &ddgamt3323 = tl_ddgamt(mask, index5)(1,2)(2,2);
const vreal &ddgamt3333 = tl_ddgamt(mask, index5)(2,2)(2,2);
const vreal &ddalpha11 = tl_ddalpha(mask, index5)(0,0);
const vreal &ddalpha12 = tl_ddalpha(mask, index5)(0,1);
const vreal &ddalpha13 = tl_ddalpha(mask, index5)(0,2);
const vreal &ddalpha22 = tl_ddalpha(mask, index5)(1,1);
const vreal &ddalpha23 = tl_ddalpha(mask, index5)(1,2);
const vreal &ddalpha33 = tl_ddalpha(mask, index5)(2,2);
const vreal &ddbeta111 = tl_ddbeta(mask, index5)(0)(0,0);
const vreal &ddbeta112 = tl_ddbeta(mask, index5)(1)(0,0);
const vreal &ddbeta113 = tl_ddbeta(mask, index5)(2)(0,0);
const vreal &ddbeta121 = tl_ddbeta(mask, index5)(0)(0,1);
const vreal &ddbeta122 = tl_ddbeta(mask, index5)(1)(0,1);
const vreal &ddbeta123 = tl_ddbeta(mask, index5)(2)(0,1);
const vreal &ddbeta131 = tl_ddbeta(mask, index5)(0)(0,2);
const vreal &ddbeta132 = tl_ddbeta(mask, index5)(1)(0,2);
const vreal &ddbeta133 = tl_ddbeta(mask, index5)(2)(0,2);
const vreal &ddbeta221 = tl_ddbeta(mask, index5)(0)(1,1);
const vreal &ddbeta222 = tl_ddbeta(mask, index5)(1)(1,1);
const vreal &ddbeta223 = tl_ddbeta(mask, index5)(2)(1,1);
const vreal &ddbeta231 = tl_ddbeta(mask, index5)(0)(1,2);
const vreal &ddbeta232 = tl_ddbeta(mask, index5)(1)(1,2);
const vreal &ddbeta233 = tl_ddbeta(mask, index5)(2)(1,2);
const vreal &ddbeta331 = tl_ddbeta(mask, index5)(0)(2,2);
const vreal &ddbeta332 = tl_ddbeta(mask, index5)(1)(2,2);
const vreal &ddbeta333 = tl_ddbeta(mask, index5)(2)(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 exKh = tl_exKh(mask, index5);
const vreal exAt11 = tl_exAt(mask, index5)(0,0);
const vreal exAt12 = tl_exAt(mask, index5)(0,1);
const vreal exAt13 = tl_exAt(mask, index5)(0,2);
const vreal exAt22 = tl_exAt(mask, index5)(1,1);
const vreal exAt23 = tl_exAt(mask, index5)(1,2);
const vreal exAt33 = tl_exAt(mask, index5)(2,2);
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 beta1 = tl_beta(mask, index5)(0);
const vreal beta2 = tl_beta(mask, index5)(1);
const vreal beta3 = tl_beta(mask, index5)(2);
const vreal dchi1 = tl_dchi(mask, index5)(0);
const vreal dchi2 = tl_dchi(mask, index5)(1);
const vreal dchi3 = tl_dchi(mask, index5)(2);
const vreal dgamt111 = tl_dgamt(mask, index5)(0,0)(0);
const vreal dgamt112 = tl_dgamt(mask, index5)(0,1)(0);
const vreal dgamt113 = tl_dgamt(mask, index5)(0,2)(0);
const vreal dgamt122 = tl_dgamt(mask, index5)(1,1)(0);
const vreal dgamt123 = tl_dgamt(mask, index5)(1,2)(0);
const vreal dgamt133 = tl_dgamt(mask, index5)(2,2)(0);
const vreal dgamt211 = tl_dgamt(mask, index5)(0,0)(1);
const vreal dgamt212 = tl_dgamt(mask, index5)(0,1)(1);
const vreal dgamt213 = tl_dgamt(mask, index5)(0,2)(1);
const vreal dgamt222 = tl_dgamt(mask, index5)(1,1)(1);
const vreal dgamt223 = tl_dgamt(mask, index5)(1,2)(1);
const vreal dgamt233 = tl_dgamt(mask, index5)(2,2)(1);
const vreal dgamt311 = tl_dgamt(mask, index5)(0,0)(2);
const vreal dgamt312 = tl_dgamt(mask, index5)(0,1)(2);
const vreal dgamt313 = tl_dgamt(mask, index5)(0,2)(2);
const vreal dgamt322 = tl_dgamt(mask, index5)(1,1)(2);
const vreal dgamt323 = tl_dgamt(mask, index5)(1,2)(2);
const vreal dgamt333 = tl_dgamt(mask, index5)(2,2)(2);
const vreal dexKh1 = tl_dexKh(mask, index5)(0);
const vreal dexKh2 = tl_dexKh(mask, index5)(1);
const vreal dexKh3 = tl_dexKh(mask, index5)(2);
const vreal dexAt111 = tl_dexAt(mask, index5)(0,0)(0);
const vreal dexAt112 = tl_dexAt(mask, index5)(0,1)(0);
const vreal dexAt113 = tl_dexAt(mask, index5)(0,2)(0);
const vreal dexAt122 = tl_dexAt(mask, index5)(1,1)(0);
const vreal dexAt123 = tl_dexAt(mask, index5)(1,2)(0);
const vreal dexAt133 = tl_dexAt(mask, index5)(2,2)(0);
const vreal dexAt211 = tl_dexAt(mask, index5)(0,0)(1);
const vreal dexAt212 = tl_dexAt(mask, index5)(0,1)(1);
const vreal dexAt213 = tl_dexAt(mask, index5)(0,2)(1);
const vreal dexAt222 = tl_dexAt(mask, index5)(1,1)(1);
const vreal dexAt223 = tl_dexAt(mask, index5)(1,2)(1);
const vreal dexAt233 = tl_dexAt(mask, index5)(2,2)(1);
const vreal dexAt311 = tl_dexAt(mask, index5)(0,0)(2);
const vreal dexAt312 = tl_dexAt(mask, index5)(0,1)(2);
const vreal dexAt313 = tl_dexAt(mask, index5)(0,2)(2);
const vreal dexAt322 = tl_dexAt(mask, index5)(1,1)(2);
const vreal dexAt323 = tl_dexAt(mask, index5)(1,2)(2);
const vreal dexAt333 = tl_dexAt(mask, index5)(2,2)(2);
const vreal dtrGt11 = tl_dtrGt(mask, index5)(0)(0);
const vreal dtrGt12 = tl_dtrGt(mask, index5)(1)(0);
const vreal dtrGt13 = tl_dtrGt(mask, index5)(2)(0);
const vreal dtrGt21 = tl_dtrGt(mask, index5)(0)(1);
const vreal dtrGt22 = tl_dtrGt(mask, index5)(1)(1);
const vreal dtrGt23 = tl_dtrGt(mask, index5)(2)(1);
const vreal dtrGt31 = tl_dtrGt(mask, index5)(0)(2);
const vreal dtrGt32 = tl_dtrGt(mask, index5)(1)(2);
const vreal dtrGt33 = tl_dtrGt(mask, index5)(2)(2);
const vreal dTheta1 = tl_dTheta(mask, index5)(0);
const vreal dTheta2 = tl_dTheta(mask, index5)(1);
const vreal dTheta3 = tl_dTheta(mask, index5)(2);
const vreal dalpha1 = tl_dalpha(mask, index5)(0);
const vreal dalpha2 = tl_dalpha(mask, index5)(1);
const vreal dalpha3 = tl_dalpha(mask, index5)(2);
const vreal dbeta11 = tl_dbeta(mask, index5)(0)(0);
const vreal dbeta12 = tl_dbeta(mask, index5)(1)(0);
const vreal dbeta13 = tl_dbeta(mask, index5)(2)(0);
const vreal dbeta21 = tl_dbeta(mask, index5)(0)(1);
const vreal dbeta22 = tl_dbeta(mask, index5)(1)(1);
const vreal dbeta23 = tl_dbeta(mask, index5)(2)(1);
const vreal dbeta31 = tl_dbeta(mask, index5)(0)(2);
const vreal dbeta32 = tl_dbeta(mask, index5)(1)(2);
const vreal dbeta33 = tl_dbeta(mask, index5)(2)(2);
const vreal ddchi11 = tl_ddchi(mask, index5)(0,0);
const vreal ddchi12 = tl_ddchi(mask, index5)(0,1);
const vreal ddchi13 = tl_ddchi(mask, index5)(0,2);
const vreal ddchi22 = tl_ddchi(mask, index5)(1,1);
const vreal ddchi23 = tl_ddchi(mask, index5)(1,2);
const vreal ddchi33 = tl_ddchi(mask, index5)(2,2);
const vreal ddgamt1111 = tl_ddgamt(mask, index5)(0,0)(0,0);
const vreal ddgamt1112 = tl_ddgamt(mask, index5)(0,1)(0,0);
const vreal ddgamt1113 = tl_ddgamt(mask, index5)(0,2)(0,0);
const vreal ddgamt1122 = tl_ddgamt(mask, index5)(1,1)(0,0);
const vreal ddgamt1123 = tl_ddgamt(mask, index5)(1,2)(0,0);
const vreal ddgamt1133 = tl_ddgamt(mask, index5)(2,2)(0,0);
const vreal ddgamt1211 = tl_ddgamt(mask, index5)(0,0)(0,1);
const vreal ddgamt1212 = tl_ddgamt(mask, index5)(0,1)(0,1);
const vreal ddgamt1213 = tl_ddgamt(mask, index5)(0,2)(0,1);
const vreal ddgamt1222 = tl_ddgamt(mask, index5)(1,1)(0,1);
const vreal ddgamt1223 = tl_ddgamt(mask, index5)(1,2)(0,1);
const vreal ddgamt1233 = tl_ddgamt(mask, index5)(2,2)(0,1);
const vreal ddgamt1311 = tl_ddgamt(mask, index5)(0,0)(0,2);
const vreal ddgamt1312 = tl_ddgamt(mask, index5)(0,1)(0,2);
const vreal ddgamt1313 = tl_ddgamt(mask, index5)(0,2)(0,2);
const vreal ddgamt1322 = tl_ddgamt(mask, index5)(1,1)(0,2);
const vreal ddgamt1323 = tl_ddgamt(mask, index5)(1,2)(0,2);
const vreal ddgamt1333 = tl_ddgamt(mask, index5)(2,2)(0,2);
const vreal ddgamt2211 = tl_ddgamt(mask, index5)(0,0)(1,1);
const vreal ddgamt2212 = tl_ddgamt(mask, index5)(0,1)(1,1);
const vreal ddgamt2213 = tl_ddgamt(mask, index5)(0,2)(1,1);
const vreal ddgamt2222 = tl_ddgamt(mask, index5)(1,1)(1,1);
const vreal ddgamt2223 = tl_ddgamt(mask, index5)(1,2)(1,1);
const vreal ddgamt2233 = tl_ddgamt(mask, index5)(2,2)(1,1);
const vreal ddgamt2311 = tl_ddgamt(mask, index5)(0,0)(1,2);
const vreal ddgamt2312 = tl_ddgamt(mask, index5)(0,1)(1,2);
const vreal ddgamt2313 = tl_ddgamt(mask, index5)(0,2)(1,2);
const vreal ddgamt2322 = tl_ddgamt(mask, index5)(1,1)(1,2);
const vreal ddgamt2323 = tl_ddgamt(mask, index5)(1,2)(1,2);
const vreal ddgamt2333 = tl_ddgamt(mask, index5)(2,2)(1,2);
const vreal ddgamt3311 = tl_ddgamt(mask, index5)(0,0)(2,2);
const vreal ddgamt3312 = tl_ddgamt(mask, index5)(0,1)(2,2);
const vreal ddgamt3313 = tl_ddgamt(mask, index5)(0,2)(2,2);
const vreal ddgamt3322 = tl_ddgamt(mask, index5)(1,1)(2,2);
const vreal ddgamt3323 = tl_ddgamt(mask, index5)(1,2)(2,2);
const vreal ddgamt3333 = tl_ddgamt(mask, index5)(2,2)(2,2);
const vreal ddalpha11 = tl_ddalpha(mask, index5)(0,0);
const vreal ddalpha12 = tl_ddalpha(mask, index5)(0,1);
const vreal ddalpha13 = tl_ddalpha(mask, index5)(0,2);
const vreal ddalpha22 = tl_ddalpha(mask, index5)(1,1);
const vreal ddalpha23 = tl_ddalpha(mask, index5)(1,2);
const vreal ddalpha33 = tl_ddalpha(mask, index5)(2,2);
const vreal ddbeta111 = tl_ddbeta(mask, index5)(0)(0,0);
const vreal ddbeta112 = tl_ddbeta(mask, index5)(1)(0,0);
const vreal ddbeta113 = tl_ddbeta(mask, index5)(2)(0,0);
const vreal ddbeta121 = tl_ddbeta(mask, index5)(0)(0,1);
const vreal ddbeta122 = tl_ddbeta(mask, index5)(1)(0,1);
const vreal ddbeta123 = tl_ddbeta(mask, index5)(2)(0,1);
const vreal ddbeta131 = tl_ddbeta(mask, index5)(0)(0,2);
const vreal ddbeta132 = tl_ddbeta(mask, index5)(1)(0,2);
const vreal ddbeta133 = tl_ddbeta(mask, index5)(2)(0,2);
const vreal ddbeta221 = tl_ddbeta(mask, index5)(0)(1,1);
const vreal ddbeta222 = tl_ddbeta(mask, index5)(1)(1,1);
const vreal ddbeta223 = tl_ddbeta(mask, index5)(2)(1,1);
const vreal ddbeta231 = tl_ddbeta(mask, index5)(0)(1,2);
const vreal ddbeta232 = tl_ddbeta(mask, index5)(1)(1,2);
const vreal ddbeta233 = tl_ddbeta(mask, index5)(2)(1,2);
const vreal ddbeta331 = tl_ddbeta(mask, index5)(0)(2,2);
const vreal ddbeta332 = tl_ddbeta(mask, index5)(1)(2,2);
const vreal ddbeta333 = tl_ddbeta(mask, index5)(2)(2,2);

vreal detinvgamt
=
Expand Down Expand Up @@ -1349,13 +1349,14 @@ dttrGt1.store(mask, index2,
7*ddbeta121*invgamt12 + ddbeta222*invgamt12 + ddbeta233*invgamt12 +
7*ddbeta131*invgamt13 + ddbeta232*invgamt13 + ddbeta333*invgamt13 +
3*ddbeta221*invgamt22 + 6*ddbeta231*invgamt23 + 3*ddbeta331*invgamt33 -
2*dbeta11*trGtd1 - 2*dbeta22*trGtd1 - 2*dbeta33*trGtd1 +
dbeta11*trGtd1 + 2*dbeta22*trGtd1 + 2*dbeta33*trGtd1 +
2*alpha*(3*exAtUU11*Gt111 + 6*exAtUU12*Gt112 + 6*exAtUU13*Gt113 +
3*exAtUU22*Gt122 + 6*exAtUU23*Gt123 + 3*exAtUU33*Gt133 -
2*dexKh1*invgamt11 - dTheta1*invgamt11 - 2*dexKh2*invgamt12 -
dTheta2*invgamt12 - 2*dexKh3*invgamt13 - dTheta3*invgamt13 -
24*cpi*invgam11*Sm1 - 24*cpi*invgam12*Sm2 - 24*cpi*invgam13*Sm3 -
3*ckappa1*trGt1 + 3*ckappa1*trGtd1))/3.
24*cpi*invgamt11*Sm1 - 24*cpi*invgamt12*Sm2 -
24*cpi*invgamt13*Sm3 - 3*ckappa1*trGt1 + 3*ckappa1*trGtd1) -
3*dbeta21*trGtd2 - 3*dbeta31*trGtd3)/3.
);

dttrGt2.store(mask, index2,
Expand All @@ -1365,13 +1366,15 @@ dttrGt2.store(mask, index2,
ddbeta133*invgamt12 + 6*ddbeta132*invgamt13 + ddbeta121*invgamt22 +
4*ddbeta222*invgamt22 + ddbeta233*invgamt22 + ddbeta131*invgamt23 +
7*ddbeta232*invgamt23 + ddbeta333*invgamt23 + 3*ddbeta332*invgamt33 -
2*dbeta11*trGtd2 - 2*dbeta22*trGtd2 - 2*dbeta33*trGtd2 +
2*alpha*(3*exAtUU11*Gt211 + 6*exAtUU12*Gt212 + 6*exAtUU13*Gt213 +
3*dbeta12*trGtd1 + 2*dbeta11*trGtd2 - dbeta22*trGtd2 +
2*dbeta33*trGtd2 + 2*alpha*
(3*exAtUU11*Gt211 + 6*exAtUU12*Gt212 + 6*exAtUU13*Gt213 +
3*exAtUU22*Gt222 + 6*exAtUU23*Gt223 + 3*exAtUU33*Gt233 -
2*dexKh1*invgamt12 - dTheta1*invgamt12 - 2*dexKh2*invgamt22 -
dTheta2*invgamt22 - 2*dexKh3*invgamt23 - dTheta3*invgamt23 -
24*cpi*invgam12*Sm1 - 24*cpi*invgam22*Sm2 - 24*cpi*invgam23*Sm3 -
3*ckappa1*trGt2 + 3*ckappa1*trGtd2))/3.
24*cpi*invgamt12*Sm1 - 24*cpi*invgamt22*Sm2 -
24*cpi*invgamt23*Sm3 - 3*ckappa1*trGt2 + 3*ckappa1*trGtd2) -
3*dbeta32*trGtd3)/3.
);

dttrGt3.store(mask, index2,
Expand All @@ -1381,12 +1384,13 @@ dttrGt3.store(mask, index2,
ddbeta122*invgamt13 + 7*ddbeta133*invgamt13 + 3*ddbeta223*invgamt22 +
ddbeta121*invgamt23 + ddbeta222*invgamt23 + 7*ddbeta233*invgamt23 +
ddbeta131*invgamt33 + ddbeta232*invgamt33 + 4*ddbeta333*invgamt33 -
2*dbeta11*trGtd3 - 2*dbeta22*trGtd3 - 2*dbeta33*trGtd3 +
3*dbeta13*trGtd1 - 3*dbeta23*trGtd2 + 2*dbeta11*trGtd3 +
2*dbeta22*trGtd3 - dbeta33*trGtd3 +
2*alpha*(3*exAtUU11*Gt311 + 6*exAtUU12*Gt312 + 6*exAtUU13*Gt313 +
3*exAtUU22*Gt322 + 6*exAtUU23*Gt323 + 3*exAtUU33*Gt333 -
2*dexKh1*invgamt13 - dTheta1*invgamt13 - 2*dexKh2*invgamt23 -
dTheta2*invgamt23 - 2*dexKh3*invgamt33 - dTheta3*invgamt33 -
24*cpi*invgam13*Sm1 - 24*cpi*invgam23*Sm2 - 24*cpi*invgam33*Sm3 -
24*cpi*invgamt13*Sm1 - 24*cpi*invgamt23*Sm2 - 24*cpi*invgamt33*Sm3 -
3*ckappa1*trGt3 + 3*ckappa1*trGtd3))/3.
);

Expand Down

0 comments on commit ee6ffbf

Please sign in to comment.