From 3b9fc72451acd01d54213203385b782d06dd438e Mon Sep 17 00:00:00 2001
From: Patrick Oare <patrickoare@gmail.com>
Date: Thu, 14 Aug 2025 14:06:32 -0400
Subject: [PATCH] modified debug output slightly

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

diff --git a/Grid/algorithms/iterative/KrylovSchur.h b/Grid/algorithms/iterative/KrylovSchur.h
index 8cbd8b03..824067c6 100644
--- a/Grid/algorithms/iterative/KrylovSchur.h
+++ b/Grid/algorithms/iterative/KrylovSchur.h
@@ -288,6 +288,8 @@ class KrylovSchur {
 
         std::cout << GridLogMessage << "*** TRUNCATING FOR RESTART *** " << std::endl;
 
+        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);
         // evecs = std::vector<Field> (evecs.begin(), evecs.begin() + Nk);
@@ -298,8 +300,7 @@ class KrylovSchur {
         checkKSDecomposition();
 
         // Compute eigensystem of Rayleigh. Note the eigenvectors correspond to the sorted eigenvalues.
-        std::cout << GridLogDebug << "S: " << std::endl << S << std::endl;
-        // computeEigensystem(S);
+        std::cout << GridLogDebug << "Rayleigh after truncation: " << std::endl << Rayleigh << std::endl;
         computeEigensystem(Rayleigh);
         std::cout << GridLogMessage << "Eigenvalues (first Nk sorted): " << std::endl << evals << std::endl;