<div dir="ltr"><div>Hi Jed,<br><br></div>Thank you very much!<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 2, 2014 at 3:02 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"><span class="">Sharp Stone <<a href="mailto:thronesf@gmail.com">thronesf@gmail.com</a>> writes:<br>
<br>
> Hi all,<br>
><br>
> I'm sorry I'm still confusing that if I have calculated the Jacobian<br>
> elements, how can I go ahead to solve the linear system with dof>1 at each<br>
> node? The example code ex19 solves the linear system with Guass-Seidel<br>
> method,<br>
<br>
</span>As we just said, the same code (ex19.c) can use many different methods.<br>
The default method is Newton-Krylov with block Jacobi preconditioning.<br>
Run with -snes_view to see details.<br>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div>
</div>