[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