<p>&gt; Thanks for a useful hint. The smallest problem I was able to generate<br>
&gt; has around 1000 unknowns but it finally outputs something. I got some<br>
&gt; matrix, with the biggest absolute entry order 0.1 the smallest<br>
&gt; absolute entry order 1e-14, average values seem to be order 1e-8 or<br>
&gt; so. End line of the output says: Norm of matrix ratio 0.579314<br>
&gt; difference 10.9951 (constant state 1.0)<br>
&gt;<br>
&gt; Dominik</p>
<p>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.</p>

<p>Dominik<br>
</p>