CAMLprim value ml_gsl_linalg_PTLQ_solve_T (value QR, value TAU, value P, value B, value X)
{
_DECLARE_MATRIX(QR);
_DECLARE_VECTOR3(TAU,B,X);
GSL_PERMUT_OF_BIGARRAY(P);
_CONVERT_MATRIX(QR);
_CONVERT_VECTOR3(TAU,B,X);
gsl_linalg_PTLQ_solve_T (&m_QR, &v_TAU, &perm_P, &v_B, &v_X);
return Val_unit;
}
/**
* C++ version of gsl_linalg_PTLQ_solve_T().
* @param QR A QR decomposition matrix
* @param tau A vector
* @param p A permutation
* @param b A vector
* @param x A vector
* @return Error code on failure
*/
inline int PTLQ_solve_T( matrix const& QR, vector const& tau, permutation const& p,
vector const& b, vector& x ){
return gsl_linalg_PTLQ_solve_T( QR.get(), tau.get(), p.get(), b.get(), x.get() ); }