On Tue, May 31, 2011 at 4:16 AM, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hello,<br>
<br>
Now I use MUMPS to solve Helmholtz equation for electromagnetic problem in frequency domain. My system matrix depends on frequency, since I have more than one frequency I need to do factorization for every one.<br>
Within one MPI communicator it is implemented as a loop over frequencies, where I assemble new system matrix at each iteration, factorize it and solve system for many RHS. The problem is that at each iteration previous factorization has been overwritten by new one, but now I came to point where I need factorizations for previous frequencies later in my program and that means that they need to be stored somewhere.<br>
<br>
My question is that possible to create several factorizations using PETSc classes and MUMPS package, store them and use later any you need?<br></blockquote><div><br></div><div>The easiest way to do this is to create several KSP objects, which will have associated PC and Mat objects. You can resolve with each one later,</div>
<div>and it will preserve the existing factorization.</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;">
Regards,<br><font color="#888888">
Alexander<br>
</font></blockquote></div><br><br clear="all"><br>-- <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<br>