@MSGID:
<4d061409-7c22-4e8c-a020-6e4476e8f165n@googlegroups.com> 90ca26cc
@REPLY:
<f8466a68-55d1-4d8c-a607-5c41c85404d2n@googlegroups.com> c559ed72
@REPLYADDR gah4 <gah4@u.washington.edu>
@REPLYTO 2:5075/128 gah4
@CHRS: CP866 2
@RFC: 1 0
@RFC-References:
<f8466a68-55d1-4d8c-a607-5c41c85404d2n@googlegroups.com>
@RFC-Message-ID:
<4d061409-7c22-4e8c-a020-6e4476e8f165n@googlegroups.com>
@TZUTC: -0700
@PID: G2/1.0
@TID: FIDOGATE-5.12-ge4e8b94
On Thursday, August 31, 2023 at 7:31:38 AM UTC-7, Farzad Tatar wrote:
 
> I am trying to use LAPack libraries in Fortran to invert a square matrix. 
(snip)

> real(DP):: TEST(2,2) = TRANSPOSE(reshape((/1,2,3,4/), (/2,2/))) 
> real(DP):: TEST2(2) = (/1,2/) 
> PRINT*, `TEST , TEST2 ARE:`, TEST, TEST2 
> CALL DGESV(2, 1, TEST, 2, IPIV, TEST2, 2, INFO) 
> PRINT*, `AFTER THE INVERSION, TEST and TEST2 ARE:`, TEST, TEST2 
 
Seems to me that there is an easy way to get the inverse matrix.

Set TEST2 to an identity matrix, and so the solution will be the inverse.
--- G2/1.0
 * Origin: usenet.network (2:5075/128)
SEEN-BY: 5001/100 5005/49 5015/255 5019/40 5020/715
848 1042 4441 12000
SEEN-BY: 5030/49 1081 5058/104 5075/128
@PATH: 5075/128 5020/1042 4441