mod_rinvgj Module

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


Uses

    • mod_datatypes

Subroutines

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

License
Creative Commons License
Category
numerical-tools, KKRhost

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 .

Read more…

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