portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-05-21 17:44:16 +01:00
Code Releases Activity

Merge pull request #470 from paboyle/specflow

Browse Source 
Spectral flow, DWF/Mobius kernel measurement
This commit is contained in:
Peter Boyle
2025-04-24 14:34:33 -04:00
committed by GitHub
parent b75809ed61 a957e7bfa1
commit 6a1c64fbdd
5 changed files with 853 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
+4 -2
View File
@@ -245,9 +245,10 @@ until convergence
_HermOp(src_n,tmp);
// std::cout << GridLogMessage<< tmp<<std::endl; exit(0);
// std::cout << GridLogIRL << " _HermOp " << norm2(tmp) << std::endl;
RealD vnum = real(innerProduct(src_n,tmp)); // HermOp.
// RealD vnum = real(innerProduct(src_n,tmp)); // HermOp.
RealD vnum = real(innerProduct(tmp,tmp)); // HermOp^2.
RealD vden = norm2(src_n);
RealD na = vnum/vden;
RealD na = std::sqrt(vnum/vden);
if (fabs(evalMaxApprox/na - 1.0) < 0.0001)
i=_MAX_ITER_IRL_MEVAPP_;
evalMaxApprox = na;
@@ -255,6 +256,7 @@ until convergence
src_n = tmp;
}
}
std::cout << GridLogIRL << " Final evalMaxApprox " << evalMaxApprox << std::endl;
std::vector<RealD> lme(Nm);
std::vector<RealD> lme2(Nm);