[petsc-users] Different factorizations within one scope

[Alexander Grayver]

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>