<div dir="ltr"><div><div><div><div>I originally solve that example problem using LU. But when I solve this one:<br><br><a href="http://fenicsproject.org/documentation/dolfin/1.5.0/python/demo/documented/stokes-iterative/python/documentation.html">http://fenicsproject.org/documentation/dolfin/1.5.0/python/demo/documented/stokes-iterative/python/documentation.html</a><br><br></div>By simply running their code as is for TH and adding the one like I mentioned for MTH, I get the following outputs when I pass in -ksp_monitor -ksp_view and -log_summary<br><br></div>The latter obviously takes a greater amount of time and iterations to converge, and it was using the solver and precondition options that was originally designed for P2/P1. I haven't experimented around with this fully yet.<br><br></div>Thanks,<br></div>Justin<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 2, 2015 at 1:19 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="">Lawrence Mitchell <<a href="mailto:lawrence.mitchell@imperial.ac.uk">lawrence.mitchell@imperial.ac.uk</a>> writes:<br>
> Maybe Justin can chime in here, I don't know, I just happened to know<br>
> how the fenics implementation produces the "basis", so proffered that.<br>
<br>
</span>Thanks, Lawrence.  Unfortunately, my original questions remain<br>
unanswered and now I'm doubly curious why FEniCS appears not to fail due<br>
to the singular linear system.<br>
</blockquote></div><br></div>