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;