Function rgsl::linear_algebra::QRPT_svx

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

This function solves the square system A x = b in-place using the QRP^T decomposition of A held in (QR,tau,p). On input x should contain the right-hand side b, which is replaced by the solution on output.