On Tue, 29 Jun 2010 12:43:36 -0500, Matthew Knepley <knepley at gmail.com> wrote: > 3) What you really want is to make MatScalar to be float. Then use > doubles for the residual calculation.
I have no idea what Luke's usage is, but MatScalar != PetscScalar is not currently supported. http://lists.mcs.anl.gov/pipermail/petsc-dev/2009-June/001417.html I'm intrigued by it because my matrices are only for preconditioning purposes, but haven't caught the bug to try fixing it since the above thread. Jed
