From 16d3c9cf7580e1fbad5a7099c1285911c42fedbc Mon Sep 17 00:00:00 2001
From: Patrick Oare <patrickoare@gmail.com>
Date: Thu, 14 Aug 2025 14:37:49 -0400
Subject: [PATCH] added another debug feature

---
 Grid/algorithms/iterative/KrylovSchur.h | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/Grid/algorithms/iterative/KrylovSchur.h b/Grid/algorithms/iterative/KrylovSchur.h
index 824067c6..8750ee65 100644
--- a/Grid/algorithms/iterative/KrylovSchur.h
+++ b/Grid/algorithms/iterative/KrylovSchur.h
@@ -519,6 +519,7 @@ class KrylovSchur {
       int k = basis.size();         // number of basis vectors, also the size of Rayleigh.
 
       // rotate basis by Rayleigh to construct UR
+      // std::vector<Field> rotated;
       std::vector<Field> rotated = basis;
       Eigen::MatrixXcd Rt = Rayleigh.adjoint();
       constructUR(rotated, basis, Rayleigh, k);             // manually rotate
@@ -611,8 +612,10 @@ class KrylovSchur {
       for (int i = 0; i < N; i++) {
         tmp = Zero();
         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);
         }
+        std::cout << GridLogDebug << "rotated norm at i = " << i << " is: " << norm2(tmp) << std::endl;
         UR.push_back(tmp);
       }
       return;