diff --git a/src/svm/svc.rs b/src/svm/svc.rs
index 04b3b7b..119b812 100644
--- a/src/svm/svc.rs
+++ b/src/svm/svc.rs
@@ -247,7 +247,7 @@ impl<T: RealNumber, M: Matrix<T>, K: Kernel<T, M::RowVector>> SVC<T, M, K> {
 
 impl<T: RealNumber, M: Matrix<T>, K: Kernel<T, M::RowVector>> PartialEq for SVC<T, M, K> {
     fn eq(&self, other: &Self) -> bool {
-        if self.b != other.b
+        if (self.b - other.b).abs() > T::epsilon() * T::two()
             || self.w.len() != other.w.len()
             || self.instances.len() != other.instances.len()
         {
diff --git a/src/svm/svr.rs b/src/svm/svr.rs
index 0fcaa30..61feb80 100644
--- a/src/svm/svr.rs
+++ b/src/svm/svr.rs
@@ -204,7 +204,7 @@ impl<T: RealNumber, M: Matrix<T>, K: Kernel<T, M::RowVector>> SVR<T, M, K> {
 
 impl<T: RealNumber, M: Matrix<T>, K: Kernel<T, M::RowVector>> PartialEq for SVR<T, M, K> {
     fn eq(&self, other: &Self) -> bool {
-        if self.b != other.b
+        if (self.b - other.b).abs() > T::epsilon() * T::two()
             || self.w.len() != other.w.len()
             || self.instances.len() != other.instances.len()
         {