Test of matrix inversion on a GPU, main programm is writen in FORTRAN, GPU kernel and wrapper is in C.
We use partial pivoting, which selects an element with the largest absolute value from a column as a pivot. To understand the basic idea behind the code, subroutines related to this are not important (switchRows, createIndexVector, storeRows, createPivotVector) and can be commented out.