[petsc-users] Different factorizations within one scope

[Matthew Knepley]

On Tue, May 31, 2011 at 4:16 AM, Alexander Grayver
<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/5a943e5c/attachment.htm>