Invert a matrix using the Gauss-Jordan algorithm. The 1- matrix is not set up and use is made of its structure.
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 |