[Numpy-discussion] numpy ndarray questions

[Ryan May]

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

-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma