pub fn compare_matrices( a: &DMatrix<f64>, b: &DMatrix<f64>, tolerance: f64, ) -> ComparisonResult
Compare two matrices element-wise with configurable tolerance.