pub fn utsolve<T: Numeric<T>>(u: &Sprs<T>, x: &mut [T])Expand description
Solves U’x=b where x and b are dense.
x=b on input, solution on output.
§Example:
let u = [
vec![0.9842, 0.1720, 0.9948, 0.2766, 0.4560, 0.1462, 0.8124],
vec![0.0000, 0.6894, 0.1043, 0.4486, 0.5217, 0.7157, 0.4132],
vec![0.0000, 0.0000, -0.5500, -0.2340, 0.0822, 0.2176, -0.1996],
vec![0.0000, 0.0000, 0.0000, 0.6554, -0.1564, -0.0287, 0.2107],
vec![0.0000, 0.0000, 0.0000, 0.0000, -0.4127, -0.4652, -0.6993],
vec![0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.6881, 0.3037],
vec![0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, -0.7740]
];
let mut u_sparse = rsparse::data::Sprs::new();
u_sparse.from_vec(&u);
let mut b = [
0.444841,
0.528773,
0.988345,
0.097749,
0.996166,
0.068040,
0.844511
];
rsparse::utsolve(&u_sparse, &mut b);