On Oct 30, 2007 11:40 PM, Alan G Isaac <[EMAIL PROTECTED]> wrote: > >>> 1.0**numpy.array([1,2,3]) > array([ 1., 1., 1.]) > >>> 1.0**numpy.mat([1,2,3]) > Traceback (most recent call last): > File "<stdin>", line 1, in <module> > TypeError: unsupported operand type(s) for ** or pow(): 'float' and 'matrix' > > Why the restriction for matrices? > > Same question for matrices conformable for broadcasting: > why not produce the equivalent of the array result?
With the caveat that I don't use the matrix class: I believe it's because M**n results in the matrix power of n. It would be confusing if n**M did a broadcast element wise power. -- . __ . |-\ . . [EMAIL PROTECTED] _______________________________________________ Numpy-discussion mailing list [email protected] http://projects.scipy.org/mailman/listinfo/numpy-discussion
