CAMLprim value ml_gsl_linalg_LQ_LQsolve(value Q, value L, value B, value X) { _DECLARE_MATRIX2(Q,L); _DECLARE_VECTOR2(B,X); _CONVERT_MATRIX2(Q,L); _CONVERT_VECTOR2(B,X); gsl_linalg_LQ_LQsolve(&m_Q, &m_L, &v_B, &v_X); return Val_unit; }
/** * C++ version of gsl_linalg_LQ_LQsolve(). * @param Q A matrix * @param L A matrix * @param b A vector * @param x A vector * @return Error code on failure */ inline int LQ_LQsolve( matrix& Q, matrix& L, vector const& b, vector& x ){ return gsl_linalg_LQ_LQsolve( Q.get(), L.get(), b.get(), x.get() ); }