Time for action – comparing using maxulp of 2

Let's do the same comparisons as in the previous Time for action section, but specify a maxulp of 2 when necessary:

  1. Determine the machine epsilon with the finfo() function:
    eps = np.finfo(float).eps
print("EPS", eps)

    The epsilon would be as follows:

    EPS 2.22044604925e-16
  2. Do the comparisons as done in the previous Time for action section, but use the assert_array_max_ulp() function with the appropriate maxulp value:
    print("1", np.testing.assert_array_max_ulp(1.0, 1.0 + eps))
print("2", np.testing.assert_array_max_ulp(1.0, 1 + 2 * eps, maxulp=2))

    The output is as follows:

    1 1.0
2 2.0

What just happened?

We compared the same values as the previous Time for action section, but specified a maxulp of 2 in the second comparison. Using the assert_array_max_ulp() function with the appropriate maxulp value, these tests passed with a return value of the number of ULPs.

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset