<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Feb 25, 2014 at 4:12 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> writes:<br>
<br>
> For the MatSetValuesStencil, should it always work in pair with<br>
> DMCreateMatrix?<br>
<br>
</div>Yes, or MatSetStencil.<br>
<div class=""><br>
> However, in some example, e.g., src/ksp/ksp/examples/tutorials/ex25.c, the<br>
> DMCreateMatrix is never called. Is it because KspSetComputeOperators is<br>
> used and DMCreateMatrix is called implicitly somewhere?<br>
<br>
</div>No, it is just testing that certain functionality works even if the DM<br>
is not set.  If you have a DM, we recommend setting it.<br></blockquote><div><br></div><div>I am not clear about this. In this example, where is the format declaration (e.g.,mpiaij) of Mat jac? Why is this mat jac always distributed in a compatible way as the DM vectors?</div>
<div><br></div><div>Thanks.</div><div><br></div><div>Xiangdong</div></div><br></div></div>