[petsc-users] Keep symbolic factorization

Matthew Knepley knepley at gmail.com
Thu Sep 25 10:51:24 CDT 2014

On Thu, Sep 25, 2014 at 10:46 AM, Mathieu MORLIGHEM <
Mathieu.Morlighem at uci.edu> wrote:

> Hi,
> I am using an augmented Lagrangian approach that requires to solve
> multiple times a system with the same Stiffness matrix but different right
> hand sides:
> K U1 = F1
> K U2 = F2
>> K UN = FN
> I am using MUMPS right now, and I would like to do the analysis and
> factorization only once (to compute U1), and then only call the back
> substitution in order to speed up the convergence.
> Is there any option in PETSc to do that? If yes, would it work with other
> direct solvers (cholmod, Pardiso)?

In general, this works automatically, and will work with all
factorizations. If
there is a problem with one, let us know.



> Thanks a lot for your help
> --
> Mathieu Morlighem
> Assistant Professor
> Department of Earth System Science
> University of California, Irvine
> 3218 Croul Hall, Irvine, CA 92697-3100
> Mathieu.Morlighem at uci.edu

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/20140925/d0ed5898/attachment.html>

More information about the petsc-users mailing list