Function sparse_upper_triangular_solve 

pub fn sparse_upper_triangular_solve(u: &CsrMatrix, b: &[f64]) -> Vec<f64>

Back-substitution solve U x = b where U is upper-triangular (CSR).

Assumes U has non-zero diagonal entries.