<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Here is ts_view.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 2, 2015 at 2:40 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
    Run for one time-step with -ts_view so we can see exactly what solver options you are using.<br>
<br>
    Then run with the additional options -ksp_pc_side right -ts_max_snes_failures -1<br>
<br>
<br>
> On Jun 2, 2015, at 12:26 PM, Italo Tasso <<a href="mailto:italo@tasso.com.br">italo@tasso.com.br</a>> wrote:<br>
><br>
> I made a code to solve the Navier-Stokes equations, incompressible, non-linear, all coupled, finite differences, staggered grid.<br>
><br>
> I am running the code with:<br>
><br>
> -ts_monitor -snes_monitor -ksp_monitor_true_residual -snes_converged_reason -ksp_converged_reason -pc_type fieldsplit -pc_fieldsplit_type schur -pc_fieldsplit_detect_saddle_point<br>
><br>
> It works very well most of the time. But in some cases, the solver halts for a long time then KSP does not converge.<br>
><br>
> See output1.txt. It seems that the residual is already very small, close to machine zero, but KSP doesn't stop.<br>
><br>
> So I added -ksp_atol 1e-10. See output2.txt. Now it fails on a different time step.<br>
><br>
> I also tried -ksp_norm_type unpreconditioned. It works for this case (grid size), but fail for other cases.<br>
><br>
> I also tried building the Jacobian and including null space. It fixes some cases but causes others that worked before to fail. Seems really random.<br>
><br>
> It feels like this is related to the PC, because the code halts for a long time at the first KSP step, then diverges.<br>
><br>
> Any suggestions?<br>
> <output1.txt><output2.txt><br>
<br>
</blockquote></div><br></div></div>