Hi,
On 16 August 2012 19:38, Sebastian Berg sebast...@sipsolutions.net wrote:
Hello,
Just to throw it in, if you do not mind useing scipy, you can use its
multidimensional correlate method instead:
Well, I don't need to count the number of neighbours as in the linked
example. What I wanted
Hi,
I've just come across Travis Oliphant's array version of the game of life
here https://gist.github.com/3353411. It's a really nice example of how
to efficiently find neighbours using numpy and family. However, I wanted to
apply these ideas to higher order arrays (eg GRID being three
Hello,
Just to throw it in, if you do not mind useing scipy, you can use its
multidimensional correlate method instead:
stamp = np.ones((3,3,3))
stamp[1,1,1] = 0
num_neighbours = scipy.ndimage.correlate(x, stamp, mode='wrap'))
In the link np.roll is used to implement periodic boundaries
I have a pull request for a neighborhood function at
https://github.com/numpy/numpy/pull/303 . I think IMHO it handles these
problems quite handily. It does rely on my pad routine that is in Numpy
1.7, so you would need to get the 1.7 beta installed or install the
development branch.
For your