Thank you all kindly for your responses! Based on your encouragement, I
will pursue an ndarray subclass / __array_ufunc__ implementation. I had
been toying with np.set_numeric_ops, which is less than ideal (for example,
np.ndarray.around segfaults if I use set_numeric_ops in any way).
A second que
Right before 1.12, I arranged an API around an np.ndarray subclass, making
use of __array_ufunc__ to customize behavior based on structured dtype (we
come from c++ and really like operator overloading). Having seen
__array_ufunc__ featured in Travis Oliphant's Guide to NumPy: 2nd Edition,
I assumed