Function LU_solve (o2scl_linalg)

O2scl : Function List

template<class mat_t, class vec_t, class vec2_t>
int o2scl_linalg::LU_solve(const size_t N, const mat_t &LU, const o2scl::permutation &p, const vec_t &b, vec2_t &x)

Solve a linear system after LU decomposition.

This function solve the square system A x = b using the LU decomposition of A into (LU, p) given by gsl_linalg_LU_decomp or gsl_linalg_complex_LU_decomp.