[Numpy-discussion] untenable matrix behavior in SVN

[Charles R Harris]

On Tue, Apr 29, 2008 at 12:20 AM, Alan G Isaac <aisaac at american.edu> wrote:

> On Mon, 28 Apr 2008, Christopher Barker apparently wrote:
> > I'm going to try to put down what I think are the key,
> > very simple, questions:
>
> For useful reference, these are now here:
> <URL:http://www.scipy.org/MatrixIndexing#guidelines>
>
> Cheers,
> Alan
>
>
May I add that if I edit defmatrix.py to act like an array for scalar
indexing, then the following works.

In [1]: a = matrix(eye(2))

In [2]: array([a,a])
Out[2]:
array([[[ 1.,  0.],
        [ 0.,  1.]],

       [[ 1.,  0.],
        [ 0.,  1.]]])

This generates an error with the current version of matrix and, frankly, I
am not going to be bothered going all through the numpy c sources to special
case matrices to fix that. Someone else can do it if they wish. There are
recursive routines that expect the dimensions to decrease on each call.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20080429/579116b6/attachment.html>