> 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. Dominik -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120426/7c438bb8/attachment.htm>
