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>
