import numpy as np eps = np.finfo(float).eps print "EPS", eps print "1", np.testing.assert_array_almost_equal_nulp(1.0, 1.0 + eps) print "2", np.testing.assert_array_almost_equal_nulp(1.0, 1.0 + 2 * eps)