Trying to see why Femto is producing different results

primme · Mar 22, 2024 · 722a4a7 · 722a4a7
1 parent 97a5096
commit 722a4a7
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git 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
@@ -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;