portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-03-21 19:56:09 +00:00
Code Releases Activity

may have found bug

Browse Source
This commit is contained in:
Patrick Oare
2025-08-15 12:13:01 -04:00
parent a18b0d496c
commit 6fd71aea9d
1 changed files with 11 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
Grid/algorithms/iterative/KrylovSchur.h
View File

@@ -290,11 +290,19 @@ class KrylovSchur {
std::cout << GridLogDebug << "Rayleigh before truncation: " << std::endl << Rayleigh << std::endl;
Rayleigh = Rayleigh(Eigen::seqN(0, Nk), Eigen::seqN(0, Nk));
basis = std::vector<Field> (basis.begin(), basis.begin() + Nk);
// Rayleigh = Rayleigh(Eigen::seqN(0, Nk), Eigen::seqN(0, Nk));
Eigen::MatrixXcd RayTmp = Rayleigh(Eigen::seqN(0, Nk), Eigen::seqN(0, Nk));
Rayleigh = RayTmp;
// basis = std::vector<Field> (basis.begin(), basis.begin() + Nk);
std::vector<Field> basisTmp = std::vector<Field> (basis.begin(), basis.begin() + Nk);
basis = basisTmp;
// evecs = std::vector<Field> (evecs.begin(), evecs.begin() + Nk);
// littleEvecs = littleEvecs(Eigen::seqN(0, Nk), Eigen::seqN(0, Nk));
Eigen::VectorXcd btmp = b.head(Nk); // careful about how you slice Eigen::VectorXcd!
// b = b.head(Nk);
Eigen::VectorXcd btmp = b.head(Nk);
b = btmp;
std::cout << GridLogDebug << "Rayleigh after truncation: " << std::endl << Rayleigh << std::endl;