Scythe-1.0.3
template<matrix_order RO, matrix_style RS, typename T , matrix_order PO1, matrix_style PS1, matrix_order PO2, matrix_style PS2>
Matrix<T,RO,RS> scythe::lu_solve ( Matrix< T, PO1, PS1 >  A,
const Matrix< T, PO2, PS2 > &  b 
)

Solve $Ax=b$ for x via forward and backward substitution, using LU decomposition.

This function solves the system of equations $Ax = b$ via forward and backward substitution and LU decomposition. A must be a non-singular square matrix for this method to work.

Parameters:
AA non-singular square Matrix to decompose.
bA column vector with as many rows as A.
See also:
lu_solve (const Matrix<T,PO1,PS1>&, const Matrix<T,PO2,PS2>&, const Matrix<T,PO3,PS3>&, const Matrix<T,PO4,PS4>&, const Matrix<unsigned int, PO5, PS5>&)
lu_decomp(Matrix<T,PO1,PS1>, Matrix<T,PO2,Concrete>&, Matrix<T,PO3,Concrete>&, Matrix<unsigned int, PO4, Concrete>&)
chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &)
chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &, const Matrix<T,PO3,PS3> &)
Exceptions:
scythe_null_error(Level 1)
scythe_dimension_error(Level 1)
scythe_conformation_error(Level 1)
scythe_type_error(Level 2)