<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#ffffff" text="#000000">
    Thank you Matthew, that's what I wanted to hear.<br>
    <br>
    Regards,<br>
    Alexander<br>
    <br>
    On 31.05.2011 15:13, Matthew Knepley wrote:
    <blockquote
      cite="mid:BANLkTimMd-g2O5cemy3W8ZF2zKpjgpVVsg@mail.gmail.com"
      type="cite">On Tue, May 31, 2011 at 4:16 AM, Alexander Grayver <span
        dir="ltr">&lt;<a moz-do-not-send="true"
          href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>&gt;</span>
      wrote:<br>
      <div class="gmail_quote">
        <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
          0.8ex; border-left: 1px solid rgb(204, 204, 204);
          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>&nbsp; &nbsp;Matt</div>
        <div>&nbsp;</div>
        <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
          0.8ex; border-left: 1px solid rgb(204, 204, 204);
          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>
    </blockquote>
    <br>
  </body>
</html>