<div dir="ltr"><div><div>Hi all,<br><br></div>In the ex19, specifically from Lines 581-587, about the "invert the system" part, <br><pre width="80"><a name="line581">581: </a>            <font color="#B22222">/* invert the system:</font>
<a name="line582">582: </a><font color="#B22222">             [ dfu / du     0        0        0    ][yu] = [fu]</font>
<a name="line583">583: </a><font color="#B22222">             [     0    dfv / dv     0        0    ][yv]   [fv]</font>
<a name="line584">584: </a><font color="#B22222">             [ dfo / du dfo / dv dfo / do     0    ][yo]   [fo]</font>
<a name="line585">585: </a><font color="#B22222">             [ dft / du dft / dv     0    dft / dt ][yt]   [ft]</font>
<a name="line586">586: </a><font color="#B22222">             by simple back-substitution</font>
<a name="line587">587: </a><font color="#B22222">           */</font></pre> I realized u and v are independent variables, but even so, the first two rows should be both containing dfu/dv and dfv/du terms just like the 3rd and 4th rows, do I miss some points?<br><br></div>Thank you in advance!<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>