| Line 479... |
Line 479... |
#endif
|
#endif
|
if (info) SETERRQ1(PETSC_ERR_LIB,"Error in Lapack xTREVC %i",info);
|
if (info) SETERRQ1(PETSC_ERR_LIB,"Error in Lapack xTREVC %i",info);
|
}
|
}
|
|
|
/* pX[i] <- pX[i] / ||pX[i]|| */
|
/* pX[i] <- pX[i] / ||pX[i]|| */
|
ierr = SlepcDenseNorm(pX, d->nconv, d->nconv, d->nconv); CHKERRQ(ierr);
|
ierr = SlepcDenseNorm(pX, d->nconv, d->nconv, d->nconv, d->eigi);
|
|
CHKERRQ(ierr);
|
|
|
/* V <- cX * pX */
|
/* V <- cX * pX */
|
ierr = SlepcUpdateVectorsZ(eps->V, 0.0, 1.0, d->cX, d->size_cX, pX,
|
ierr = SlepcUpdateVectorsZ(eps->V, 0.0, 1.0, d->cX, d->size_cX, pX,
|
d->nconv, d->nconv, d->nconv); CHKERRQ(ierr);
|
d->nconv, d->nconv, d->nconv); CHKERRQ(ierr);
|
|
|