
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 13 Oct 2019, at 3:48 PM, Mark Adams <<a href="mailto:mfadams@lbl.gov" class="">mfadams@lbl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Oh I see you have established that it is the overlap and got something that looks like the load on a homogenous BC and the processor boundary ... very odd.<div class=""><br class=""></div><div class="">Does the bad version converge to the correct solution if you iterate? It looks like your BC is wrong.</div></div></div></blockquote><div><br class=""></div><div>Yeah, of course: I switched back from P_3 to P_1 and now the first iterate is OK.</div><div>I am indeed most definitely messing up with the RHS at some point before the KSPSolve.</div><div><br class=""></div><div>Thanks!</div><div>Pierre</div><br class=""><blockquote type="cite" class=""><div class=""><div dir="ltr" class=""><div class="">It sure looks like a bug. Can you valgrind this?</div><div class=""><br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 13, 2019 at 9:34 AM Mark Adams <<a href="mailto:mfadams@lbl.gov" class="">mfadams@lbl.gov</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr" class="">Try <span style="font-family: Times; font-size: inherit;" class="">-pc_asm_overlap 0 with ASM.</span><div class=""><span style="font-family: Times; font-size: inherit;" class=""><br class=""></span></div><div class=""><span style="font-family: Times; font-size: inherit;" class="">And I trust the KO run works with 1 processor (direct solve)</span></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 13, 2019 at 3:41 AM Pierre Jolivet via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov" target="_blank" class="">petsc-dev@mcs.anl.gov</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="">Hello,<div class="">I’m struggling to understand the following weirdness with PCASM with exact subdomain solvers.</div><div class="">I’m dealing with a very simple Poisson problem with Dirichlet + Neumann BCs.</div><div class="">If I use PCBJACOBI + KSPPREONLY or 1 iteration of GMRES either preconditioned on the right or on the left, I get the expected result, cf. attached screenshot.</div><div class="">If I use PCASM + KSPPREONLY or 1 iteration of GMRES preconditioned on the left, I get the expected result as well.</div><div class="">However, with PCASM + 1 iteration of GMRES preconditioned on the right, I don’t get what I should (I believe).</div><div class="">Furthermore, this problem is specific to -pc_asm_type restrict,none (I get the expected result with basic,interpolate).</div><div class=""><br class=""></div><div class="">Any hint?</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Pierre</div><div class=""><br class=""></div><div class="">$ -sub_pc_type lu -ksp_max_it 1 -ksp_type gmres -pc_type bjacobi -ksp_pc_side right -> bjacobi_OK</div><div class="">$ -sub_pc_type lu -ksp_max_it 1 -ksp_type gmres -pc_type asm -ksp_pc_side left -> asm_OK</div><div class="">$ -sub_pc_type lu -ksp_max_it 1 -ksp_type gmres -pc_type asm -ksp_pc_side right -> asm_KO</div><div class=""><br class=""></div><div class=""><span id="cid:16dc550492b2dc4f06a1"><bjacobi_OK.png></span><span id="cid:16dc550492b91e8295b2"><asm_OK.png></span><span id="cid:16dc550492b8b4d24e33"><asm_KO.png></span></div></div></blockquote></div>

</blockquote></div>

</div></blockquote></div><br class=""></body></html>