<div dir="ltr">Hi all,<div><br></div><div>I'm using shift-invert with EPS to solve for eigenvalues. I find that if I do only</div><div><br></div><div><font face="monospace">...</font></div><div><font face="monospace">  ierr = EPSGetST(eps,&st);CHKERRQ(ierr);</font></div><div><font face="monospace">  ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);</font></div><div><font face="monospace">...</font></div><div><font face="monospace"><br></font></div><div>in my code I get correct eigenvalues. But if I do </div><div><font face="monospace"><br></font></div><div><font face="monospace">...</font></div><div><div style="font-family:monospace">  ierr = EPSGetST(eps,&st);CHKERRQ(ierr);</div><div style="font-family:monospace">  ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);</div><div style="font-family:monospace">  ierr = STGetKSP(st,&ksp);CHKERRQ(ierr);</div><div style="font-family:monospace">  ierr = KSPGetPC(ksp,&pc);CHKERRQ(ierr);</div><div style="font-family:monospace">  ierr = KSPSetType(ksp,KSPPREONLY);CHKERRQ(ierr);</div><div style="font-family:monospace">  ierr = PCSetType(pc,PCCHOLESKY);CHKERRQ(ierr);</div><div style="font-family:monospace">...</div><div style=""><br></div><div style="">the eigenvalues found by EPS are completely wrong! Somehow I thought I was supposed to do the latter, from the examples etc, but I guess that was not correct? I attach the full piece of test code and a test matrix.</div><div style=""><br></div><div style="">Best,</div><div style="">Greg</div></div></div>