[petsc-users] How to get the eigenvectors in Slepc

[Hong Zhang]

Sonya :
Test slpec example first, e.g.,
slepc-3.3-p3/src/eps/examples/tutorials>grep GetEigenvector *.c
ex12.c:      ierr = EPSGetEigenvector(eps,i,X[i],PETSC_NULL);CHKERRQ(ierr);
ex12.c:      ierr =
EPSGetEigenvectorLeft(eps,i,Y[i],PETSC_NULL);CHKERRQ(ierr);
ex7.c:      ierr = EPSGetEigenvector(eps,i,xr,xi);CHKERRQ(ierr);

i.e., ex7.c and ex12.c call EPSGetEigenvector().

Hong

Dear all,
>
> I experience a difficulty with retrieving the eigenvectors
> obtained after eigen solution of generalized eigenvalue A*X=lambda*B*X
> problem.The Slepc nicely finds the eigenvalues but doesn't retrieve the
> eigenvectors after calling the EPSGetEigenvector(eps,1,vec,vec1). I also
> tried the EPSGetEigenpairs function but neither of them
>
> I attach the main program codes.
>
> Your help will be appreciated.
>
> Regards,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130413/94c39790/attachment.html>