<div dir="ltr">Mark, the problem was that the Inode structure was not learning about the diagonal change in MatScale(). Fixed now:<div><br></div><div><a href="https://bitbucket.org/petsc/petsc-dev/commits/44a71c6d4b84406696d7c21549916099ada394b8">https://bitbucket.org/petsc/petsc-dev/commits/44a71c6d4b84406696d7c21549916099ada394b8</a><br>
<br>I also modified the example to get rid of all the preprocessor crud with the same functionality (now at run-time instead of compile time).<br><br><a href="https://bitbucket.org/petsc/petsc-dev/commits/a4ec749136a25a44f7afdd293fac7c273adfaed0">https://bitbucket.org/petsc/petsc-dev/commits/a4ec749136a25a44f7afdd293fac7c273adfaed0</a><br>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Dec 18, 2012 at 4:02 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I'm getting a floating point exception in the eigen estimator for cheby on the 2nd solve of ex56 (you need to uncomment the define of TwoSolve) when I use sor PC (runex56 uses jacobi so you need to change that).  I see this on Cray and my Mac.<br>

<br>
I also found that the eigen estimator was called before the first solve, and worked fine, but I set the eigen values so I don't want it computing eigen values.<br>
<br>
I can understand cheby being the eager beaver and recomputing eigen values when I set a new matrix, but this not the semantics that I prefer.<br>
<br>
So there seem to be three issues:<br>
<br>
1) seems the floating point exception.  This only happens with a non-debug build.  I'm also using 64 bit indices.<br>
<br>
2) computing eigenvalues even though I set them.<br>
<br>
3) I would like to tell cheby not to recompute eigenvalues.<br>
<span class="HOEnZb"><font color="#888888"><br>
Mark<br>
<br>
<br>
</font></span></blockquote></div><br></div>