Thank you Matthew, that's what I wanted to hear. Regards, Alexander
On 31.05.2011 15:13, Matthew Knepley wrote: > On Tue, May 31, 2011 at 4:16 AM, Alexander Grayver > <agrayver at gfz-potsdam.de <mailto:agrayver at gfz-potsdam.de>> wrote: > > Hello, > > 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. > 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. > > My question is that possible to create several factorizations > using PETSc classes and MUMPS package, store them and use later > any you need? > > > 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, > and it will preserve the existing factorization. > > Matt > > Regards, > Alexander > > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which > their experiments lead. > -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110531/52a69c6e/attachment.htm>
