pub fn solve_pivoted<const C: usize>(
a: Vec<Vec<Scalar>>,
rhs: Vec<Vector<C>>,
) -> GeopResult<Vec<Vector<C>>>Expand description
Solves a x = rhs for the n x n matrix a and C right-hand sides
by Gaussian elimination, pivoting on the entry of largest magnitude in
each column. Which row to pivot on is a free choice — every pivot that
cannot be zero gives the same enclosed solution — and the largest is the
well-conditioned one. Fails if no entry of a column definitely differs
from zero.