<div dir="ltr">Each time I call my Jacobian function, it generates the Jacoiban matrix from scratch by inserting entries. The result of this assembling process gives the same nonzero pattern as in the previous iteration. Can SNES/KSP track this too?<div><br></div><div>Thanks.</div><div><br></div><div>Xiangdong </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 9, 2015 at 5:01 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Mat now keeps track of its nonzero pattern. If you do not change the pattern by inserting values then SNES/KSP will automatically keep the same pattern and optimize for it.<br>
<span class="HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
> On Sep 9, 2015, at 3:57 PM, Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> wrote:<br>
><br>
> Hello everyone,<br>
><br>
> As mentioned in the 3.5 changes, MatStructure argument is removed from SNESSetJacobian(). If each time, my Jacobian always has same_nonzero_pattern. How do I pass this information to SNES or the KSP? I only find reuse same preconditioner, but not reuse same pattern.<br>
><br>
> Thank you.<br>
><br>
> Xiangdong<br>
<br>
</div></div></blockquote></div><br></div>