| Line 189... |
Line 189... |
/* Reduce S to Hessenberg form, S <- Q S Q' */
|
/* Reduce S to Hessenberg form, S <- Q S Q' */
|
ierr = EPSDenseHessenberg(n,eps->nconv,S,lds,Q);CHKERRQ(ierr);
|
ierr = EPSDenseHessenberg(n,eps->nconv,S,lds,Q);CHKERRQ(ierr);
|
/* Reduce S to (quasi-)triangular form, S <- Q S Q' */
|
/* Reduce S to (quasi-)triangular form, S <- Q S Q' */
|
ierr = EPSDenseSchur(n,eps->nconv,S,lds,Q,eps->eigr,eps->eigi);CHKERRQ(ierr);
|
ierr = EPSDenseSchur(n,eps->nconv,S,lds,Q,eps->eigr,eps->eigi);CHKERRQ(ierr);
|
/* Sort the remaining columns of the Schur form */
|
/* Sort the remaining columns of the Schur form */
|
ierr = EPSSortDenseSchur(eps,n,eps->nconv,S,PETSC_NULL,lds,Q,PETSC_NULL,eps->eigr,eps->eigi);CHKERRQ(ierr);
|
ierr = EPSSortDenseSchur(eps,n,eps->nconv,S,lds,Q,eps->eigr,eps->eigi);CHKERRQ(ierr);
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "EPSSolve_KRYLOVSCHUR_HARMONIC"
|
#define __FUNCT__ "EPSSolve_KRYLOVSCHUR_HARMONIC"
|