[petsc-users] Inverse matrix of MATMPIDENSE matrix

[Joon Hee Choi]

Thank you for your reply.
However, I followed the manual using MatLUFactor() and MatMatSolve(). However, MatLUFactor() did not support MATMPIDENSE. If my code is wrong, then please let me know. My code is following:

MatGetOrdering(invM2T, MATORDERINGNATURAL, &isr, &isc);
MatFactorInfoInitialize(&info);
MatLUFactor(invM2T, isr, isc, &info);
MatMatSolve(invM2T, M1T, *M3T);

If not so, then please let me know which function I have to use. Thank you.

Sincerely,
Joon

----- Original Message -----
From: "Barry Smith" <bsmith at mcs.anl.gov>
To: "Joon Hee Choi" <choi240 at purdue.edu>
Cc: petsc-users at mcs.anl.gov
Sent: Tuesday, September 17, 2013 7:59:42 AM
Subject: Re: [petsc-users] Inverse matrix of MATMPIDENSE matrix

http://www.mcs.anl.gov/petsc/documentation/faq.html#invertmatrix


On Sep 17, 2013, at 1:07 AM, Joon Hee Choi <choi240 at purdue.edu> wrote:

> Hello,
> 
> I want to get the inverse matrix of 10 x 10 MATMPIDENSE matrix. I know I cannot use MatLUFactor() for MATMPIDENSE. Is there any other method? Could you let me know the fastest way to get the inverse? Also, I don't want to rebuild or configure something for the inverse because I use petsc in my school server and I don't have permission to configure.
> 
> Thanks,
> Joon