In the case of a matrix, shouldn't these checks go the first place we encounter them (MatSetValues)?
Running with -fp_trap is really good at finding these kind of mistakes. On Wed, Jun 13, 2012 at 8:50 AM, Barry Smith <bsmith at mcs.anl.gov> wrote: > > On return from user FormFunction and FormJacobian in debug mode should > we check if all user output is sane? That is numerical values in the matrix > and vector are sane, no 1.e-307 or 1.e85 and stuff like that? > > Barry > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120614/3876f2dd/attachment.html>
