Subversion Repositories slepc-dev

Rev

Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 2370 Rev 2394
Line 274... Line 274...
    }
    }
 
 
    if (qep->reason == QEP_CONVERGED_ITERATING) {
    if (qep->reason == QEP_CONVERGED_ITERATING) {
      if (breakdown) {
      if (breakdown) {
        /* Stop 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;
        qep->reason = QEP_DIVERGED_BREAKDOWN;
      } else {
      } else {
        /* Prepare the Rayleigh quotient for restart */
        /* Prepare the Rayleigh quotient for restart */
        for (i=k;i<k+l;i++) {
        for (i=k;i<k+l;i++) {
          S[i*qep->ncv+k+l] = Q[(i+1)*nv-1]*beta;
          S[i*qep->ncv+k+l] = Q[(i+1)*nv-1]*beta;