<div dir="ltr"><div>Hi all,<br><br></div><div>I'm sorry I'm still confusing that if I have calculated the Jacobian elements, how can I go ahead to solve the linear system with dof>1 at each node? The example  code ex19 solves the linear system with Guass-Seidel method, and now I do want to use Jacobian matrix because I've found the Jacobian elements. Any ideas? Many thanks!<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 2, 2014 at 1:26 PM, Sharp Stone <span dir="ltr"><<a href="mailto:thronesf@gmail.com" target="_blank">thronesf@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matt and Jed,<div><br></div><div>Thank you very much for your replies. In source code ex19, the routine SNESSetNGS() is used, and no other options or Jacobi routines are used. How can we change the iteration method? The routine SNESSetFromOptions(snes) could do this for the user? If I understand correctly, petsc would use Gauss-Seidel or Jacobi method as designated with even the same source code?<br></div><div><br></div><div>Thank you in advance!</div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Thu, Oct 2, 2014 at 1: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"><span>Sharp Stone <<a href="mailto:thronesf@gmail.com" target="_blank">thronesf@gmail.com</a>> writes:<br>
<br>
> Matt,<br>
><br>
> Thank you very much!<br>
><br>
> The ex19 used the Gauss-Seidel method. But I'm a little bit wonder whether<br>
> we have any examples of Jacobi method to solve such linear systems "(mat)A<br>
> dot (vec)X = (vec)B" with dof>1? Thanks in advance!<br>
<br>
</span>ex19 does whatever you want: -pc_type jacobi and the like.<br>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div>
</div>