[Numpy-discussion] numpy ndarray questions

[Jochen]

On Mon, 2009-01-26 at 19:25 -0600, Ryan May wrote:
> Jochen wrote:
> > Hi all,
> > 
> > I just wrote ctypes bindings to fftw3 (see
> > http://projects.scipy.org/pipermail/scipy-user/2009-January/019557.html
> > for the post to scipy). 
> > Now I have a couple of numpy related questions:
> > 
> > In order to be able to use simd instructions I 
> > create an ndarray subclass, which uses fftw_malloc to allocate the
> > memory and fftw_free to free the memory when the array is deleted. This
> > works fine for inplace operations however if someone does something like
> > this:
> > 
> > a = fftw3.AlignedArray(1024,complex)
> > 
> > a = a+1
> > 
> > a.ctypes.data points to a different memory location (this is actually an
> > even bigger problem when executing fftw plans), however 
> > type(a) still gives me <class 'fftw3.planning.AlignedArray'>.
> 
> This might help some:
> 
> http://www.scipy.org/Subclasses
> 
> Ryan
> 
Thanks, I had read about __array_finalize__, but not about
__array_wrap__. I'm not quite sure if I understand how I can use this to
get around my problem, can you explain?

Cheers
Jochen