portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-11 14:40:46 +01:00
Code Releases Activity

started on baryon flavour sums

Browse Source
This commit is contained in:
ferben 2019-02-15 15:51:46 +00:00
parent 668b1e77c7
commit bfd2770657
1 changed files with 18 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
Hadrons/Modules/MDistil/BContraction.hpp
View File

@ -178,7 +178,24 @@ void TBContraction<FImpl>::execute(void)
BaryonTensorSet BField3(Nmom,4,Nt,N_1,N_2,N_3); BaryonTensorSet BField3(Nmom,4,Nt,N_1,N_2,N_3);
Eigen::Tensor<Complex, 3> corr(Nmom,4,Nt); Eigen::Tensor<Complex, 3> corr(Nmom,4,Nt);
//Needs more work - but this is important for contraction
/* int Npairs = 0;
char left[] = "uud";
char right[] = "uud";
std::vector<std::vector<int>> pairs;
for (int il=0, i=0 ; il < 3 ; il++){
for (int ir=0 ; ir < 3 ; ir++){
if (il>ir) continue;
if (left[il] != right[il]) continue;
pairs[i][0]=il;
pairs[i][1]=ir;
i++;
Npairs = i;
}
std::cout << "pairs: " << pairs << std::endl;
}
*/
Complex diquark2; Complex diquark2;
for (int i1=0 ; i1 < N_1 ; i1++){ for (int i1=0 ; i1 < N_1 ; i1++){
for (int i2=0 ; i2 < N_2 ; i2++){ for (int i2=0 ; i2 < N_2 ; i2++){