| Line 271... |
Line 271... |
|
|
/* 5. Sort diagonal elements in T and accumulate rotations on U */
|
/* 5. Sort diagonal elements in T and accumulate rotations on U */
|
ierr = EPSSortDenseSchur(nv,eps->nconv,T,ncv,U,eps->eigr,eps->eigi,eps->which);CHKERRQ(ierr);
|
ierr = EPSSortDenseSchur(nv,eps->nconv,T,ncv,U,eps->eigr,eps->eigi,eps->which);CHKERRQ(ierr);
|
|
|
/* 6. AV(:,idx) = AV * U(:,idx) */
|
/* 6. AV(:,idx) = AV * U(:,idx) */
|
ierr = EPSUpdateVectors(nv,eps->AV,eps->nconv,nv,U,nv,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SlepcUpdateVectors(nv,eps->AV,eps->nconv,nv,U,nv,PETSC_FALSE);CHKERRQ(ierr);
|
|
|
/* 7. V(:,idx) = V * U(:,idx) */
|
/* 7. V(:,idx) = V * U(:,idx) */
|
ierr = EPSUpdateVectors(nv,eps->V,eps->nconv,nv,U,nv,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SlepcUpdateVectors(nv,eps->V,eps->nconv,nv,U,nv,PETSC_FALSE);CHKERRQ(ierr);
|
|
|
/* Compute residuals */
|
/* Compute residuals */
|
for (i=0;i<nv;i++) { rsdold[i] = rsd[i]; }
|
for (i=0;i<nv;i++) { rsdold[i] = rsd[i]; }
|
|
|
ierr = EPSSchurResidualNorms(eps,eps->V,eps->AV,T,eps->nconv,nv,ncv,rsd);CHKERRQ(ierr);
|
ierr = EPSSchurResidualNorms(eps,eps->V,eps->AV,T,eps->nconv,nv,ncv,rsd);CHKERRQ(ierr);
|