lu_solve Function

public function lu_solve(m, b) result(x)

Solve a linear equation for b with LU-decomposed matrix M Return the solution for b as a new vector x b : row vector, not column vector, that is, just 1D array of n

Arguments

Type	Intent	Optional	Attributes		Name
type(lu_matrix),	intent(in)			::	m
real(kind=long),	intent(in),		dimension(:)	::	b

Return Value real(kind=long), dimension(m%n)