Tolerance check of absolute value in 'matrix_eq'

1 files changed, 1 insertions, 1 deletions

diff --git a/matrix.c b/matrix.c
index 46302f7..1e1a4dd 100644
--- a/matrix.c
+++ b/matrix.c
@@ -247,7 +247,7 @@ int matrix_eq(const Matrix *A, const Matrix *B, const double tol)
     for (size_t i = 0; i < A->m; ++i)
     {
         for (size_t j = 0; j < A->n; ++j)
-            if (matrix_at(A, i, j) - matrix_at(B, i, j) > tol) return FALSE;   
+            if (ABS(matrix_at(A, i, j) - matrix_at(B, i, j)) > tol) return FALSE;
     }
     return TRUE;
 }