portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Baryon revert sign

Browse Source
This commit is contained in:
Raoul Hodgson 2020-06-26 10:08:42 +01:00
parent 102089798c
commit 77af9a3ddc
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Grid/qcd/utils/BaryonUtils.h
View File

@ -257,7 +257,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
auto GAf_D1_GAi_P_rr_cc = GAf_D1_GAi_P()(rho,rho)(c_f,c_i);
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() -= ee * GAf_D1_GAi_P_rr_cc
result()()() += ee * GAf_D1_GAi_P_rr_cc
* D2_GBi()(alpha_f,beta_i)(a_f,a_i)
* GBf_D3 ()(alpha_f,beta_i)(b_f,b_i);
}}
@ -269,7 +269,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
auto D1_GAi_P_ar_ac = D1_GAi_P()(alpha_f,rho)(a_f,c_i);
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() -= ee * D1_GAi_P_ar_ac
result()()() += ee * D1_GAi_P_ar_ac
* GBf_D2_GBi ()(alpha_f,beta_i) (b_f,a_i)
* GAf_D3 ()(rho,beta_i) (c_f,b_i);
}
@ -281,7 +281,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
auto GBf_D1_GAi_P_ar_bc = GBf_D1_GAi_P()(alpha_f,rho)(b_f,c_i);
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() -= ee * GBf_D1_GAi_P_ar_bc
result()()() += ee * GBf_D1_GAi_P_ar_bc
* GAf_D2_GBi ()(rho,beta_i) (c_f,a_i)
* D3 ()(alpha_f,beta_i) (a_f,b_i);
}
@ -293,7 +293,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
auto GAf_D1_GAi_P_rr_cc = GAf_D1_GAi_P()(rho,rho)(c_f,c_i);
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() += ee * GAf_D1_GAi_P_rr_cc
result()()() -= ee * GAf_D1_GAi_P_rr_cc
* GBf_D2_GBi ()(alpha_f,beta_i) (b_f,a_i)
* D3 ()(alpha_f,beta_i) (a_f,b_i);
}
@ -305,7 +305,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
auto GBf_D1_GAi_P_ar_bc = GBf_D1_GAi_P()(alpha_f,rho)(b_f,c_i);
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() += ee * GBf_D1_GAi_P_ar_bc
result()()() -= ee * GBf_D1_GAi_P_ar_bc
* D2_GBi ()(alpha_f,beta_i) (a_f,a_i)
* GAf_D3 ()(rho,beta_i) (c_f,b_i);
}
@ -317,7 +317,7 @@ void BaryonUtils<FImpl>::baryon_site(const mobj &D1,
for (int alpha_f=0; alpha_f<Ns; alpha_f++){
auto D1_GAi_P_ar_ac = D1_GAi_P()(alpha_f,rho)(a_f,c_i);
for (int beta_i=0; beta_i<Ns; beta_i++){
result()()() += ee * D1_GAi_P_ar_ac
result()()() -= ee * D1_GAi_P_ar_ac
* GAf_D2_GBi ()(rho,beta_i) (c_f,a_i)
* GBf_D3 ()(alpha_f,beta_i) (b_f,b_i);
}