Re: [Numpy-discussion] Github migration?

2010-09-01 Thread Charles R Harris
Hi Jason, On Tue, Aug 31, 2010 at 2:56 PM, Jason McCampbell wrote: > Hi Chuck (and anyone else interested), > > I updated the refactoring page on the NumPy developer wiki (seems to be > down or I'd paste in the link). It certainly isn't complete, but there are > a lot more details about the dat

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread David Cournapeau
On Thu, Sep 2, 2010 at 10:03 AM, Charles R Harris wrote: > > > On Wed, Sep 1, 2010 at 3:30 PM, Michael Gilbert > wrote: >> >> On Wed, 1 Sep 2010 21:15:22 + (UTC), Pauli Virtanen wrote: >> > Wed, 01 Sep 2010 16:26:59 -0400, Michael Gilbert wrote: >> > > I've been using numpy's float96 class la

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Charles R Harris
On Wed, Sep 1, 2010 at 3:30 PM, Michael Gilbert wrote: > On Wed, 1 Sep 2010 21:15:22 + (UTC), Pauli Virtanen wrote: > > Wed, 01 Sep 2010 16:26:59 -0400, Michael Gilbert wrote: > > > I've been using numpy's float96 class lately, and I've run into some > > > strange precision errors. > > [clip]

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Colin Macdonald
On 09/01/10 22:30, Michael Gilbert wrote: > Interesting. float96( '0.0001' ) also seems to evaluate to the first > result. I assume that it also does a float64( '0.0001' ) conversion > first. I understand that you can't change how python passes in floats, > but wouldn't it be better to exactly han

Re: [Numpy-discussion] kron produces F-contiguous?

2010-09-01 Thread Lisandro Dalcin
On 1 September 2010 19:24, Neal Becker wrote: > It seems if I call kron with 2 C-contiguous arrays, it returns an F- > contiguous array.  Any reason for this (it's not what I wanted)? > Try numpy.linalg.inv ... -- Lisandro Dalcin --- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa F

[Numpy-discussion] kron produces F-contiguous?

2010-09-01 Thread Neal Becker
It seems if I call kron with 2 C-contiguous arrays, it returns an F- contiguous array. Any reason for this (it's not what I wanted)? ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] Array slices and number of dimensions

2010-09-01 Thread Anne Archibald
On 1 September 2010 17:54, Thomas Robitaille wrote: > Hi, > > I'm trying to extract sub-sections of a multidimensional array while keeping > the number of dimensions the same. If I just select a specific element along > a given direction, then the number of dimensions goes down by one: > im

Re: [Numpy-discussion] Array slices and number of dimensions

2010-09-01 Thread Pauli Virtanen
Wed, 01 Sep 2010 17:54:26 -0400, Thomas Robitaille wrote: > >>> a[[0],:,:].shape > (1, 10, 10) [clip] > so I can get around the issue, but just wanted to check whether the > issue with a[[0],:,[5]] is a bug? No. The syntax does not mean what you think it means, see http://docs.scipy.org/doc/numpy

Re: [Numpy-discussion] Array slices and number of dimensions

2010-09-01 Thread Warren Weckesser
Thomas Robitaille wrote: > Hi, > > I'm trying to extract sub-sections of a multidimensional array while keeping > the number of dimensions the same. If I just select a specific element along > a given direction, then the number of dimensions goes down by one: > > > > In fact, I can get what

[Numpy-discussion] Array slices and number of dimensions

2010-09-01 Thread Thomas Robitaille
Hi, I'm trying to extract sub-sections of a multidimensional array while keeping the number of dimensions the same. If I just select a specific element along a given direction, then the number of dimensions goes down by one: >>> import numpy as np >>> a = np.zeros((10,10,10)) >>> a.shape (10, 1

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Michael Gilbert
On Wed, 1 Sep 2010 21:15:22 + (UTC), Pauli Virtanen wrote: > Wed, 01 Sep 2010 16:26:59 -0400, Michael Gilbert wrote: > > I've been using numpy's float96 class lately, and I've run into some > > strange precision errors. > [clip] > > >>> x = numpy.array( [0.01] , numpy.float96 ) > [clip] > > I

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Pauli Virtanen
Wed, 01 Sep 2010 16:26:59 -0400, Michael Gilbert wrote: > I've been using numpy's float96 class lately, and I've run into some > strange precision errors. [clip] > >>> x = numpy.array( [0.01] , numpy.float96 ) [clip] > I would expect the float96 calculation to also produce 0.0 exactly as > found

Re: [Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Charles R Harris
On Wed, Sep 1, 2010 at 2:26 PM, Michael Gilbert wrote: > Hi, > > I've been using numpy's float96 class lately, and I've run into some > strange precision errors. See example below: > > >>> import numpy > >>> numpy.version.version > '1.5.0' > >>> sys.version > '3.1.2 (release31-maint, Jul 8

[Numpy-discussion] Unexpected float96 precision loss

2010-09-01 Thread Michael Gilbert
Hi, I've been using numpy's float96 class lately, and I've run into some strange precision errors. See example below: >>> import numpy >>> numpy.version.version '1.5.0' >>> sys.version '3.1.2 (release31-maint, Jul 8 2010, 01:16:48) \n[GCC 4.4.4]' >>> x = numpy.array( [0.01] , numpy.

[Numpy-discussion] Developer Job Openings at Space Telescope Science Institute

2010-09-01 Thread Perry Greenfield
We are advertising for two different positions at the Space Telescope Science Institute (located on the Johns Hopkins University Campus in Baltimore, Md). STScI is seeking Senior Systems Software Engineers to develop applications to calibrate and analyze data from the Hubble and the James W

[Numpy-discussion] Parsing PyArrays

2010-09-01 Thread Babak Ahmadi
Hi, I'm having a little Problem that might be trivial. I'm trying to use a nupmy array within a c++ method and having some trouble. The PyArray is not parsed correctly and PyArg_ParseTuple returns 0. Appreciate any help. #include "Python.h" #include void fromMatrix(PyObject *args){ _import_a

Re: [Numpy-discussion] Github migration?

2010-09-01 Thread Charles R Harris
On Tue, Aug 31, 2010 at 2:56 PM, Jason McCampbell wrote: > Hi Chuck (and anyone else interested), > > I updated the refactoring page on the NumPy developer wiki (seems to be > down or I'd paste in the link). It certainly isn't complete, but there are > a lot more details about the data structure

Re: [Numpy-discussion] Where is the dev version of numpydoc?

2010-09-01 Thread Ralf Gommers
On Wed, Sep 1, 2010 at 11:31 AM, John Salvatier wrote: > Hello, > > I would like to update my numpydoc so it works with sphinx 1.0, but I am > not sure where the dev version is; can someone point me in the right > direction? > > In numpy trunk, under doc/sphinxext/. That works *only* with sphinx 1

Re: [Numpy-discussion] inversion of large matrices

2010-09-01 Thread Sebastian Walter
is it really the covariance matrix you want to invert? Or do you want to compute something like x^T C^{-1} x, where x is an array of size N and C an array of size (N,N)? It would also be interesting to know how the covariance matrix gets computed and what its condition number is, at least approxim