Subversion Repositories slepc-dev

Rev

Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 2797 Rev 2804
Line 205... Line 205... 
    ierr = PSSort(eps->ps,eps->eigr,eps->eigi,eps->which_func,eps->which_ctx);CHKERRQ(ierr);
 
    ierr = PSSort(eps->ps,eps->eigr,eps->eigi,eps->which_func,eps->which_ctx);CHKERRQ(ierr);
 
 
 
 
 
    /* Check convergence */
 
    /* Check convergence */
 
    ierr = PSGetArray(eps->ps,PS_MAT_A,&S);CHKERRQ(ierr);
 
    ierr = PSGetArray(eps->ps,PS_MAT_A,&S);CHKERRQ(ierr);
 
    ierr = PSGetArray(eps->ps,PS_MAT_Q,&Q);CHKERRQ(ierr);
 
    ierr = PSGetArray(eps->ps,PS_MAT_Q,&Q);CHKERRQ(ierr);
 
    ierr = EPSKrylovConvergence(eps,PETSC_FALSE,eps->trackall,eps->nconv,nv-eps->nconv,S,ld,Q,ld,eps->V,nv,beta,gamma,&k,work);CHKERRQ(ierr);
 
    ierr = EPSKrylovConvergence(eps,PETSC_FALSE,PETSC_FALSE,eps->nconv,nv-eps->nconv,S,ld,Q,ld,eps->V,nv,beta,gamma,&k,work);CHKERRQ(ierr);
 
    ierr = PSRestoreArray(eps->ps,PS_MAT_A,&S);CHKERRQ(ierr);
 
    ierr = PSRestoreArray(eps->ps,PS_MAT_A,&S);CHKERRQ(ierr);
 
    ierr = PSRestoreArray(eps->ps,PS_MAT_Q,&Q);CHKERRQ(ierr);
 
    ierr = PSRestoreArray(eps->ps,PS_MAT_Q,&Q);CHKERRQ(ierr);
 
    if (eps->its >= eps->max_it) eps->reason = EPS_DIVERGED_ITS;
 
    if (eps->its >= eps->max_it) eps->reason = EPS_DIVERGED_ITS;
 
    if (k >= eps->nev) eps->reason = EPS_CONVERGED_TOL;
 
    if (k >= eps->nev) eps->reason = EPS_CONVERGED_TOL;