Mark, the problem was that the Inode structure was not learning about the
diagonal change in MatScale(). Fixed now:
https://bitbucket.org/petsc/petsc-dev/commits/44a71c6d4b84406696d7c21549916099ada394b8
I also modified the example to get rid of all the preprocessor crud with
the same functionalit
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.
I also found that the eigen estimator was called be