diff --git a/tests2/fn_spsolve.cpp b/tests2/fn_spsolve.cpp
index 077a3ced..16c981b3 100644
--- a/tests2/fn_spsolve.cpp
+++ b/tests2/fn_spsolve.cpp
@@ -54,7 +54,7 @@ TEST_CASE("fn_spsolve_sparse_test")
       {
       for (uword j = 0; j < dX.n_rows; ++j)
         {
-        REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)) );
+        REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)).epsilon(0.01) );
         }
       }
     }
@@ -94,7 +94,7 @@ TEST_CASE("fn_spsolve_sparse_nonsymmetric_test")
       {
       for (uword j = 0; j < dX.n_rows; ++j)
         {
-        REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)) );
+        REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)).epsilon(0.01) );
         }
       }
     }
@@ -135,7 +135,7 @@ TEST_CASE("fn_spsolve_sparse_float_test")
       {
       for (size_t j = 0; j < dX.n_rows; ++j)
         {
-        REQUIRE( (float) X(j, i) == Approx((float) dX(j, i)) );
+        REQUIRE( (float) X(j, i) == Approx((float) dX(j, i)).epsilon(0.01) );
         }
       }
     }
@@ -175,7 +175,7 @@ TEST_CASE("fn_spsolve_sparse_nonsymmetric_float_test")
       {
       for (uword j = 0; j < dX.n_rows; ++j)
         {
-        REQUIRE( (float) X(j, i) == Approx((float) dX(j, i)) );
+        REQUIRE( (float) X(j, i) == Approx((float) dX(j, i)).epsilon(0.01) );
         }
       }
     }
@@ -217,7 +217,7 @@ TEST_CASE("fn_spsolve_sparse_complex_float_test")
       for (uword j = 0; j < dX.n_rows; ++j)
         {
         REQUIRE( (float) std::abs((cx_float) X(j, i)) ==
-                 Approx((float) std::abs((cx_float) dX(j, i))) );
+                 Approx((float) std::abs((cx_float) dX(j, i))).epsilon(0.01) );
         }
       }
     }
@@ -258,7 +258,7 @@ TEST_CASE("fn_spsolve_sparse_nonsymmetric_complex_float_test")
       for (uword j = 0; j < dX.n_rows; ++j)
         {
         REQUIRE( (float) std::abs((cx_float) X(j, i)) ==
-                 Approx((float) std::abs((cx_float) dX(j, i))) );
+                 Approx((float) std::abs((cx_float) dX(j, i))).epsilon(0.01) );
         }
       }
     }
@@ -300,7 +300,7 @@ TEST_CASE("fn_spsolve_sparse_complex_test")
       for (uword j = 0; j < dX.n_rows; ++j)
         {
         REQUIRE( (double) std::abs((cx_double) X(j, i)) ==
-                 Approx((double) std::abs((cx_double) dX(j, i))) );
+                 Approx((double) std::abs((cx_double) dX(j, i))).epsilon(0.01) );
         }
       }
     }
@@ -341,7 +341,7 @@ TEST_CASE("fn_spsolve_sparse_nonsymmetric_complex_test")
       for (uword j = 0; j < dX.n_rows; ++j)
         {
         REQUIRE( (double) std::abs((cx_double) X(j, i)) ==
-                 Approx((double) std::abs((cx_double) dX(j, i))) );
+                 Approx((double) std::abs((cx_double) dX(j, i))).epsilon(0.01) );
         }
       }
     }
@@ -376,7 +376,7 @@ TEST_CASE("fn_spsolve_delayed_sparse_test")
     {
     for (uword j = 0; j < dX.n_rows; ++j)
       {
-      REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)) );
+      REQUIRE( (double) X(j, i) == Approx((double) dX(j, i)).epsilon(0.01) );
       }
     }
   }
@@ -420,7 +420,7 @@ TEST_CASE("fn_spsolve_superlu_solve_test")
     {
     for (uword j = 0; j < x.n_rows; ++j)
       {
-      REQUIRE( (double) x(j, i) == Approx(dx(j, i)) );
+      REQUIRE( (double) x(j, i) == Approx(dx(j, i)).epsilon(0.01) );
       }
     }
   }
@@ -451,7 +451,7 @@ TEST_CASE("fn_spsolve_random_superlu_solve_test")
       {
       for (uword j = 0; j < x.n_rows; ++j)
         {
-        REQUIRE( x(j, i) == Approx((double) trueX(j, i)) );
+        REQUIRE( x(j, i) == Approx((double) trueX(j, i)).epsilon(0.01) );
         }
       }
     }
@@ -496,7 +496,7 @@ TEST_CASE("fn_spsolve_float_superlu_solve_test")
     {
     for (uword j = 0; j < x.n_rows; ++j)
       {
-      REQUIRE( (float) x(j, i) == Approx(dx(j, i)) );
+      REQUIRE( (float) x(j, i) == Approx(dx(j, i)).epsilon(0.01) );
       }
     }
   }
@@ -754,7 +754,7 @@ TEST_CASE("fn_spsolve_function_test")
     {
     for (uword j = 0; j < x.n_rows; ++j)
       {
-      REQUIRE( (double) trueX(j, i) == Approx(x(j, i)) );
+      REQUIRE( (double) trueX(j, i) == Approx(x(j, i)).epsilon(0.01) );
       }
     }
   }