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

Merge branch 'feature/gparity_HMC' of https://github.com/paboyle/Grid into feature/gparity_HMC

Browse Source
This commit is contained in:
Peter Boyle 2021-05-04 13:45:11 -04:00
commit 7b02acb2bd
3 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Grid/algorithms/iterative/ConjugateGradientMultiShift.h
View File

@ -182,6 +182,9 @@ public:
for(int s=0;s<nshift;s++) { for(int s=0;s<nshift;s++) {
axpby(psi[s],0.,-bs[s]*alpha[s],src,src); axpby(psi[s],0.,-bs[s]*alpha[s],src,src);
} }
std::cout << GridLogIterative << "ConjugateGradientMultiShift: initial rn (|src|^2) =" << rn << " qq (|MdagM src|^2) =" << qq << " d ( dot(src, [MdagM + m_0]src) ) =" << d << " c=" << c << std::endl;
/////////////////////////////////////// ///////////////////////////////////////
// Timers // Timers

2
HMC/DWF2p1fIwasakiGparity.cc
View File

@ -354,7 +354,7 @@ int main(int argc, char **argv) {
//Setup the BCs //Setup the BCs
FermionActionD::ImplParams Params; FermionActionD::ImplParams Params;
for(int i=0;i<Nd-1;i++) Params.twists = user_params.GparityDirs[i]; //G-parity directions for(int i=0;i<Nd-1;i++) Params.twists[i] = user_params.GparityDirs[i]; //G-parity directions
Params.twists[Nd-1] = 1; //APBC in time direction Params.twists[Nd-1] = 1; //APBC in time direction
std::vector<int> dirs4(Nd); std::vector<int> dirs4(Nd);

2
HMC/DWF2p1fIwasakiGparityRHMCdouble.cc
View File

@ -354,7 +354,7 @@ int main(int argc, char **argv) {
//Setup the BCs //Setup the BCs
FermionActionD::ImplParams Params; FermionActionD::ImplParams Params;
for(int i=0;i<Nd-1;i++) Params.twists = user_params.GparityDirs[i]; //G-parity directions for(int i=0;i<Nd-1;i++) Params.twists[i] = user_params.GparityDirs[i]; //G-parity directions
Params.twists[Nd-1] = 1; //APBC in time direction Params.twists[Nd-1] = 1; //APBC in time direction
std::vector<int> dirs4(Nd); std::vector<int> dirs4(Nd);