[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: matrix inversion
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: matrix inversion
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Fri, 21 Sep 2007 10:13:01 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=beta; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=zO8fLKAfy56NDTfAI/kDhbQie2ddwENUjnQ++Ii1spA=; b=k9Y8fi5HBYh5r8ktNbLR1Bz+E8Mlp8nxHOAMi1Vr//6XPVpvLM4/jbHcZMglRGXClt9g41KgIPEfX1p0cq/DjYNHiezs6OTRLJ1SSGqRfPyC6VXS/66H1pHgG7zULy4ukox33psRR6Up0kAfWTziSxYPKVhAzNSC2xMwsBn/rwA=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=SiVA5/XT2Cnt1y+ieLPV+6AUFbtJtaVok9Fb8ofUyltg9HIcJ68VC5MS85gtoW9n4rpiDvvmzpYfCZrwc/hYbJoxnJAEs3Bb4a+fyQgXJV2fDraYtWTD68FcH0o3weV+SF2VUDmDoKbbeu4U/HfkJtRUnHd7zSN3ZIaMwgupdB4=
- In-reply-to: <46F3DEDF.9080307@geosystem.us>
- References: <46F3DC02.20503@geosystem.us> <a9f269830709210806j3b7d9f52tfe8caaab28cd0f99@mail.gmail.com> <46F3DEDF.9080307@geosystem.us>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
On 9/21/07, Randall Mackie <randy@xxxxxxxxxxxx> wrote:
> Matt,
>
> There may be very valid reasons why one would want a matrix inverse,
> say in the context of solving optimizations where you need the
> inverse of the covariance matrix, and you need those values to
> embed in another system.
I do not agree. What are you going to do with the inverse that cannot be
done with LU? Backsolves against a matrix are commonplace.
Matt
> In this case, just getting the action of the inverse times a vector
> is not what is needed.
>
> Randy
>
>
> Matthew Knepley wrote:
> > On 9/21/07, Randall Mackie <randy@xxxxxxxxxxxx> wrote:
> >> I have a collegue who is considering to use PETSc for programming
> >> his problem. He would like to know if you can get the actual
> >> inverse of a matrix, because this inverse is embedded in another
> >> larger system that is then solved.
> >>
> >> I suppose that if you solve the system directly, say with SuperLU
> >> or MUMPS, you could get the actual inverse? I looked briefly through
> >> the manual and documentation, but it didn't jump out at me.
> >
> > This makes several time this week people have asked for inverses. Never
> > never never use the inverse. Its not stable. The storage blows up. That is
> > why factorizations were invented. Numerical analysis courses have failed
> > the country entirely. Just use a KSP and tell it -ksp_type preonly -pc_type lu
> > and then you can customize the matrix type for different solvers. KSPSolve()
> > will apply the matrix. If you want to get the action on an entire matrix, this
> > is dense and we would have to put in the BLAS3 triangular solve calls.
> >
> > Matt
> >
> >> Any advice here?
> >>
> >> Thanks, Randy
> >>
> >> --
> >> Randall Mackie
> >> GSY-USA, Inc.
> >> PMB# 643
> >> 2261 Market St.,
> >> San Francisco, CA 94114-1600
> >> Tel (415) 469-8649
> >> Fax (415) 469-5044
> >>
> >> California Registered Geophysicist
> >> License No. GP 1034
> >>
> >>
> >
> >
>
> --
> Randall Mackie
> GSY-USA, Inc.
> PMB# 643
> 2261 Market St.,
> San Francisco, CA 94114-1600
> Tel (415) 469-8649
> Fax (415) 469-5044
>
> California Registered Geophysicist
> License No. GP 1034
>
>
--
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