<div>I am doing a transient computation, solving one linear problem per timestep, so naturally I want to exploit &#39;x&#39; from the previous time step to be the initial value for the next solve (<span class="Apple-style-span" style="font-family: monospace; white-space: pre; ">KSPSetInitialGuessNonzero)</span>.</div>
<div>For the longest time, however, I was getting wrong results, unless I was resetting &#39;x&#39; each time step (to some constant value, pure zero caused bcgs to break down).</div><div>After hours of debugging I was unable to find any errors in my coefficients, I experimentally found out, however, that changing the solver from bcgs to gmres or fgmres removes the problem: I no longer need to clear the solution vector.</div>
<div>Now I am a bit worried, if this is still some time bomb in my code or is a known phenomenon. Thanks for any hints.</div><div><br></div><div>Regards, Dominik</div>