| Line 154... |
Line 154... |
if (eps->ispositive) {
|
if (eps->ispositive) {
|
/* Purify eigenvectors */
|
/* Purify eigenvectors */
|
for (i=0;i<eps->nconv;i++) {
|
for (i=0;i<eps->nconv;i++) {
|
ierr = VecCopy(eps->V[i],w);CHKERRQ(ierr);
|
ierr = VecCopy(eps->V[i],w);CHKERRQ(ierr);
|
ierr = STApply(eps->OP,w,eps->V[i]);CHKERRQ(ierr);
|
ierr = STApply(eps->OP,w,eps->V[i]);CHKERRQ(ierr);
|
ierr = IPNorm(eps->ip,eps->V[i],&norm);CHKERRQ(ierr);
|
ierr = VecNormalize(eps->V[i],&norm);CHKERRQ(ierr);
|
ierr = VecScale(eps->V[i],1.0/norm);CHKERRQ(ierr);
|
|
}
|
}
|
}
|
}
|
|
|
/* left eigenvectors */
|
/* left eigenvectors */
|
if (eps->solverclass == EPS_TWO_SIDE) {
|
if (eps->solverclass == EPS_TWO_SIDE) {
|