pub fn matrix_approx_eq_float<T: PartialEq + NumAssign + Copy + Default + Float + Debug, const ROW: usize, const COL: usize>( a: &Matrix<T, ROW, COL>, b: &Matrix<T, ROW, COL>, tol: T, ) -> bool