@MSGID: <klcia5FihqfU3@mid.individual.net>
92d534d1
@REPLY:
<f8466a68-55d1-4d8c-a607-5c41c85404d2n@googlegroups.com> c559ed72
@REPLYADDR Giorgio Pastore <pastgio@units.it>
@REPLYTO 2:5075/128 Giorgio Pastore
@CHRS: CP866 2
@RFC: 1 0
@RFC-Message-ID:
<klcia5FihqfU3@mid.individual.net>
@RFC-References:
<f8466a68-55d1-4d8c-a607-5c41c85404d2n@googlegroups.com>
@TZUTC: 0200
@PID: Mozilla/5.0 (Macintosh; Intel Mac OS X
10.15; rv:102.0) Gecko/20100101 Thunderbird/102.14.0
@TID: FIDOGATE-5.12-ge4e8b94
Il 31/08/23 16:31, Farzad Tatar ha scritto:
>   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

You should be aware that dgesv doesn`t invert a matrix. It solves a 
linear system where test is the matrix of coefficients, and test2 the 
array of the known terms. In aoutput test2 is the array containing the 
solution.

Giorgio
 --- Mozilla/5.0 (Macintosh; Intel Mac OS X 10.15; rv:102.0)
Gecko/20100101 Thunderbird/102.14.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