<div dir="ltr">Anton,<div>I pushed a fix </div><div><a href="https://bitbucket.org/petsc/petsc/commits/28865de08051eb99557d70672c208e14da23c8b1">https://bitbucket.org/petsc/petsc/commits/28865de08051eb99557d70672c208e14da23c8b1</a></div><div>in branch hzhang/fix-superlu_dist-reuse-factornumeric.</div><div>Can you give it a try to see if it works?</div><div>I do not have an example which produces your problem.</div><div><br></div><div><span style="font-size:12.8px">In your email, you asked "Setting Options.Fact = DOFACT for all factorizations is currently impossible via PETSc interface.</span><br style="font-size:12.8px"><span style="font-size:12.8px">The user is expected to choose some kind of reuse model.</span><br style="font-size:12.8px"><span style="font-size:12.8px">If you could add it, I (and other users probably too) would really appreciate that."</span><br></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">We do not allow user set superlu' </span><span style="font-size:12.8px">Options.Fact = DOFACT. If user changes matrix structure, then user must call</span></div><div><span style="font-size:12.8px">KSPSetOperators() -> call symbolic matrix factorization again, in which we set </span><span style="font-size:12.8px">Options.Fact = DOFACT.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Hong</span></div></div>