[petsc-dev] Computing inverse of a small matrix: MatLUFactor() usage

[Jed Brown]

On Wed, Nov 21, 2012 at 9:03 PM, Gautam Bisht <gbisht at lbl.gov> wrote:

> I need to compute an inverse of M, a small matrix (3x3), and am coding in
> Fortran 90. Following the suggestion under FAQ on the website, I'm trying
> to use MatMatSolve to obtain the inverse of matrix. I did the following:
>
> - Created three sequential dense matrices: M, I, Minv
> - Filled values in M and I using MatSetValues
> - Tried factoring Matrix M before calling MatMatSolve(M,I,Minv,ierr)
>
> I'm not sure what arguments needs to be passed as arguments in
> MatLUFactor() corresponding to IS row, IS col, and MatFactorInfo. I noticed
> the following note on website that states "*Developer Note: fortran
> interface is not autogenerated as the f90 interface defintion cannot be
> generated correctly [due to MatFactorInfo]*".
> Is there a Fortran90 example demonstrating the use of MatMatSolve()?
>

For a 3x3 matrix, you should just use the explicit formula for the inverse.
Going through the PETSc interface doesn't make sense at all.


>
> Additionally, I have been trying to update my copy of petsc-dev, and am
> getting timed out from the server.
>

The server died so it has moved to https://bitbucket.org/petsc/petsc-dev.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121121/8531ac95/attachment.html>