Switch to numexpr3 when / if available
Numexpr3 adds many new functions and a more general, pythonic syntax. It would be great to move to it. For example, bitmasks can be handled properly: https://github.com/pydata/numexpr/issues/242
It's not clear to me the status of numexpr3 is though.