<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, May 14, 2018 at 5:45 PM, Salazar De Troya, Miguel <span dir="ltr"><<a href="mailto:salazardetro1@llnl.gov" target="_blank">salazardetro1@llnl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="m_5444543013086775529WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hello,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Up until now, I have been solving my elasticity problem using the gamg preconditioner with the following options:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-ksp_monitor_true_residual <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-ksp_converged_reason<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-ksp_type cg<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-log_view<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-mg_levels_esteig_ksp_type cg<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-mg_levels_ksp_chebyshev_<wbr>esteig_steps 10<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-mg_levels_ksp_type chebyshev<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-mg_levels_pc_type sor<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-pc_type gamg<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-pc_gamg_verbose 7<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-pc_gamg_type agg<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-pc_gamg_agg_nsmooths 1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-pc_gamg_threshold 0.001<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-snes_linesearch_type basic<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="ES" style="font-size:11.0pt">-snes_atol 1e-6 <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="ES" style="font-size:11.0pt">-ksp_atol 1e-7 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-ksp_rtol 1e-9<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-ksp_norm_type unpreconditioned <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">It has worked great so far, even for my problem where I have two materials with more than 1e3 ratio.
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Now, when I add robin boundary conditions, the solver fails with DIVERGED_INDEFINITE_PC. Using -ksp_type bcgs makes the solver converge. I imagine it is because it can handle the indefinite PC. Why is the
 PC indefinite for the Robin problem?</span></p></div></div></blockquote><div><br></div><div>It matters how you are enforcing the conditions. Just write out the problem completely for 1 vertex constrained. It should be clear whether its symmetric or not.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_5444543013086775529WordSection1"><p class="MsoNormal"><span style="font-size:11.0pt"> Is there any way to make it positive definite?</span></p></div></div></blockquote><div><br></div><div>It should be possible if the BC are imposed in the right way.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_5444543013086775529WordSection1"><p class="MsoNormal"><span style="font-size:11.0pt"> I am also passing the rigid body modes of my mesh to the solver. I am not sure if these modes change with the Robin boundary conditions, they should not, should they?</span></p></div></div></blockquote><div><br></div><div>No, these are precisely the null modes of the symbol, and not affected by BC.</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"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_5444543013086775529WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Thanks<span class="HOEnZb"><font color="#888888"><u></u><u></u></font></span></span></p><span class="HOEnZb"><font color="#888888">
<p class="MsoNormal"><span style="font-size:11.0pt">Miguel<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Times New Roman",serif">-- <u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
</font></span></div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>