decomplu.go

func MatrixLUDecomposition(matrix Matrix) Matrix {
	columns := matrix.cols()
	rows := matrix.rows()

	for i := 0; i < columns-1; i++ {
		val, pivot_index := PivotIndex(matrix, i, i)
		MatrixExchangeRows(matrix, i, pivot_index)

		for j := i + 1; j < rows; j++ {
			multiplier := -matrix.at(j, i) / val
			MatrixScaleAndSumRows(matrix, i, j, multiplier)
			matrix.set(j, i, -(multiplier))
		}

	}

	return matrix
}