On 10/12/06, David Novakovic <[EMAIL PROTECTED]> wrote: > Johannes Loehnert wrote: > This is very nice, exactly what i want, but it doesnt work for mxn > matricies: > > >>> x = zeros((5,3)) > >>> x > array([[0, 0, 0], > [0, 0, 0], > [0, 0, 0], > [0, 0, 0], > [0, 0, 0]]) > >>> index = arange(min(x.shape[0],x.shape[1])) > >>> x[index,index] += 1 > >>> x > array([[1, 0, 0], > [0, 1, 0], > [0, 0, 1], > [0, 0, 0], > [0, 0, 0]])
Exactly what output are you expecting? That is the definition of the 'diagonal' for a non-square matrix. If you're expecting something else then what you want is not the diagonal. > Just for reference, this is the line of perl i'm trying to port: > > > like: > > for index in diag_iter(matrix,*axes): > matrix[index] +=1 That's not going to change the mathematical definition of the diagonal of a non-square matrix. > PS: If anyone would care to link me to the subscription page for the > mailing list so you dont have to CC me all the time :) Check the bottom of this message. > _______________________________________________ > Numpy-discussion mailing list > Numpy-discussion@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > --bb ------------------------------------------------------------------------- Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion