Re: [Numpy-discussion] Recurrence relationships

2009-05-06 Thread Talbot, Gerry
To: Discussion of Numerical Python Subject: Re: [Numpy-discussion] Recurrence relationships On 5/6/2009 10:00 AM Talbot, Gerry apparently wrote: > for n in xrange(1,N): > y[n] = A*x[n] + B*y[n-1] So, x is known before you start? How big is N? Also, is y.shape (N,)? Do you need al

Re: [Numpy-discussion] Recurrence relationships

2009-05-06 Thread Talbot, Gerry
lf Of Keith Goodman Sent: Wednesday, May 06, 2009 9:53 AM To: Discussion of Numerical Python Subject: Re: [Numpy-discussion] Recurrence relationships On Wed, May 6, 2009 at 6:44 AM, Talbot, Gerry wrote: > Does anyone know how to efficiently implement a recurrence relationship in > nu

[Numpy-discussion] Recurrence relationships

2009-05-06 Thread Talbot, Gerry
Does anyone know how to efficiently implement a recurrence relationship in numpy such as: y[n] = A*x[n] + B*y[n-1] Thanks, Gerry ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://mail.scipy.org/mailman/listinfo/

Re: [Numpy-discussion] Indexing

2009-04-17 Thread Talbot, Gerry
/4/17 Talbot, Gerry : > I stumbled across this anomaly, the assignment to ncycle[i][b] fails, yet > ncycle[:][b] works and the final print works.  Is this a known bug or > limitation ?  I am running numpy-1.2.1 with python 2.5.4 The fancy indexing (indexing with an array) in ncycle[i] mak

[Numpy-discussion] Indexing

2009-04-17 Thread Talbot, Gerry
Hi, I stumbled across this anomaly, the assignment to ncycle[i][b] fails, yet ncycle[:][b] works and the final print works. Is this a known bug or limitation ? I am running numpy-1.2.1 with python 2.5.4 Gerry from numpy import * ncycle = arange(10) i = array([0, 1, 2, 3, 4, 5