<div dir="ltr"><div dir="ltr">On Thu, May 30, 2019 at 6:08 AM Pierre Jolivet via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello,<br>
I’m doing a continuation loop with something like:<br>
KSP ksp;<br>
KSPSetUp(ksp);<br>
while(cond) {<br>
        SNES snes;<br>
        SNESCreate(&snes);<br>
        SNESSetKSP(snes, &ksp);<br>
        SNESSolve(snes);<br>
        SNESDestroy(&snes);<br>
}<br>
<br>
For the first iteration, everything is OK.<br>
After that, it looks like SNESDestroy is messing up with my outer-defined KSP(*), i.e., for the code to work I need to comment out the SNESDestroy.<br>
Is this the expected behavior? (if not, I’ll provide a MWE)<br>
How to fix this properly, without having leaks by not destroying the inner-defined SNESes?<br>
<br>
Thanks in advance,<br>
Pierre<br>
<br>
(*) to the outer-defined KSP is attached a PCFIELDSPLIT, first SNESSolve, everything OK:<br>
PC Object: 1 MPI processes<br>
  type: fieldsplit<br>
    FieldSplit with MULTIPLICATIVE composition: total splits = 2<br>
<br>
Second SNESSolve:<br>
[0]PETSC ERROR: PCFieldSplitSetDefaults() Unhandled case, must have at least two fields, not 1</blockquote></div><br clear="all"><div>I believe SNESDestroy() calls SNESReset(), which calls KSPReset(), which is wiping out your information.</div><div><br></div><div>This pattern is in a lot of places. Destroy calls Reset because it is natural to reuse the code which tears down</div><div>the internal structures. Reset calls reset on subobjects because we also want to tear down their internal structures (usually).</div><div>Reset wouldn't do what we want on resizing if we left subobjects unreset. However, this chain is causing subobjects to</div><div>be wiped clean on Destroy, so that we cannot reuse subobjects.</div><div><br></div><div>It would not even work if you change the Destroy to Reset in the continuation loop. However, I think you can hoist all</div><div>the SNES setup outside the loop. If you are changing sizes, then the KSP needs to be wiped out anyway. If not, then</div><div>you do not need to wipe out the SNES.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>