Exception when checking if two bfloat16 arrays are close

[UtkarshKunwar]

It fails on doing np.isclose(arr, arr) but works properly when doing np.isclose(arr, arr.astype(np.float64))

[sousoux]

This is intrinsic to the type system in numpy and there is no way around it currently.