Function rgsl::linear_algebra::QRPT_Rsvx

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

This function solves the triangular system R P^T x = b in-place for the N-by-N matrix R contained in QR. On input x should contain the right-hand side b, which is replaced by the solution on output.