Function rgsl::linear_algebra::QRPT_solve

pub fn QRPT_solve(
    qr: &MatrixF64, 
    tau: &VectorF64, 
    p: &Permutation, 
    b: &VectorF64, 
    x: &mut VectorF64
) -> Value

This function solves the square system A x = b using the QRP^T decomposition of A held in (QR, tau, p) which must have been computed previously by QRPT_decomp.