rinvgj Subroutine

public subroutine rinvgj(ainv, a, arraydim, n)

Invert a matrix using the Gauss-Jordan algorithm. The 1- matrix is not set up and use is made of its structure.

Arguments

Type IntentOptional 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