Re: [Numpy-discussion] unable to build atlas and lapack correctly

Alexandre Lacoste wrote: > > Thanks for your answer > > > You can simply use the sse2 version of atlas: SSE2 is available on any > x86_64 CPU. > > > Are you sure ? I can't find it... > http://packages.ubuntu.com/search?arch=amd64&searchon=names&keywords=atlas%20sse2 >

Re: [Numpy-discussion] help to speed up the python code

frank wang wrote: > Hi, Bob, > > The problem is that I want to resample my data with another sampling > rate. the two rates is very close. I use the formula: > > s(t)=sum(a_k*sinc(t-kTs)). > > the new sampling rate is Ts', so I have > s(nTs')=sum(a_k*sinc(nTs'-kTs)). The sum index k is over t

Re: [Numpy-discussion] help to speed up the python code

frank wang wrote: > Hi, Bob, > > The problem is that I want to resample my data with another sampling > rate. the two rates is very close. I use the formula: > > s(t)=sum(a_k*sinc(t-kTs)). > > the new sampling rate is Ts', so I have > s(nTs')=sum(a_k*sinc(nTs'-kTs)). The sum index k is over t

Re: [Numpy-discussion] help to speed up the python code

Hi, Bob, The problem is that I want to resample my data with another sampling rate. the two rates is very close. I use the formula: s(t)=sum(a_k*sinc(t-kTs)). the new sampling rate is Ts', so I have s(nTs')=sum(a_k*sinc(nTs'-kTs)). The sum index k is over the (-P, P), Centered at n. The n

Re: [Numpy-discussion] unable to build atlas and lapack correctly

Thanks for your answer > You can simply use the sse2 version of atlas: SSE2 is available on any > x86_64 CPU. > Are you sure ? I can't find it... http://packages.ubuntu.com/search?arch=amd64&searchon=names&keywords=atlas%20sse2 > > Note that on Ubuntu, you can link to libatals3fg-base (the gfo

Re: [Numpy-discussion] (no subject)

frank wang wrote: > Hi, > > In my work, I want to implement a fir filter with an input array. > Since performing the filter on each input sample is slow, are there > fast way to perform the fir filter operation? Are there ways to > convert input into an array and perform the array multipication?

[Numpy-discussion] (no subject)

Hi, In my work, I want to implement a fir filter with an input array. Since performing the filter on each input sample is slow, are there fast way to perform the fir filter operation? Are there ways to convert input into an array and perform the array multipication? Thanks Frank

Re: [Numpy-discussion] unable to build atlas and lapack correctly

Alexandre Lacoste wrote: > Hi I'm trying to install numpy on a x86_64 ubuntu hardy and I'm having > a hard time to get atlas to be linked correctly. > > Installing it through apt-get install with the astraw repository > worked on my intel dual core. But now I need to install it in on a > x86_64 and

[Numpy-discussion] unable to build atlas and lapack correctly

Hi I'm trying to install numpy on a x86_64 ubuntu hardy and I'm having a hard time to get atlas to be linked correctly. Installing it through apt-get install with the astraw repository worked on my intel dual core. But now I need to install it in on a x86_64 and the astraw binaries will link again

Re: [Numpy-discussion] [Fwd: Large array using 4 times as much memory as it should]

On Thu, Oct 30, 2008 at 6:48 PM, Rick Giuly <[EMAIL PROTECTED]> wrote: > Please disregard the last message. > > I just realised 8*4=32, very sorry about this. Apparently I need some sleep. > > -Rick > > Original Message > Subject: Large array using 4 times as much memory as it sho

[Numpy-discussion] [Fwd: Large array using 4 times as much memory as it should]

Please disregard the last message. I just realised 8*4=32, very sorry about this. Apparently I need some sleep. -Rick Original Message Subject: Large array using 4 times as much memory as it should Date: Thu, 30 Oct 2008 18:41:44 -0700 From: Rick Giuly <[EMAIL PROTECTED]> To: n

Re: [Numpy-discussion] Large array using 4 times as much memory as it should

On Thu, Oct 30, 2008 at 20:41, Rick Giuly <[EMAIL PROTECTED]> wrote: > Hello All, > > I find that python is using about four times as much memory as it should > need for arrays. This is problematic as I need to use all available > memory for large 3D imaging datasets. Is there a way to get around t

[Numpy-discussion] Large array using 4 times as much memory as it should

Hello All, I find that python is using about four times as much memory as it should need for arrays. This is problematic as I need to use all available memory for large 3D imaging datasets. Is there a way to get around this problem? Am I making a mistake? Is it a bug? (I'm running windowsXP 32

Re: [Numpy-discussion] Numpy matrix multiplication slow even though ATLAS linked

On Thursday 30 October 2008 22:32:22 Robert Kern wrote: > On Thu, Oct 30, 2008 at 17:19, Jan-Willem van de Meent > <[EMAIL PROTECTED]> wrote: > > I posted a full build log earlier today, but I think it is still awaiting > > approval of moderators because of file size limitations. > > There are no m

Re: [Numpy-discussion] Numpy matrix multiplication slow even though ATLAS linked

On Thu, Oct 30, 2008 at 17:19, Jan-Willem van de Meent <[EMAIL PROTECTED]> wrote: > I posted a full build log earlier today, but I think it is still awaiting > approval of moderators because of file size limitations. There are no moderators to approve it. The file size limitation is just a limitat

Re: [Numpy-discussion] Numpy matrix multiplication slow even though ATLAS linked

On Thursday 30 October 2008 18:41:51 Charles R Harris wrote: > On Thu, Oct 30, 2008 at 5:19 AM, Jan-Willem van de Meent < > > [EMAIL PROTECTED]> wrote: > > Dear all, > > > > This is my first post to this list. I am having perfomance issues with > > with numpy/atlas. Doing dot(a,a) for a 2000x2000 m

[Numpy-discussion] Complete LAPACK needed

Dear all, I need to use functions in scipy which depend on having a complete lapack library. However, I am having a bit of trouble installing numpy and referencing a complete lapack library that I built. I have a few questions that I am hoping someone can help me answer: 1) This machine is a cl

Re: [Numpy-discussion] Numpy matrix multiplication slow even though ATLAS linked

On Thu, Oct 30, 2008 at 5:19 AM, Jan-Willem van de Meent < [EMAIL PROTECTED]> wrote: > Dear all, > > This is my first post to this list. I am having perfomance issues with with > numpy/atlas. Doing dot(a,a) for a 2000x2000 matrix takes about 1m40s, even > though numpy is appears to link to my atl

Re: [Numpy-discussion] passing a C array to embedded Python fromC code

Hi Chris, > Matthieu Brucher wrote: > > If you can > > follow a French tutorial, you can go on > > > http://matthieu-brucher.developpez.com/tutoriels/python/swig-numpy/#LV > > to have a skeletton for your issue. > > That looks very useful -- any chance of an English > translation? My one > ye

Re: [Numpy-discussion] passing a C array to embedded Python from C code

2008/10/30 Christopher Barker <[EMAIL PROTECTED]>: > Matthieu Brucher wrote: >> If you can >> follow a French tutorial, you can go on >> http://matthieu-brucher.developpez.com/tutoriels/python/swig-numpy/#LV >> to have a skeletton for your issue. > > That looks very useful -- any chance of an Engl

Re: [Numpy-discussion] passing a C array to embedded Python from C code

Matthieu Brucher wrote: > If you can > follow a French tutorial, you can go on > http://matthieu-brucher.developpez.com/tutoriels/python/swig-numpy/#LV > to have a skeletton for your issue. That looks very useful -- any chance of an English translation? My one year of high school French is provi

Re: [Numpy-discussion] passing a C array to embedded Python from C code

A few comments: Chris LeBlanc wrote: > I'm thinking I could create a > NumPy array that uses the same C array (1d, 2d, or 3d) that our > program is using instead of copying the memory. yes, you can do that. > At this point, I'm thinking the python code wouldn't need to return > any objects bec

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

Francesc Alted wrote: > A Thursday 30 October 2008, Robert Kern escrigué: > [clip] > OTOH, now that I think about it, I don't think there is really any coherent way to mix field selection with any other indexing operations. At least, not within the same brackets. Hmm. So maybe >>>

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

A Thursday 30 October 2008, Robert Kern escrigué: [clip] > >> OTOH, now that I think about it, I don't think there is really any > >> coherent way to mix field selection with any other indexing > >> operations. At least, not within the same brackets. Hmm. So maybe > >> the link to fancy indexing ca

Re: [Numpy-discussion] Automatic differentiation (was Re: second-order gradient)

> In your experience, is this functionality enough to start a separate > package, or should we try to include it somewhere else? Otherwise we > could think of a new SciKit. I confess to knowing no details about scikits so I don't know what the difference really is between a "new package" and a sc

Re: [Numpy-discussion] Automatic differentiation (was Re: second-order gradient)

2008/10/30 Rob Clewley <[EMAIL PROTECTED]>: >> http://en.wikipedia.org/wiki/Automatic_differentiation > > That would be great, but wouldn't that be best suited as a utility > requiring Sympy? You'll want to take advantage of all sorts of > symbolic classes, especially for any source code transforma

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

On Thu, Oct 30, 2008 at 04:33, Francesc Alted <[EMAIL PROTECTED]> wrote: > A Thursday 30 October 2008, Robert Kern escrigué: >> On Wed, Oct 29, 2008 at 19:05, Travis E. Oliphant >> >> <[EMAIL PROTECTED]> wrote: >> > Hi all, >> > >> > I'd like to add to NumPy the ability to clone a data-type object

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

On Thu, Oct 30, 2008 at 08:27, Travis E. Oliphant <[EMAIL PROTECTED]> wrote: > >> I'm not sure what this accomplishes. Would the dummy fields that fill >> in the space be inaccessible? E.g. tuple(subarr[i,j,k]) gives a tuple >> with no numpy.void scalars? That would be a novel feature, but I'm not

Re: [Numpy-discussion] any interest in including asecond-ordergradient?

> 2008/10/29 Fernando Perez <[EMAIL PROTECTED]>: >> I think it's fine to ask for functions that compute higher order >> derivatives of n-d arrays: we already have diff(), which operates on a >> single direction, and a hessian could make sense (with the caveats >> David points out). But with highe

[Numpy-discussion] Automatic differentiation (was Re: second-order gradient)

> Maybe we should focus on writing a decent 'deriv' function then. I > know Konrad Hinsen's Scientific had a derivatives package > (Scientific.Functions.Derivatives) that implemented automatic > differentiation: > > http://en.wikipedia.org/wiki/Automatic_differentiation That would be great, but w

Re: [Numpy-discussion] numpy array change notifier?

On Oct 30 10:06 -0400, Pierre GM wrote: > Steve, > Right, there's not a lot of visibility for this one. Yes, wouldn't it be better placed in the Cookbook section (or any other suitable place that's visible to anyone entering scipy.org)? Other pages are unfortunately also not visible, like http:

Re: [Numpy-discussion] numpy array change notifier?

Steve, Right, there's not a lot of visibility for this one. You could try google 'numpy subclass'. there's also an entry on subclassing in the numpy user guide (through docs.scipy.org). Cheers P. On Thursday 30 October 2008 09:57:03 Steve Schmerler wrote: > On Oct 27 16:43 -0400, Pierre GM wrot

Re: [Numpy-discussion] numpy array change notifier?

On Oct 27 16:43 -0400, Pierre GM wrote: > > Erik, may be you could try the trick presented here : > http://www.scipy.org/Subclasses > in the __array_wrap__ section. Stupid question: How do I find pages like http://www.scipy.org/Subclasses on scipy.org? I can find them with the search function,

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

> I'm not sure what this accomplishes. Would the dummy fields that fill > in the space be inaccessible? E.g. tuple(subarr[i,j,k]) gives a tuple > with no numpy.void scalars? That would be a novel feature, but I'm not > sure it fits the problem. On the contrary: > Yes, that was the idea. You

Re: [Numpy-discussion] numpy array change notifier?

Erik Tollerud gmail.com> writes: > > Is there any straightforward way of notifying on change of a numpy > array that leaves the numpy arrays still efficient? > You can instantiate the following ndarray derived class with a callable argument. Any change to the data will call the callback funct

[Numpy-discussion] Numpy matrix multiplication slow even though ATLAS linked

Dear all, This is my first post to this list. I am having perfomance issues with with numpy/atlas. Doing dot(a,a) for a 2000x2000 matrix takes about 1m40s, even though numpy is appears to link to my atlas libraries: I did a quick benchmark by running the following script: #! /usr/bin/env

Re: [Numpy-discussion] Adding the ability to "clone" a few fields from a data-type

A Thursday 30 October 2008, Robert Kern escrigué: > On Wed, Oct 29, 2008 at 19:05, Travis E. Oliphant > > <[EMAIL PROTECTED]> wrote: > > Hi all, > > > > I'd like to add to NumPy the ability to clone a data-type object so > > that only a view fields are copied over but that it retains the > > same t

Re: [Numpy-discussion] passing a C array to embedded Python from C code

> Does this approach make sense? Is there a better way to go about it? > Maybe calling a custom module from the Python code that does the C > array to NumPy translation using Cython/pyrex/swig/etc. Would it be > possible to use the same C arrays from here without copying them? Hi, Your case see