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

bugfix

Browse Source
This commit is contained in:
ferben 2019-02-21 11:13:10 +00:00
parent f70c5b004a
commit 5d6462b706
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Hadrons/Modules/MDistil/Baryon2pt.hpp
View File

@ -97,8 +97,8 @@ void TBaryon2pt<FImpl>::execute(void)
const std::string &inputL{par().inputL};
const std::string &inputR{par().inputR};
const std::string &inputL{par().quarksL};
const std::string &inputR{par().quarksR};
const std::string &quarksL{par().quarksL};
const std::string &quarksR{par().quarksR};
const std::string &output{par().output};
int Nmom=1;
@ -166,8 +166,8 @@ void TBaryon2pt<FImpl>::execute(void)
Eigen::Tensor<Complex,4> B4L = B5L.chip(t,0);
Eigen::Tensor<Complex,4> B4R = B5R.chip(tsrc,0);
for (int is=0 ; is < 4 ; is++){
Eigen::Tensor<Complex,4> B3L = B4L.chip(is,0);
Eigen::Tensor<Complex,4> B3R = B4R.chip(is,0);
Eigen::Tensor<Complex,3> B3L = B4L.chip(is,0);
Eigen::Tensor<Complex,3> B3R = B4R.chip(is,0);
Eigen::Tensor<Complex,0> C2 = B3L.contract(B3R,product_dims);
corr(imom,t) += (double)epsilon_sgn[pairs[ipair]]*C2(0);
}