diff --git a/src/psimpy/inference/bayes_inference.py b/src/psimpy/inference/bayes_inference.py
index 183025fac627a66dcfc862e09892f1ff5b74b9da..ce5a21c5d023947ae0d4a3dbf119f6c68cee5f7a 100644
--- a/src/psimpy/inference/bayes_inference.py
+++ b/src/psimpy/inference/bayes_inference.py
@@ -10,7 +10,7 @@ from psimpy.utility import check_bounds
 
 _min_float = 10**(sys.float_info.min_10_exp)
 _max_float = 10**(sys.float_info.max_10_exp)
-_max_e = np.log(np.array(_max_float, dtype=np.float128))
+_max_e = np.log(np.array(_max_float, dtype=float))
 
 class BayesInferenceBase(ABC):
 
diff --git a/tests/test_run_mass_point_model.py b/tests/test_run_mass_point_model.py
index 26d5a79580d3c96dc8b48fad37f605ff13957ac7..a406be5be1c758c0b3d591d6d22733cb0b1dd118 100644
--- a/tests/test_run_mass_point_model.py
+++ b/tests/test_run_mass_point_model.py
@@ -47,6 +47,6 @@ def test_run_mass_point_model():
 
     assert len(serial_output) == len(parallel_output)
     for i in range(len(serial_output)):
-        assert np.array_equal(serial_output[i], parallel_output[i])
+        assert np.max(np.abs(serial_output[i] - parallel_output[i])) < 1e-5
     
     assert serial_time > parallel_time