SPAM-LOW: Re: incrementing along a diagonal
David Novakovic
david at distip.com.au
Wed Oct 11 19:35:35 CDT 2006
David Novakovic wrote:
> Thanks for the help, i've learnt a lot and also figured out something
> that does what I want, i'll paste an interactive session below:
>
> x = zeros((4,7))
> x
> array([[0, 0, 0, 0, 0, 0, 0],
> [0, 0, 0, 0, 0, 0, 0],
> [0, 0, 0, 0, 0, 0, 0],
> [0, 0, 0, 0, 0, 0, 0]])
> index = arange(min(x.shape[0], x.shape[1]))
> index2 = copy.deepcopy(index) #deep copy may be overkill
> for a,b in enumerate(index):
> ... index2[a] += a
>
Turns out this is not good at all. I guess I'm still open to suggestions
then :(
Dave
> ...
> if len(x[:,0]) > len(x[0]):
> ... x[index2,index] +=1
> ... else:
> ... x[index,index2] +=1
> ...
> x
> array([[1, 0, 0, 0, 0, 0, 0],
> [0, 0, 1, 0, 0, 0, 0],
> [0, 0, 0, 0, 1, 0, 0],
> [0, 0, 0, 0, 0, 0, 1]])
>
>
> Thanks for the tips
>
> Dave Novakovic
>
> P.S. subscribed to the list now
>
> Bill Baxter wrote:
>
>> Forgot to CC you...
>>
>> ---------- Forwarded message ----------
>> From: Bill Baxter <wbaxter at gmail.com>
>> Date: Oct 12, 2006 8:58 AM
>> Subject: Re: [Numpy-discussion] incrementing along a diagonal
>> To: Discussion of Numerical Python
>> <numpy-discussion at lists.sourceforge.net>
>>
>>
>> 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
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>
>
-------------------------------------------------------------------------
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