<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jan 20, 2013 at 1:06 PM, Hui Zhang <span dir="ltr"><<a href="mailto:mike.hui.zhang@hotmail.com" target="_blank">mike.hui.zhang@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":2ox">Thanks very much! I am taking a look of this example.  But MATNEST does not support factorization.<br>
I understand now I must do it in a basic way using MatSetValuesLocal.</div></blockquote></div><br>With that example, you run with -dm_mat_type nest if you want to use MATNEST, acknowledging that it means you'll be using a fieldsplit preconditioner. Note that the term "nest" does not appear anywhere in the code. Without that option, an AIJ format will be used, so you can use a direct solver or anything else. MATNEST here is a run-time optimization that is primarily only useful for fieldsplit (or shell) preconditioning. It's a very bad idea to make your code depend on MATNEST because it severely limits your choices of algorithms.</div>
</div>