incrementing along a diagonal

Bill Baxter wbaxter at gmail.com
Wed Oct 11 18:58:45 CDT 2006


On 10/12/06, David Novakovic <david at distip.com.au> 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 at 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




More information about the Numpy-discussion mailing list