[petsc-users] EPSGetConverge problem
Eleni Siampli
siaeleni at hotmail.com
Thu Aug 23 17:30:17 CDT 2012
Hello,
I am new to petsc and I am trying to calculate the eigenpairs.
My code is like the following:
ierr = EPSCreate(PETSC_COMM_WORLD, &eps); CHKERRQ(ierr);
ierr = EPSSetOperators(eps, S, D); CHKERRQ(ierr);
ierr = EPSSetProblemType(eps, EPS_GHEP); CHKERRQ(ierr);
PetscInt nev = 300;
ierr = EPSSetDimensions(eps, nev, PETSC_DECIDE, PETSC_DECIDE); CHKERRQ(ierr);
ierr = EPSSetWhichEigenpairs(eps, EPS_SMALLEST_MAGNITUDE);CHKERRQ(ierr);
ierr = EPSSolve(eps); CHKERRQ(ierr);
ierr = EPSGetConverged(eps, &m_nconv);CHKERRQ(ierr);
What I don't understand is why m_conv is greater than nev.
when I try to print the m_conv out:
PetscPrintf(PETSC_COMM_WORLD," Number of converged eigenpairs: %D\n\n",m_nconv);
it gives me 399 instead of 300 (that I want).
Thank you in advance,
Eleni
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120824/988ac54e/attachment.html>
More information about the petsc-users
mailing list