From 722a4a723eee4e7a10aba98b175fc32fed9817b5 Mon Sep 17 00:00:00 2001
From: Heather M  Switzer <hmswitzer@email.wm.edu>
Date: Fri, 22 Mar 2024 10:26:14 -0400
Subject: [PATCH] Trying to see why Femto is producing different results

---
 src/eigs/convergence.c      | 1 +
 src/eigs/solve_projection.c | 4 ++--
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/eigs/convergence.c b/src/eigs/convergence.c
index c1ce35fa..63b51171 100644
--- a/src/eigs/convergence.c
+++ b/src/eigs/convergence.c
@@ -149,6 +149,7 @@ int check_convergence_Sprimme(SCALAR *X, PRIMME_INT ldX, int givenX, SCALAR *R,
       CHKERR(convTestFun_Sprimme(hVals[i], X ? &X[ldX * (i - left)] : NULL,
             givenX, blockNorms[i - left], &isConv, ctx));
 
+      
       if (isConv) {
          flags[i] = CONVERGED;
       }
diff --git a/src/eigs/solve_projection.c b/src/eigs/solve_projection.c
index 3bc41a16..7cb6a453 100644
--- a/src/eigs/solve_projection.c
+++ b/src/eigs/solve_projection.c
@@ -1024,7 +1024,7 @@ int map_vecs_Sprimme(HSCALAR *V, int m, int nV, int ldV, HSCALAR *W, int n0,
       Vnorms[i] = sqrt(REAL_PART(
             Num_dot_SHprimme(m, &V[ldV * i], 1, &V[ldV * i], 1, ctx)));
    }
-      
+
    /* Compute V'*W[n0:n-1] */
 
    HSCALAR *ip = NULL;
@@ -1056,7 +1056,7 @@ int map_vecs_Sprimme(HSCALAR *V, int m, int nV, int ldV, HSCALAR *W, int n0,
       if (jmax < 0) {
          jmax = i;
       }
-
+      //printf("i = %d, jmax = %d, ipmax = %lf\n", i, jmax, ipmax);
       /* Assign the map */
 
       p[i] = jmax;