portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Peter's bugfix in ImplicitlyRestartedLanczos.h

Browse Source 
My bugfix in MomentumPhase.hpp
This commit is contained in:
Michael Marshall 2019-10-07 17:36:11 +01:00
parent 803329af99
commit 88d6ff8f1d
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
View File

@ -64,7 +64,7 @@ void basisRotate(std::vector<Field> &basis,Eigen::MatrixXd& Qt,int j0, int j1, i
thread_region thread_region
{ {
std::vector < vobj , commAllocator<vobj> > B(Nm); // Thread private std::vector < vobj > B(Nm); // Thread private
thread_for_in_region(ss, grid->oSites(),{ thread_for_in_region(ss, grid->oSites(),{
for(int j=j0; j<j1; ++j) B[j]=0.; for(int j=j0; j<j1; ++j) B[j]=0.;

2
Hadrons/Modules/MSource/MomentumPhase.hpp
View File

@ -140,7 +140,7 @@ void TMomentumPhase<FImpl>::execute(void)
envGetTmp(LatticeComplex, coor); envGetTmp(LatticeComplex, coor);
p = strToVec<Real>(par().mom); p = strToVec<Real>(par().mom);
ph = zero; ph = Zero();
for(unsigned int mu = 0; mu < env().getNd(); mu++) for(unsigned int mu = 0; mu < env().getNd(); mu++)
{ {
LatticeCoordinate(coor, mu); LatticeCoordinate(coor, mu);