Fix numpy's version of power operator

Test & implementation