portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-19 00:07:05 +01:00
Code Releases Activity

thread_for caused the problems - slow for loop for now

Browse Source
This commit is contained in:
Felix Erben
2019-09-26 13:33:49 +01:00
parent 8415e23fc6
commit bf62ec163d
2 changed files with 11 additions and 38 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -96,10 +96,9 @@ void BaryonUtils<FImpl>::ContractBaryons(const PropagatorField &q1_src,
auto v2 = q2_src.View();
auto v3 = q3_src.View();
std::cout << "wick contract " << wick_contraction << std::endl;
int count=0;
// accelerator_for(ss, grid->oSites(), grid->Nsimd(), {
thread_for(ss,grid->oSites(),{
//thread_for(ss,grid->oSites(),{
for(int ss=0; ss < grid->oSites(); ss++){
auto D1 = v1[ss];
auto D2 = v2[ss];
@ -111,10 +110,6 @@ void BaryonUtils<FImpl>::ContractBaryons(const PropagatorField &q1_src,
auto gD3 = GammaB * D3;
vobj result=Zero();
/* if (count<10){
std::cout << "outside epsilon " << count << std::endl;
count++;
}*/
for (int ie_src=0; ie_src < 6 ; ie_src++){
int a_src = epsilon[ie_src][0]; //a
@ -182,6 +177,6 @@ void BaryonUtils<FImpl>::ContractBaryons(const PropagatorField &q1_src,
}
vbaryon_corr[ss] = result;
} );//end loop over lattice sites
} // );//end loop over lattice sites
}
NAMESPACE_END(Grid);