pub fn row_reduce3<S>(m: Matrix3<S>) -> Matrix3<S>where S: OrderedField + RelativeEq<Epsilon = S> + Debug,
Reduce a 3x3 matrix to row echelon form using Gaussian elimination