numpy uses: allclose(a, b, rtol=1e-05, atol=1e-08, equal_nan=False)