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

[Gautam Bisht]

Hi All,

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()?

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

Thanks.
-Gautam.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121121/a969a2b8/attachment.html>