Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
}
if (qep->reason == QEP_CONVERGED_ITERATING) {
if (breakdown) {
/* Stop if breakdown */
PetscInfo2(qep,"Breakdown Quadratic Arnoldi method (it=%i norm=%g)\n",qep->its,beta);
PetscInfo2(qep,"Breakdown Quadratic Arnoldi method (it=%D norm=%G)\n",qep->its,beta);
qep->reason = QEP_DIVERGED_BREAKDOWN;
} else {
/* Prepare the Rayleigh quotient for restart */
for (i=k;i<k+l;i++) {
S[i*qep->ncv+k+l] = Q[(i+1)*nv-1]*beta;