<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 20, 2014 at 7:37 AM, Massimiliano Leoni <span dir="ltr"><<a href="mailto:leoni.massimiliano1@gmail.com" target="_blank">leoni.massimiliano1@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Matt, thanks for your reply.<br>
<br>
In data mercoledì 19 novembre 2014 06:19:05, Matthew Knepley ha scritto:<br>
> It looks to me like you have memory corruption in a different part of the<br>
> code. I<br>
> would build a debugging version, and consider running under valgrind.<br>
I think I got what's happening.<br>
I found out this morning that if I set ksp_type multiplicative in the source<br>
code, with PCFieldSplitSetType, it works!<br>
This suggests that I am messing in the code with the relative positions of<br>
{KSP|PC}SetFromOptions, SetUp and the rest.<br>
<br>
Ksp_view confirms that PETSc actually does what I want, and ksp_monitor shows<br>
different iterations, which means that I am actually influencing something.<br></blockquote><div><br></div><div>Great, I am glad something is working. This is tricky, and it is why in our own code</div><div>we have moved almost completely to an options driven configuration, with SetUp()</div><div>called automatically by the solver. I have not been able to come up with a better option.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
> Here is a sample run<br>
><br>
>   cd src/snes/examples/tutorials<br>
>   make ex62<br>
>   ./ex62 -run_type full -refinement_limit 0.00625 -bc_type dirichlet<br>
> -interpolate 1 -vel_petscspace_order 2 -pres_petscspace_order 1 -ksp_type<br>
> fgmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -pc_type fieldsplit<br>
> -pc_fieldsplit_type schur -pc_fieldsplit_schur_factorization_type full<br>
> -fieldsplit_pressure_ksp_rtol 1e-10 -fieldsplit_velocity_ksp_type gmres<br>
> -fieldsplit_velocity_pc_type lu -fieldsplit_pressure_pc_type jacobi<br>
> -snes_monitor_short -ksp_monitor_short -snes_converged_reason<br>
> -ksp_converged_reason -snes_view -show_solution 0<br>
<br>
I can confirm this works for me as well, of course.<br>
<br>
>   Thanks,<br>
><br>
>      Matt<br>
<br>
Thanks for the help!<br>
<span class="HOEnZb"><font color="#888888">Massimiliano<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>