On Tue, Apr 29, 2008 at 12:20 AM, Alan G Isaac <[EMAIL PROTECTED]> 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
_______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
