Add __array_ufunc__ to handle ufuncs functions.

__array_function__ and @implements decorator are not necessary but will probably be useful later with the numpy 1.16 protocol (optional in version 1.16 but enabled by default in 1.17).

https://numpy.org/doc/stable/reference/arrays.classes.html?highlight=array_function#numpy.class.__array_function__