Re: [Numpy-discussion] Scipy 2011 Convore thread now open

2011-07-12 Thread Long Duong
Does anybody know if there are there videos of the conference this year? Best regards, Long Duong UC Irvine Biomedical Engineering l...@studioart.org On Tue, Jul 12, 2011 at 1:00 PM, Peter Wang wrote: > Hi folks, > > I have gone ahead and created a Convore group for the SciPy 2011 > confere

[Numpy-discussion] named ndarray axes

2011-07-12 Thread Craig Yoshioka
I brought up a while ago about how it would be nice if numpy arrays could have their axes 'labeled'.= I got an implementation that works pretty well for me and in the process learned quite a few things, and was hoping to foster some more discussion on this topic, as I think I have found a si

[Numpy-discussion] Scipy 2011 Convore thread now open

2011-07-12 Thread Peter Wang
Hi folks, I have gone ahead and created a Convore group for the SciPy 2011 conference: https://convore.com/scipy-2011/ I have already created threads for each of the tutorial topics, and once the conference is underway, we'll create threads for each talk, so that audience can interact and post q

Re: [Numpy-discussion] object scalars

2011-07-12 Thread Johann Hibschman
Olivier Delalleau writes: > 2011/7/12 Johann Hibschman > > Is there any way to wrap a sequence (in particular a python list) as a > numpy object scalar, without it being promoted to an object array? > I found a workaround but it's a bit ugly: > def some_call(x): >   rval = numpy.array(N

Re: [Numpy-discussion] New arrays in 1.6 not always C-contiguous

2011-07-12 Thread Frédéric Bastien
Hi, On Tue, Jul 12, 2011 at 12:48 PM, Mark Wiebe wrote: [...] > > This only added to the C-API, pre-existing API remained the same for API/ABI > compatibility reasons. C code already had to deal with the possibility of > differing memory layouts, for example if someone passes in carr.T, something

Re: [Numpy-discussion] New arrays in 1.6 not always C-contiguous

2011-07-12 Thread Mark Wiebe
2011/7/12 Frédéric Bastien > Hi, > > We depend highly on numpy, but don't have the time to follow all the > mailing lists regularly of all our tools. Having some information on > the release note about this would have been useful to many people I > think. > You're absolutely right, not including

Re: [Numpy-discussion] New arrays in 1.6 not always C-contiguous

2011-07-12 Thread Frédéric Bastien
Hi, We depend highly on numpy, but don't have the time to follow all the mailing lists regularly of all our tools. Having some information on the release note about this would have been useful to many people I think. Also, did this affect the C-API? Do the default value of newlly created ndarray

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Daniel Wheeler
On Tue, Jul 12, 2011 at 11:19 AM, Sturla Molden wrote: > Den 11.07.2011 23:01, skrev Daniel Wheeler: > To make the loop over N matrices fast, there is nothing that beats a > loop in C or Fortran (or Cython) if you have a 3D array. And that brings > us to the second issue, which is that it would be

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Gael Varoquaux
On Tue, Jul 12, 2011 at 12:16:29PM -0400, greg whittier wrote: > Gael, your code addresses inverses, but I take it something similar for > eigenvalues of a matrix bigger than 5x5 doesn't exists since a > closed-form solution doesn't exist for finding polynomials roots for > order > 5? I guess so :

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread greg whittier
On Tue, Jul 12, 2011 at 11:30 AM, Gael Varoquaux wrote: > On Mon, Jul 11, 2011 at 05:01:07PM -0400, Daniel Wheeler wrote: >> Hi, I am trying to find the eigenvalues and eigenvectors as well as >> the inverse for a large number of small matrices. The matrix size >> (MxM) will typically range from 2

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Daniel Wheeler
On Tue, Jul 12, 2011 at 10:52 AM, Dag Sverre Seljebotn wrote: > On 07/12/2011 04:10 PM, Daniel Wheeler wrote: >> On Tue, Jul 12, 2011 at 3:51 AM, Dag Sverre Seljebotn >> Thanks for the heads up. Looks like an option. Presumably, it would >> still have to use "map" even with more direct access to B

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Gael Varoquaux
On Mon, Jul 11, 2011 at 05:01:07PM -0400, Daniel Wheeler wrote: > Hi, I am trying to find the eigenvalues and eigenvectors as well as > the inverse for a large number of small matrices. The matrix size > (MxM) will typically range from 2x2 to 8x8 at most. If you really care about speed, for matri

Re: [Numpy-discussion] object scalars

2011-07-12 Thread Olivier Delalleau
I found a workaround but it's a bit ugly: def some_call(x): rval = numpy.array(None, dtype='object') rval.fill(x) return rval -=- Olivier 2011/7/12 Johann Hibschman > Is there any way to wrap a sequence (in particular a python list) as a > numpy object scalar, without it being promoted to

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Sturla Molden
Den 11.07.2011 23:01, skrev Daniel Wheeler: > The above uses "map" to fake a vector solution, but this is heinously > slow. Are there any better ways to do this without resorting to cython > or weave (would it even be faster (or possible) to use "np.linalg.eig" > and "np.linalg.inv" within cython)?

[Numpy-discussion] object scalars

2011-07-12 Thread Johann Hibschman
Is there any way to wrap a sequence (in particular a python list) as a numpy object scalar, without it being promoted to an object array? In particular, np.object_([1, 2]).shape == (2,) np.array([1,2], dtype='O').shape == (2,) while I want some_call([1,2]).shape = () Thanks, Johann

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Dag Sverre Seljebotn
On 07/12/2011 04:10 PM, Daniel Wheeler wrote: > On Tue, Jul 12, 2011 at 3:51 AM, Dag Sverre Seljebotn > wrote: >> On 07/11/2011 11:01 PM, Daniel Wheeler wrote: >>> Hi, I am trying to find the eigenvalues and eigenvectors as well as >>> the inverse for a large number of small matrices. The matrix

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Daniel Wheeler
On Tue, Jul 12, 2011 at 3:51 AM, Dag Sverre Seljebotn wrote: > On 07/11/2011 11:01 PM, Daniel Wheeler wrote: >> Hi, I am trying to find the eigenvalues and eigenvectors as well as >> the inverse for a large number of small matrices. The matrix size > If you want to go the Cython route, here's a s

Re: [Numpy-discussion] Are .M and .H removed in NumPy 1.6?

2011-07-12 Thread Sturla Molden
Den 12.07.2011 15:55, skrev Charles R Harris: On Tue, Jul 12, 2011 at 7:36 AM, Sturla Molden > wrote: After upgrading EPD, I just discovered that my ndarrays no longer have .M and .H attributes. Were they deprectated, or is my NumPy not working correctly?

Re: [Numpy-discussion] Are .M and .H removed in NumPy 1.6?

2011-07-12 Thread Charles R Harris
On Tue, Jul 12, 2011 at 7:36 AM, Sturla Molden wrote: > After upgrading EPD, I just discovered that my ndarrays no longer have > .M and .H attributes. > > Were they deprectated, or is my NumPy not working correctly? > > I thought they were long gone: http://mail.scipy.org/pipermail/numpy-discussi

[Numpy-discussion] Are .M and .H removed in NumPy 1.6?

2011-07-12 Thread Sturla Molden
After upgrading EPD, I just discovered that my ndarrays no longer have .M and .H attributes. Were they deprectated, or is my NumPy not working correctly? Sturla ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman

Re: [Numpy-discussion] inverting and calculating eigenvalues for many small matrices

2011-07-12 Thread Dag Sverre Seljebotn
On 07/11/2011 11:01 PM, Daniel Wheeler wrote: > Hi, I am trying to find the eigenvalues and eigenvectors as well as > the inverse for a large number of small matrices. The matrix size > (MxM) will typically range from 2x2 to 8x8 at most. The number of > matrices (N) can be from 100 up to a million