pub(super) fn linsolve( n: usize, a: &mut [f64], stride: usize, b: &mut [f64], x: &mut [f64], ) -> bool
Solves Ax = b, where x and b are column vectors of size nx1 and A is nxn