Re: [Numpy-discussion] Any interest in a 'heaviside' ufunc?

2015-02-03 Thread Aron Ahmadia
> This is not correct. The discrete form of the Heaviside step function has > the value 1 for x == 0. > Yeah, I was looking at it and wondering if I'd misremembered the definition. Assuming you're implementing the discrete Heaviside function, H[0] = 1 as Sturla notes. ___

Re: [Numpy-discussion] Any interest in a 'heaviside' ufunc?

2015-02-03 Thread Aron Ahmadia
That seems useful to me. On Tue, Feb 3, 2015 at 3:58 PM, Warren Weckesser wrote: > I have an implementation of the Heaviside function as numpy ufunc. Is > there any interest in adding this to numpy? The function is simply: > > 0if x < 0 > heaviside(x) = 0.5 if x =

Re: [Numpy-discussion] Sport Rock Friday?

2014-09-25 Thread Aron Ahmadia
Sorry for the SPAM folks. Itchy smartphone :) ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

[Numpy-discussion] Fwd: Sport Rock Friday?

2014-09-25 Thread Aron Ahmadia
-- Forwarded message -- From: *Brian Joseph* Date: Thursday, September 25, 2014 Subject: Sport Rock Friday? To: Aron Ahmadia Hey Aron want to hit up Sport Rock tomorrow night? Meg and I would like to go. Meg has to work until 8, so how does 8:30 sound

Re: [Numpy-discussion] Reporting a bug to Apple.

2014-06-10 Thread Aron Ahmadia
On Mon, Jun 9, 2014 at 8:52 PM, Sturla Molden wrote: > > If we report a segfault Apple will probably think we are to blame. 99.9 % > of bug reports come from idiots, and those who screen bug reports are > basically payed to hit the ignore button. Indeed, good luck on reporting > it. I've report

Re: [Numpy-discussion] 100 Numpy exercices

2014-05-27 Thread Aron Ahmadia
Very cool! On Tue, May 27, 2014 at 12:57 AM, Nicolas Rougier wrote: > > Hi all, > > I've updated the numpy exercices collection and made it available on > github at: > > https://github.com/rougier/numpy-100 > > > These exercices mainly comes from this mailing list and also from stack > overflow.

Re: [Numpy-discussion] High-quality memory profiling for numpy in python 3.5 / volunteers needed

2014-04-17 Thread Aron Ahmadia
code) for some notes on this. A On Thu, Apr 17, 2014 at 10:18 AM, Nathaniel Smith wrote: > On 17 Apr 2014 15:09, "Aron Ahmadia" wrote: > > > > > On the one hand it would be nice to actually know whether > posix_memalign is important, before making api decisions

Re: [Numpy-discussion] High-quality memory profiling for numpy in python 3.5 / volunteers needed

2014-04-17 Thread Aron Ahmadia
> On the one hand it would be nice to actually know whether posix_memalign is important, before making api decisions on this basis. FWIW: On the lightweight IBM cores that the extremely popular BlueGene machines were based on, accessing unaligned memory raised system faults. The default behavior

Re: [Numpy-discussion] Wiki page for building numerical stuff onWindows

2014-04-15 Thread Aron Ahmadia
Fernando is referring to this: http://lightningpython.org/ I have to admit that as a Pythonist supporting Windows builds, Cygwin and/or MinGW have always been the most appealing because so much of the rest of the numerical computing infrastructure is primarily supported on UNIXy platforms. A On

Re: [Numpy-discussion] Wiki page for building numerical stuff on Windows

2014-04-11 Thread Aron Ahmadia
Thanks Matthew for putting this page together. The OpenBLAS guys have been accepting/merging pull requests (their GitHub tree shows 26 contributors and no open pull requests), and I know that several people from the Python and Julia community have gotten pull requests merged. I modified your comm

Re: [Numpy-discussion] [RFC] should we argue for a matrix power operator, @@?

2014-03-17 Thread Aron Ahmadia
On Mon, Mar 17, 2014 at 7:53 AM, Nathaniel Smith wrote: > The thread so far, it sounds like the consensus answer is "meh, > whatever". So I'm thinking we should just drop @@ from the PEP, and if > it turns out that this is a problem we can always revisit it in the > ~3.6/3.7 timeframe. > +1 from

Re: [Numpy-discussion] It looks like Py 3.5 will include a dedicated infix matrix multiply operator

2014-03-14 Thread Aron Ahmadia
That's the best news I've had all week. Thanks for all your work on this Nathan. -A On Fri, Mar 14, 2014 at 8:51 PM, Nathaniel Smith wrote: > Well, that was fast. Guido says he'll accept the addition of '@' as an > infix operator for matrix multiplication, once some details are ironed > out:

Re: [Numpy-discussion] GSoC project: draft of proposal

2014-03-12 Thread Aron Ahmadia
Hi Leo, Out of curiosity, which vector math libraries did you have in mind as likely candidates for inclusion? How are you planning on selecting the library to integrate? Cheers, Aron On Wed, Mar 12, 2014 at 12:52 PM, Leo Mao wrote: > Hi, > The attachment is my draft of proposal. The project

Re: [Numpy-discussion] Code review request: PrintOptions

2013-10-27 Thread Aron Ahmadia
On Sun, Oct 27, 2013 at 9:42 PM, Neil Girdhar wrote: > Yeah, I realized that I missed that and figured it wouldn't matter since > it was my own master and I don't plan on making other changes to numpy. If > you don't mind, how do I move my changelist into a branch? I'm really > worried I'm goin

Re: [Numpy-discussion] RFC: is it worth giving a lightning talk at PyCon 2014 on numpy vbench-marking?

2013-10-15 Thread Aron Ahmadia
FWIW, I think it's a really cool project. I'm not the target demographic at PyCon, but I would be very surprised if this isn't a well-received talk. On Tue, Oct 15, 2013 at 12:13 PM, Yaroslav Halchenko wrote: > Hi Guys, > > PyCon 2014 will be just around the corner from where I am, so I decided

Re: [Numpy-discussion] Equivalent to IDL's help function

2013-10-07 Thread Aron Ahmadia
If you're using IPython (strongly recommended) a? a?? If you're just using Python help(a) A On Mon, Oct 7, 2013 at 1:15 PM, Siegfried Gonzi wrote: > Hi all > > What is the equivalent to IDL its help function, e.g. > > == > IDL> a = make_array(23,23,) > > IDL> help,a > > will result in: > > A

[Numpy-discussion] Handling warnings when solving nearly-singular matrices

2013-09-17 Thread Aron Ahmadia
Hey folks, The mathematicians over at Walking Randomly have expressed some dissatisfaction with the lack of warnings when solving nearly-singular matrices: http://www.walkingrandomly.com/?p=5092 An issue has been raised here: https://github.com/numpy/numpy/issues/3755 The Julia developers are a

Re: [Numpy-discussion] Correct way to query NumPy for linktime BLAS and LAPACK

2013-08-08 Thread Aron Ahmadia
BLAS project, > like different MKL version. Maybe just on the them will get imported if the > library name is the same. > > HTH > > Fred > > > On Mon, Aug 5, 2013 at 5:06 PM, Aron Ahmadia wrote: > >> Dear NumPy Developers, >> >> In the Clawpack/* reposi

[Numpy-discussion] Correct way to query NumPy for linktime BLAS and LAPACK

2013-08-05 Thread Aron Ahmadia
Dear NumPy Developers, In the Clawpack/* repositories [1], we use a mixture of Fortran and Python source, currently glued together using f2py. Occasionally, we'll need to link the Fortran code directly against LAPACK. In particular, we're using dgeev and dgesv to solve several different Riemann

Re: [Numpy-discussion] Memory leak

2013-06-13 Thread Aron Ahmadia
Hi Petro, What version of numpy are you running? A On Thu, Jun 13, 2013 at 3:50 PM, Pietro Bonfa' wrote: > Dear Numpy users, > > I have a memory leak in my code. A simple way to reproduce my problem is: > > import numpy > > class test(): > def __init__(self): > pass > > def t(s

Re: [Numpy-discussion] f2py and object libraries

2013-04-25 Thread Aron Ahmadia
Hi Jon, I have personally never used f2py to link to library code, but according to the documentation: http://cens.ioc.ee/projects/f2py2e/usersguide/index.html#command-f2py if you are building a module (that is, invoking f2py -c), then you can either include the path to the .a file directly at t

[Numpy-discussion] (@Pat Marion) Re: Yes, this one again "ImportError: No module named multiarray"

2013-03-12 Thread Aron Ahmadia
Pat Marion at Kitware did some work on this, I'm pinging him on the thread. A On Tue, Mar 12, 2013 at 2:05 PM, Dinesh B Vadhia wrote: > ** > Does that mean numpy won't work with freeze/create_executable type of > tools or is there a workaround? > > > *From:* Aro

Re: [Numpy-discussion] Yes, this one again "ImportError: No module named multiarray"

2013-03-12 Thread Aron Ahmadia
multiarray is an extension module that lives within numpy/core, that is, when, "import multiarray" is called, (and it's the first imported extension module in numpy), multiarray.ext (ext being dll on Windows I guess), gets dynamically loaded. "No module named multiarray" is indicating problems wit

Re: [Numpy-discussion] Remove interactive setup

2013-03-04 Thread Aron Ahmadia
I've built numpy on many different machines, including supercomputers, and I have never used interactive setup. I agree with the proposal to remove it. A On Mon, Mar 4, 2013 at 11:10 PM, Charles R Harris wrote: > In distutils there are three files that provide some interactive setup: > > >

Re: [Numpy-discussion] polyfit with fixed points

2013-03-04 Thread Aron Ahmadia
Interesting, that question would probably have gotten a different response on scicomp, it is a pity we are not attracting more questions there! I know there are two polyfit modules in numpy, one in numpy.polyfit, the other in numpy.polynomial, the functionality you are suggesting seems to fit in e

Re: [Numpy-discussion] MKL licenses for core scientific Python projects

2012-12-16 Thread Aron Ahmadia
All open source software and research projects with numpy in the stack, including PyClaw and petsc4py. A On Saturday, December 15, 2012, Ralf Gommers wrote: > > > > On Sat, Dec 15, 2012 at 2:21 PM, Aron Ahmadia > > > wrote: > >> Ralf, >> >> Does

Re: [Numpy-discussion] MKL licenses for core scientific Python projects

2012-12-15 Thread Aron Ahmadia
All open source software and research projects with numpy in the stack, including PyClaw and petsc4py. A On Saturday, December 15, 2012, Ralf Gommers wrote: > > > > On Sat, Dec 15, 2012 at 2:21 PM, Aron Ahmadia > > > wrote: > >> Ralf, >> >> Does

Re: [Numpy-discussion] MKL licenses for core scientific Python projects

2012-12-15 Thread Aron Ahmadia
Ralf, Does "performance testing" come under building/testing? If so, Aron Ahmadia OS X.8 Thanks, A On Sat, Dec 15, 2012 at 8:13 AM, Ralf Gommers wrote: > > > > On Sat, Dec 15, 2012 at 5:06 AM, Chris Barker - NOAA Federal < > chris.bar...@noaa.gov> wrote: >

Re: [Numpy-discussion] Z-ordering (Morton ordering) for numpy

2012-11-24 Thread Aron Ahmadia
Todd, I am optimistic and I think it would be a good idea to put this in. A couple previous studies [1] haven't found any useful speedups from in-core applications for Morton-order, and if you have results for real scientific applications using numpy this would not only be great, but the resultin

Re: [Numpy-discussion] error of install numpy on linux redhat.

2012-10-27 Thread Aron Ahmadia
Dear Jack, You are starting to ask questions that should be easily answered by 5 minutes of your own time using a search engine. Speaking for myself, I'm a volunteer here, and I want to help people learn and use Python, but I also expect that you will put some of your own effort into tracking dow

Re: [Numpy-discussion] error of install numpy on linux redhat.

2012-10-25 Thread Aron Ahmadia
I'm repeating myself, but did you try following the instructions here: http://stackoverflow.com/a/622810/122022? ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] error of install numpy on linux redhat.

2012-10-25 Thread Aron Ahmadia
Hi Jack, import os > Traceback (most recent call last): > File "", line 1, in > ImportError: No module named os You've got a broken local python install. Did you try following the instructions here: http://stackoverflow.com/a/622810/122022? It is probably easiest to ask your system admi

Re: [Numpy-discussion] NumPy to CPU+GPU compiler, looking for tests

2012-10-21 Thread Aron Ahmadia
Hi Rahul, Very cool! I'm looking forward to seeing some performance results! Anders Logg posted a computational challenge to G+ about a month ago, and we got entries in Octave, Fortran, Python, and Julia (all implementing the same solution from Jed Brown). The challenge is here: https://plus.g

Re: [Numpy-discussion] Long-standing issue with using numpy in embedded CPython

2012-09-25 Thread Aron Ahmadia
Can you expand a bit? Are you trying to disable threads at compile-time or at run-time? Which threaded functionality are you trying to disable? Are you using numpy as a computational library with multiple threads making calls into its functions? I think NPY_ALLOW_THREADS is for interacting with

Re: [Numpy-discussion] FWIW: "regressions" of dependees of numpy 1.7.0b1

2012-09-06 Thread Aron Ahmadia
Are you running the valgrind test with the Python suppression file: http://svn.python.org/projects/python/trunk/Misc/valgrind-python.supp ? Cheers, A On Wed, Sep 5, 2012 at 10:14 PM, Yaroslav Halchenko wrote: > and another, quite weird one -- initially it was crashing with the same > error on >

Re: [Numpy-discussion] blaze lib announcement

2012-08-29 Thread Aron Ahmadia
The Eigen3 project (http://eigen.tuxfamily.org/index.php?title=Main_Page) is more mature, but it's good to see a new contender in the field. A On Wed, Aug 29, 2012 at 1:59 PM, Neal Becker wrote: > This looks interesting: > > http://code.google.com/p/blaze-lib/ > > __

Re: [Numpy-discussion] minor threads-related issue in numpy-release 1.6.2 and numpy-dev

2012-08-14 Thread Aron Ahmadia
> Some parts of the numpy code release the GIL, so it needs to be reacquired > in order to call the Python C_API. A quick look shows a call to > PyObject_Print in the function and there may be other such calls. The right > fix is probably to use a couple of ifdefs if there is an easy way to > deter

[Numpy-discussion] minor threads-related issue in numpy-release 1.6.2 and numpy-dev

2012-08-14 Thread Aron Ahmadia
Hi all, Installing numpy 1.6.2 against a Python interpreter built with the --without threads currently fails due to missing references to PyGILState_Ensure and PyGILState_Release. The references appear to be coming from the following code in nditer.c.src: NPY_NO_EXPORT void NpyIter_DebugPrint(Np

Re: [Numpy-discussion] building numpy 1.6.2 on OSX 10.6 / Python2.7.3

2012-08-08 Thread Aron Ahmadia
FWIW, on my OS X.6.8 system, using a brew-installed python, and installing into /usr/local/lib, the symbols appear to be present: ❯ nm /usr/local/lib/python2.7/site-packages/numpy/core/multiarray.so | grep ceil U _ceil U _ceilf U _ceill 00

Re: [Numpy-discussion] building numpy 1.6.2 on OSX 10.6 / Python2.7.3

2012-08-08 Thread Aron Ahmadia
> > `python setup.py install --user` should install bumpy in a ~/.local > directory, you'll just have to update your PYTHONPATH > As of Python 2.6/3.x, ~/.local is searched after is added before the system site directories but after Python's search paths and PYTHONPATH. See PEP 370 for more detail

Re: [Numpy-discussion] Moving away from using accelerate framework on mac os x ?

2012-08-04 Thread Aron Ahmadia
Responding to both Dag and David, > OpenMP is pretty widespread in supercomputing, and so is OpenMP + > multiple threads using LAPACK at the same time. This does NOT sound like > any issue with multiprocessing to me. > These are built on dynamic thread-level (pthreads), not process-level (multip

Re: [Numpy-discussion] Moving away from using accelerate framework on mac os x ?

2012-08-04 Thread Aron Ahmadia
Hi David, Apple's response here is somewhat confusing, but I will add that on the supercomputing side of things we rarely fork, as this is not well-supported from the vendors or the hardware (it's hard enough to performantly spawn 500,000 processes statically, doing this dynamically becomes even m

Re: [Numpy-discussion] Pull Requests I'm planning to merge

2012-07-17 Thread Aron Ahmadia
I (or somebody else) needs to fix 325, Ralf pointed out that I was slightly too aggressive preferring gfortran on vendor operating systems. I think I should be able to fix this in the next couple of hours, I've been delinquent in getting back to it. A On Tue, Jul 17, 2012 at 3:56 PM, Travis Olip

Re: [Numpy-discussion] Numpy test failure - How to fix

2012-07-06 Thread Aron Ahmadia
Disclosure: I'm not a numpy developer You've certainly got a functional numpy installation if those are the only two tests failing. Those two tests are related to the distance between floating-point numbers using long double precision (128 bits). If you're not using long doubles, then you probab

Re: [Numpy-discussion] Numpy Installation Problem on Redhat Linux

2012-07-06 Thread Aron Ahmadia
I usually find these problems by searching for "error" in the output, in your case the complete problem is at the bottom of the log. The relocation errors you're seeing are happening because the build process is trying to link in Atlas libraries (located here: /home/hxc249/lib/atlas/lib/) that wer

Re: [Numpy-discussion] Preferring gfortran over g77 on OS X and other distributions?

2012-07-06 Thread Aron Ahmadia
> > Aron made this change for {mac, linux, posix, sun, irix, aix} in his PR. I > think it's a good idea for {mac, linux} and perhaps posix, and to leave the > other platforms alone. Does anyone else have an opinion on this? > Ralf, that sounds good to me. If you feel that we should leave sun, iri

Re: [Numpy-discussion] numpy and readline installation fails

2012-07-04 Thread Aron Ahmadia
How do the .dmg files work? Are these binary installers into the system Python? I would expect that these wouldn't work with a manually installed Python from python.org, but I have no experience with them. Tom, I have had very good luck with the brew Python 2.7 installer, which will give you a /

Re: [Numpy-discussion] Meta: help, devel and stackoverflow

2012-06-30 Thread Aron Ahmadia
I and Geoff are moderators on scicomp, I'm happy to invest the effort in getting the community started there. One way to use scicomp is like a blog/faq, that is, if you get a specific question a lot here on the list or elsewhere, you can ask and answer it yourself on scicomp. If others find the p

Re: [Numpy-discussion] Meta: help, devel and stackoverflow

2012-06-30 Thread Aron Ahmadia
As I mentioned before, numpy-related questions would be welcome on scicomp, and this would have the advantage of bringing in scientists and mathematicians from related fields who might be able to answer numerical questions that sit between mathematics, programming, and science that you might not ot

Re: [Numpy-discussion] Meta: help, devel and stackoverflow

2012-06-28 Thread Aron Ahmadia
Did you mean http://programmers.stackexchange.com? The meta sites on *. stackexchange.com are used (as one might guess) for meta discussions on the site. A On Thu, Jun 28, 2012 at 4:58 PM, Cera, Tim wrote: > Similar to http://scicomp.stackexchange.com there is > http://meta.programmers.stackex

Re: [Numpy-discussion] Meta: help, devel and stackoverflow

2012-06-28 Thread Aron Ahmadia
We try to support numpy questions on http://scicomp.stackexchange.com, which is a StackOverflow site dedicated towards technical computing issues that gets a fair amount of traffic from mathematicians and computational scientists. We could always use more questions and answerers :) A On Thu, Jun

Re: [Numpy-discussion] Preferring gfortran over g77 on OS X and other distributions?

2012-06-27 Thread Aron Ahmadia
r options alone. https://github.com/numpy/numpy/pull/325 A On Wed, Jun 27, 2012 at 10:46 PM, Ralf Gommers wrote: > > > On Mon, Jun 18, 2012 at 9:47 AM, Aron Ahmadia wrote: > >> f2py, by default, seems to prefer g77 (no longer maintained, deprecated, >> speedy, doesn&#

[Numpy-discussion] Preferring gfortran over g77 on OS X and other distributions?

2012-06-18 Thread Aron Ahmadia
f2py, by default, seems to prefer g77 (no longer maintained, deprecated, speedy, doesn't support Fortran 90 or Fortran 95) over gfortran (maintained, slower, Fortran 90 and Fortran 95 support). This causes problems when we try to compile Fortran 90 extensions using f2py on platforms where both g77