lu_decompose Function

public function lu_decompose(a) result(m)

Carry out LU-decomposition for A and put the decomposed matrix into M implementation of Crout's method with partial pivoting.

Arguments

Type IntentOptional Attributes Name
real(kind=long), intent(in), dimension(:, :) :: a

n*N matrix

Return Value type(lu_matrix)