Invert a matrix using the Gauss-Jordan algorithm. The 1- matrix is not set up and use is made of its structure.
I n v e r t
a
m a t r i x
\ ( A \ )
u s i n g
t h e
G a u s s - J o r d a n
a l g o r i t h m .
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(out), | dimension(arraydim, arraydim) | :: | ainv |
Inverted matrix |
|
real(kind=dp), | intent(inout), | dimension(arraydim, arraydim) | :: | a |
Matrix to be inverted |
|
integer, | intent(in) | :: | arraydim |
Dimension of the array |
||
integer, | intent(in) | :: | n |
Number of columns |