<br><br><div class="gmail_quote">On Mon, Nov 7, 2011 at 03:00, Konstantinos Kontzialis <span dir="ltr">&lt;<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
<br>
    ierr = SNESSetFromOptions(sys.snes);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = SNESGetKSP(sys.snes, &amp;sys.ksp2);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = KSPGetPC(sys.ksp2, &amp;sys.pc);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = KSPSetFromOptions(sys.ksp2);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = TSSetFromOptions(sys.ts);<br>
    CHKERRQ(ierr);<br></blockquote><div><br></div><div>You can just call TSSetFromOptions, it automatically calls *SetFromOptions on all the inner solvers.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<br>
    ierr = TSSolve(sys.ts, sys.gsv, &amp;sys.con-&gt;etime);<br>
    CHKERRQ(ierr);<br>
<br>
and I run with:<br>
<br>
mpiexec -n 8 ./hoac cylinder -llf_flux -n_out 5 -end_time 1 -implicit -implicit_type 3 -pc_type asm -sub_pc_type ilu -sub_pc_factor_mat_ordering_<u></u>type rcm 1.0e-8 -gl -sub_ksp_type fgmres -ksp_rtol 1.0e-8 -sub_pc_factor_levels 2 -dt 1.0e-1 -snes_monitor -ksp_pc_side right -snes_converged_reason -ts_type beuler -ksp_converged_reason -ts_view<br>

<br><br></blockquote><div><br></div><div>[...]</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
I observe that the solution is not updated. What am I doing wrong?<br></blockquote><div><br></div><div>From the convergence history, your &quot;matrix&quot; might be diagonal. What is the physical time scale for this problem? Perhaps your base_residual_implicit is incorrect.</div>
</div>