code: def allclose_or_equiv(arr1, arr2, rtol=1e-5, atol=1e-8):
      
      arr1, arr2 = as_like_arrays(arr1, arr2)
      if arr1.shape != arr2.shape:
          return False
      return bool(
          isclose(arr1, arr2, rtol=rtol, atol=atol, equal_nan=True).all())
docs: Like np.allclose, but also allows values to be NaN in both arrays