<div class="gmail_extra">On Wed, Nov 21, 2012 at 9:03 PM, Gautam Bisht <span dir="ltr"><<a href="mailto:gbisht@lbl.gov" target="_blank">gbisht@lbl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>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:</div>
<div><br></div><div>- Created three sequential dense matrices: M, I, Minv</div><div>- Filled values in M and I using MatSetValues</div><div>- Tried factoring Matrix M before calling MatMatSolve(M,I,Minv,ierr)</div><div><br>
</div><div>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 "<b>Developer Note: fortran interface is not autogenerated as the f90 interface defintion cannot be generated correctly [due to MatFactorInfo]</b>".</div>
<div>Is there a Fortran90 example demonstrating the use of MatMatSolve()?</div></blockquote><div><br></div><div>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.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br></div><div>Additionally, I have been trying to update my copy of petsc-dev, and am getting timed out from the server. </div>
</blockquote></div><br></div><div class="gmail_extra">The server died so it has moved to <a href="https://bitbucket.org/petsc/petsc-dev">https://bitbucket.org/petsc/petsc-dev</a>.</div>