[petsc-users] help needed for -snes_mf_operator

Matthew Knepley knepley at gmail.com
Thu Apr 26 13:02:14 CDT 2012


On Thu, Apr 26, 2012 at 1:14 PM, Dominik Szczerba <domel07 at gmail.com> wrote:

> > Thanks for a useful hint. The smallest problem I was able to generate
> > has around 1000 unknowns but it finally outputs something. I got some
> > matrix, with the biggest absolute entry order 0.1 the smallest
> > absolute entry order 1e-14, average values seem to be order 1e-8 or
> > so. End line of the output says: Norm of matrix ratio 0.579314
> > difference 10.9951 (constant state 1.0)
> >
> > Dominik
>
> So do you think this is reasonable or you still think I evaluate the
> Jacobian incorrectly? Is there still something I can do to find out where
> the problem is? Unfortunately, there is no easy way I can run the code
> serially.
>
Yes, it looks incorrect. The difference should be very small.

   Matt


> Dominik
>
-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120426/293b20ea/attachment.htm>


More information about the petsc-users mailing list