[Numpy-discussion] Is numpy.convolve based on LAPACK routine?
braingateway
braingateway at gmail.com
Thu Nov 11 03:32:44 EST 2010
Pauli Virtanen :
> On Thu, 11 Nov 2010 00:43:32 +0100, LittleBigBrain wrote:
>
>> I am wondering, is numpy.convolve based on LAPACK routine? Can it be
>> speedup by using ATLAS?
>>
>
> LAPACK and Atlas do not AFAIK have convolution routines -- that's not
> linear algebra. MKL on the other hand would have some. The implementation
> in Numpy is the straightforward one, without SIMD etc.
>
> For large datasets, scipy.signal.fftconvolve should be faster.
>
>
Thanks for the point. The fftconvolve is only fast when two input are
both long enough say >200, other wise pure convovlution will be faster.
And fftconvolve will take too more RAM than normal convolution in this
case. I were using MKL which does have convolution routines, so I
thought might ATLAS also have it.
Thanks a lot!
LittleBigBrain
More information about the NumPy-Discussion
mailing list