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.