[Numpy-discussion] untenable matrix behavior in SVN

[Anne Archibald]

On 29/04/2008, Gael Varoquaux <gael.varoquaux at normalesup.org> wrote:
> On Tue, Apr 29, 2008 at 11:03:58PM +0200, Anne Archibald wrote:
>  > I am puzzled by this. What is the rationale for x[i,:] not being a 1-d
>  > object?
>
> It breaks A*B[i, :] where A and B are matrices.

Really? How?

In [26]: A = np.matrix([[1,0],[0,1],[1,1]])

In [28]: A*np.ones(2)
Out[28]: matrix([[ 1.,  1.,  2.]])

In [29]: np.ones(3)*A
Out[29]: matrix([[ 2.,  2.]])


I guess you don't get dot products with *:

In [30]: np.ones(3).T*np.ones(3)
Out[30]: array([ 1.,  1.,  1.])

Since the whole point of matrices it to avoid having to type "dot" I
guess this should be convincing.


Anne