<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"><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 class="">Sharp Stone <<a href="mailto:thronesf@gmail.com">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>-- <br><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div>
</div>