From 6692425aa29f9874534d3881fa8facca27037ac5 Mon Sep 17 00:00:00 2001
From: Chulwoo Jung <chulwoo@ic2submit02.sdcc.bnl.gov>
Date: Wed, 26 Nov 2025 17:17:22 -0500
Subject: [PATCH] Checking in before pulling

---
 Grid/algorithms/iterative/KrylovSchur.h | 6 ++++--
 examples/Example_krylov_schur.cc        | 2 +-
 2 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/Grid/algorithms/iterative/KrylovSchur.h b/Grid/algorithms/iterative/KrylovSchur.h
index 5f8d0f4c..642f6be3 100644
--- a/Grid/algorithms/iterative/KrylovSchur.h
+++ b/Grid/algorithms/iterative/KrylovSchur.h
@@ -320,7 +320,7 @@ class KrylovSchur {
       RealD shift_=1.;
       shift = &shift_;
       if (shift) 
-          std::cout << GridLogMessage << "Shift " << approxLambdaMax << std::endl;
+          std::cout << GridLogMessage << "Shift " << *shift << std::endl;
       MaxIter = _maxIter;
       Nm = _Nm; Nk = _Nk;
       Nstop = _Nstop;
@@ -357,7 +357,7 @@ class KrylovSchur {
 	for (int m=0;m<Nm;m++) temp(m,m) -= *shift;
         Eigen::MatrixXcd RayleighS = temp.inverse();
         Eigen::MatrixXcd temp2 = RayleighS*temp;
-	std::cout << GridLogMessage  << "Shift inverse check: shift= "<<*shift <<std::endl;
+	std::cout << GridLogMessage  << "Shift inverse check: shift= "<<*shift<<" "<< temp2 <<std::endl;
 	
 
         ComplexSchurDecomposition schur (Rayleigh, false, ritzFilter);
@@ -367,6 +367,8 @@ class KrylovSchur {
         // Rearrange Schur matrix so wanted evals are on top left (like MATLAB's ordschur)
         std::cout << GridLogMessage << "Reordering Schur eigenvalues" << std::endl;
         schur.schurReorder(Nk);
+        std::cout << GridLogMessage << "Shifted Schur eigenvalues" << std::endl;
+        schurS.schurReorder(Nk);
         Eigen::MatrixXcd Q = schur.getMatrixQ();
         Qt = Q.adjoint();                           // TODO should Q be real?
         Eigen::MatrixXcd S = schur.getMatrixS();
diff --git a/examples/Example_krylov_schur.cc b/examples/Example_krylov_schur.cc
index c1267779..54e3c9ca 100644
--- a/examples/Example_krylov_schur.cc
+++ b/examples/Example_krylov_schur.cc
@@ -332,7 +332,7 @@ int main (int argc, char ** argv)
   const int Ls=16;
 
 //   GridCartesian         * UGrid   = SpaceTimeGrid::makeFourDimGrid(GridDefaultLatt(), GridDefaultSimd(Nd,vComplex::Nsimd()),GridDefaultMpi());
-  std::vector<int> lat_size {16, 16, 16, 32};
+  std::vector<int> lat_size {16, 16, 16, 16};
   std::cout << "Lattice size: " << lat_size << std::endl;
   GridCartesian * UGrid = SpaceTimeGrid::makeFourDimGrid(lat_size, 
 								          GridDefaultSimd(Nd,vComplex::Nsimd()),