Re: [Numpy-discussion] Strange numpy behaviour (bug?)

2012-01-17 Thread Sturla Molden
Never mind this, it was my own mistake as I expected :-) def __chunk(n,size): x = range(0,n,size) x.append(n) return zip(x[:-1],x[1:]) makes it a lot better :) Sturla Den 18.01.2012 06:26, skrev Sturla Molden: > While "playing" with a point-in-polygon test, I have discovered s

[Numpy-discussion] Strange numpy behaviour (bug?)

2012-01-17 Thread Sturla Molden
While "playing" with a point-in-polygon test, I have discovered some a failure mode that I cannot make sence of. The algorithm is vectorized for NumPy from a C and Python implementation I found on the net (see links below). It is written to process a large dataset in chunks. I'm rather happy w