Comment by dumah

This is because of the dispatch semantics for dunder operator methods, which the Python docs are horrible at communicating.

The types returned as you described are exactly what I’d expect and idiomatic for Python. This isn’t a problem with numpy.