> On Tue, 2009-01-27 at 14:16 +0100, Sturla Molden wrote:
>> def aligned_array(N, dtype):
>> d = dtype()
>> tmp = numpy.array(N * d.nbytes + 16, dtype=numpy.uint8)
>> address = tmp.__array_interface__['data'][0]
>> offset = (16 - address % 16) % 16
>> return tmp[offset:off
On Tue, 2009-01-27 at 14:16 +0100, Sturla Molden wrote:
> On 1/27/2009 12:37 PM, Sturla Molden wrote:
>
> > import ctypes
> > import numpy
> >
> > fftw_malloc = ctypes.cdll.fftw.fftw_malloc
> > fftw_malloc.argtypes = [ctypes.c_ulong,]
> > fftw_malloc.restype = ctypes.c_ulong
> >
> > def aligned_
On Tue, 2009-01-27 at 12:37 +0100, Sturla Molden wrote:
> On 1/27/2009 1:26 AM, Jochen wrote:
>
> > a = fftw3.AlignedArray(1024,complex)
> >
> > a = a+1
>
> = used this way is not assignment, it is name binding.
>
> It is easy to use function's like fftw_malloc with NumPy:
>
>
> import ctypes
On 1/27/2009 12:37 PM, Sturla Molden wrote:
> import ctypes
> import numpy
>
> fftw_malloc = ctypes.cdll.fftw.fftw_malloc
> fftw_malloc.argtypes = [ctypes.c_ulong,]
> fftw_malloc.restype = ctypes.c_ulong
>
> def aligned_array(N, dtype):
> d = dtype()
> address = fftw_malloc(N * d.nbyte
On 1/27/2009 12:37 PM, Sturla Molden wrote:
> It is easy to use function's like fftw_malloc with NumPy:
Besides this, if I were to write a wrapper for FFTW in Python, I would
consider wrapping FFTW's Fortran interface with f2py.
It is probably safer, as well as faster, than using ctypes. It wou
On 1/27/2009 12:37 PM, Sturla Molden wrote:
> address = fftw_malloc(N * d.nbytes) # restype = ctypes.c_ulong
> if (address = 0):
if (address == ): raise MemoryError, 'fftw_malloc returned NULL'
Sorry for the typo.
S.M.
___
Numpy-d
On 1/27/2009 1:26 AM, Jochen wrote:
> a = fftw3.AlignedArray(1024,complex)
>
> a = a+1
= used this way is not assignment, it is name binding.
It is easy to use function's like fftw_malloc with NumPy:
import ctypes
import numpy
fftw_malloc = ctypes.cdll.fftw.fftw_malloc
fftw_malloc.argtypes =
On 1/27/2009 6:03 AM, Jochen wrote:
> BTW memmap arrays have
> the same problem
> if I create a memmap array and later do something like
> a=a+1
> all later changes will not be written to the file.
= is Python's rebinding operator.
a = a + 1 rebinds a to a different object.
As for ndarray's,
On Tue, 2009-01-27 at 14:46 +0900, David Cournapeau wrote:
> Jochen wrote:
> > On Tue, 2009-01-27 at 13:54 +0900, David Cournapeau wrote:
> >
> >> Jochen wrote:
> >>
> >>> On Tue, 2009-01-27 at 13:28 +0900, David Cournapeau wrote:
> >>>
> >>>
> Jochen wrote:
>
> >>>
Jochen wrote:
> On Tue, 2009-01-27 at 13:54 +0900, David Cournapeau wrote:
>
>> Jochen wrote:
>>
>>> On Tue, 2009-01-27 at 13:28 +0900, David Cournapeau wrote:
>>>
>>>
Jochen wrote:
> On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau wrote:
>
On Tue, 2009-01-27 at 13:54 +0900, David Cournapeau wrote:
> Jochen wrote:
> > On Tue, 2009-01-27 at 13:28 +0900, David Cournapeau wrote:
> >
> >> Jochen wrote:
> >>
> >>> On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau wrote:
> >>>
> >>>
> Jochen wrote:
>
> >>>
Jochen wrote:
> On Tue, 2009-01-27 at 13:28 +0900, David Cournapeau wrote:
>
>> Jochen wrote:
>>
>>> On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau wrote:
>>>
>>>
Jochen wrote:
> Hi all,
>
> I just wrote ctypes bindings to fftw3 (see
>
On Tue, 2009-01-27 at 13:28 +0900, David Cournapeau wrote:
> Jochen wrote:
> > On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau wrote:
> >
> >> Jochen wrote:
> >>
> >>> Hi all,
> >>>
> >>> I just wrote ctypes bindings to fftw3 (see
> >>> http://projects.scipy.org/pipermail/scipy-user/200
Jochen wrote:
> On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau 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
On Tue, 2009-01-27 at 12:49 +0900, David Cournapeau 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:
> >
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 subcl
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
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 sub
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
19 matches
Mail list logo
