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

added another debug feature

Browse Source
This commit is contained in:
Patrick Oare
2025-08-14 14:37:49 -04:00
parent 3b9fc72451
commit 16d3c9cf75
1 changed files with 3 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
Grid/algorithms/iterative/KrylovSchur.h
+3
View File
@@ -519,6 +519,7 @@ class KrylovSchur {
int k = basis.size(); // number of basis vectors, also the size of Rayleigh. int k = basis.size(); // number of basis vectors, also the size of Rayleigh.
// rotate basis by Rayleigh to construct UR // rotate basis by Rayleigh to construct UR
// std::vector<Field> rotated;
std::vector<Field> rotated = basis; std::vector<Field> rotated = basis;
Eigen::MatrixXcd Rt = Rayleigh.adjoint(); Eigen::MatrixXcd Rt = Rayleigh.adjoint();
constructUR(rotated, basis, Rayleigh, k); // manually rotate constructUR(rotated, basis, Rayleigh, k); // manually rotate
@@ -611,8 +612,10 @@ class KrylovSchur {
for (int i = 0; i < N; i++) { for (int i = 0; i < N; i++) {
tmp = Zero(); tmp = Zero();
for (int j = 0; j < N; j++) { for (int j = 0; j < N; j++) {
std::cout << GridLogDebug << "Adding R(j, i) = " << R(j, i) << " to rotated" << std::endl;
tmp = tmp + U[j] * R(j, i); tmp = tmp + U[j] * R(j, i);
} }
std::cout << GridLogDebug << "rotated norm at i = " << i << " is: " << norm2(tmp) << std::endl;
UR.push_back(tmp); UR.push_back(tmp);
} }
return; return;