<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 25, 2014 at 10:46 AM, Mathieu MORLIGHEM <span dir="ltr"><<a href="mailto:Mathieu.Morlighem@uci.edu" target="_blank">Mathieu.Morlighem@uci.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi,<div><br></div><div>I am using an augmented Lagrangian approach that requires to solve multiple times a system with the same Stiffness matrix but different right hand sides:</div><div><br></div><div>K U1 = F1</div><div>K U2 = F2</div><div>…</div><div>K UN = FN</div><div><br></div><div>I am using MUMPS right now, and I would like to do the analysis and factorization only once (to compute U1), and then only call the back substitution in order to speed up the convergence.</div><div>Is there any option in PETSc to do that? If yes, would it work with other direct solvers (cholmod, Pardiso)?</div></div></blockquote><div><br></div><div>In general, this works automatically, and will work with all factorizations. If</div><div>there is a problem with one, let us know.</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 style="word-wrap:break-word"><div>Thanks a lot for your help</div><div><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">--<br>Mathieu Morlighem<br><br>Assistant Professor<br>Department of Earth System Science<br>University of California, Irvine<br>3218 Croul Hall, Irvine, CA 92697-3100<br><a href="mailto:Mathieu.Morlighem@uci.edu" target="_blank">Mathieu.Morlighem@uci.edu</a></div></div></div></div>
</div>
<br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>