[Numpy-discussion] Suggestion for a speed improved accumarray receipe

2013-02-17 Thread Michael Löffler
Hello numpy list! I had tried to contribute to the AccumarrayLike recipe cookbook, but access seems to be restricted and no registration allowed. The current recipe mimics the matlab version in most aspects, but concerning performance it's horribly slow. The following snippet handles only the mo

Re: [Numpy-discussion] calling NumPy from Julia - a plea for fewer macros

2013-02-17 Thread Steven G. Johnson
Nathaniel Smith wrote: > There are currently no plans to produce a NumPy 2.0, but everything you > suggest would be just fine as changes to numpy 1.x. PRs gratefully accepted. Thanks, just posted https://github.com/numpy/numpy/issues/2997 https://github.com/numpy/numpy/issues/299

Re: [Numpy-discussion] calling NumPy from Julia - a plea for fewer macros

2013-02-17 Thread Charles R Harris
On Sun, Feb 17, 2013 at 12:43 PM, Charles R Harris < charlesr.har...@gmail.com> wrote: > > > On Sun, Feb 17, 2013 at 9:12 AM, Steven G. Johnson > wrote: > >> Dear NumPy developers, >> >> I've been working on a glue package that allows the Julia language >> (http://julialang.org/) to call Python r

Re: [Numpy-discussion] calling NumPy from Julia - a plea for fewer macros

2013-02-17 Thread Charles R Harris
On Sun, Feb 17, 2013 at 9:12 AM, Steven G. Johnson wrote: > Dear NumPy developers, > > I've been working on a glue package that allows the Julia language > (http://julialang.org/) to call Python routines easily > https://github.com/stevengj/PyCall.jl > and I'm using NumPy to pass multidime

Re: [Numpy-discussion] FFTW bindings now implement numpy.fft interface

2013-02-17 Thread Henry Gomersall
On Sun, 2013-02-17 at 12:38 -0500, Neal Becker wrote: > The 1st example says: > >>> import pyfftw > >>> import numpy > >>> a = pyfftw.n_byte_align_empty(128, 16, 'complex128') > >>> a[:] = numpy.random.randn(128) + 1j*numpy.random.randn(128) > >>> b = pyfftw.interfaces.numpy_fft.fft(a) > > I don't

Re: [Numpy-discussion] FFTW bindings now implement numpy.fft interface

2013-02-17 Thread Neal Becker
Henry Gomersall wrote: > Some of you may be interested in the latest release of my FFTW bindings. > It can now serve as a drop in replacement* for numpy.fft and > scipy.fftpack. > > This means you can get most of the speed-up of FFTW with a one line code > change or monkey patch existing librarie

Re: [Numpy-discussion] calling NumPy from Julia - a plea for fewer macros

2013-02-17 Thread Nathaniel Smith
On 17 Feb 2013 08:13, "Steven G. Johnson" wrote: > Julia has the ability to call C functions directly (without writing C > glue), and I've been exploiting this to write PyCall purely in Julia. > (This is nice for a number of reasons; besides programming and linking > convenience, it means that I ca

[Numpy-discussion] calling NumPy from Julia - a plea for fewer macros

2013-02-17 Thread Steven G. Johnson
Dear NumPy developers, I've been working on a glue package that allows the Julia language (http://julialang.org/) to call Python routines easily https://github.com/stevengj/PyCall.jl and I'm using NumPy to pass multidimensional arrays between languages. Julia has the ability to call C fu

[Numpy-discussion] FFTW bindings now implement numpy.fft interface

2013-02-17 Thread Henry Gomersall
Some of you may be interested in the latest release of my FFTW bindings. It can now serve as a drop in replacement* for numpy.fft and scipy.fftpack. This means you can get most of the speed-up of FFTW with a one line code change or monkey patch existing libraries. Lots of other goodness too of co

Re: [Numpy-discussion] np.percentile docstring

2013-02-17 Thread Sebastian Berg
Hey, On Sun, 2013-02-17 at 11:50 +0100, Andreas Hilboll wrote: > In my numpy 1.6.1 (from Ubuntu 12.04LTS repository), the docstring of > np.percentile is wrong. I'm not just submitting a PR because I don't > understand something. > > in the "Notes" and "Examples" sections, there seems to be some

[Numpy-discussion] np.percentile docstring

2013-02-17 Thread Andreas Hilboll
In my numpy 1.6.1 (from Ubuntu 12.04LTS repository), the docstring of np.percentile is wrong. I'm not just submitting a PR because I don't understand something. in the "Notes" and "Examples" sections, there seems to be some confusion if ``q`` should be in the range [0,1] or in [0,100]. The paramet