Function rgsl::linear_algebra::QRPT_QRsolve

pub fn QRPT_QRsolve(
    q: &MatrixF64, 
    r: &MatrixF64, 
    p: &Permutation, 
    b: &VectorF64, 
    x: &mut VectorF64
) -> Value

This function solves the square system R P^T x = Q^T b for x. It can be used when the QR decomposition of a matrix is available in unpacked form as (Q, R).