CAMLprim value ml_gsl_linalg_QR_QRsolve(value Q, value R, value B, value X) { _DECLARE_MATRIX2(Q, R); _DECLARE_VECTOR2(B, X); _CONVERT_MATRIX2(Q, R); _CONVERT_VECTOR2(B, X); gsl_linalg_QR_QRsolve(&m_Q, &m_R, &v_B, &v_X); return Val_unit; }
/** * C++ version of gsl_linalg_QR_QRsolve(). * @param Q A Matrix * @param R A Matrix * @param b A vector * @param x A vector * @return Error code on failure */ inline int QR_QRsolve( matrix& Q, matrix& R, vector const& b, vector& x ){ return gsl_linalg_QR_QRsolve( Q.get(), R.get(), b.get(), x.get() ); }