Re: [Numpy-discussion] ANN: HDF5 for Python 1.1

2009-02-09 Thread Stephen Simmons
Hi Andrew, Do you have any plans to support LZO compression in h5py? I have lots of LZO-compressed datasets created with PyTables. There's a real barrier to using both h5py and PyTables if the fast decompressor options are just LZF on h5py and LZO on PyTables. Many thanks Stephen Andrew Colle

Re: [Numpy-discussion] preferred numpy build system

2009-02-09 Thread Ondrej Certik
On Mon, Feb 9, 2009 at 2:35 PM, Brian Granger wrote: >> CMake does handle this automatically. >> E.g. if include directories are changed (which you do by editing a >> CMakeLists.txt or the cmake cache), all files which are affected by the are >> rebuilt. If some library changes, everything linking

Re: [Numpy-discussion] ANN: HDF5 for Python 1.1

2009-02-09 Thread Andrew Collette
Thanks, Ondrej. For the record, h5py is designed to provide a "NumPy-like" interface to HDF5, along with a near-complete wrapping of the low-level HDF5 C API. It has none of the database-like features of PyTables. The FAQ entry has more info. Andrew Collette On Mon, Feb 9, 2009 at 1:06 PM, Ond

Re: [Numpy-discussion] Profiling with cProfile

2009-02-09 Thread Robert Kern
On Mon, Feb 9, 2009 at 15:35, Simon Palmer wrote: > > Hi, > I am trying to profile a bit of code I have written using cProfile. When I > run it I get the message: > > TypeError: unhashable type: 'numpy.ndarray' > > I am using runctx with some local variables which are ndarrays. > > I am guessing

Re: [Numpy-discussion] SVD errors

2009-02-09 Thread Robert Kern
On Mon, Feb 9, 2009 at 16:25, M Trumpis wrote: > I played around with a C translation of that test program, and found > that dgesvd (but not dgesdd) happens to converge and return all > non-negative singular values for both operators I was having trouble > with. I'm also looking at the Octave code

Re: [Numpy-discussion] preferred numpy build system

2009-02-09 Thread Brian Granger
> CMake does handle this automatically. > E.g. if include directories are changed (which you do by editing a > CMakeLists.txt or the cmake cache), all files which are affected by the are > rebuilt. If some library changes, everything linking to this library is > linked again. > If any of the files

Re: [Numpy-discussion] SVD errors

2009-02-09 Thread M Trumpis
I played around with a C translation of that test program, and found that dgesvd (but not dgesdd) happens to converge and return all non-negative singular values for both operators I was having trouble with. I'm also looking at the Octave code just now, and I think they're using dgesvd also. Any on

Re: [Numpy-discussion] preferred numpy build system

2009-02-09 Thread Ondrej Certik
On Sun, Feb 8, 2009 at 7:56 PM, David Cournapeau wrote: > Ondrej Certik wrote: >>> That's exactly what I don't like about cmake - it means you can't >>> produce accurate builds (you need to rerun cmake everytime you change >>> the configuration or dependencies, whereas this is automatic with >>>

[Numpy-discussion] Profiling with cProfile

2009-02-09 Thread Simon Palmer
Hi, I am trying to profile a bit of code I have written using cProfile. When I run it I get the message: TypeError: unhashable type: 'numpy.ndarray' I am using runctx with some local variables which are ndarrays. I am guessing that this is a lmitation of either cProfile or numpy or the mix of

Re: [Numpy-discussion] ANN: HDF5 for Python 1.1

2009-02-09 Thread Ondrej Certik
On Mon, Feb 9, 2009 at 12:15 PM, Andrew Collette wrote: > = > Announcing HDF5 for Python (h5py) 1.1 > = > > What is h5py? > - > > HDF5 for Python (h5py) is a general-purpose Python interface to the > Hierarchical D

[Numpy-discussion] ANN: HDF5 for Python 1.1

2009-02-09 Thread Andrew Collette
= Announcing HDF5 for Python (h5py) 1.1 = What is h5py? - HDF5 for Python (h5py) is a general-purpose Python interface to the Hierarchical Data Format library, version 5. HDF5 is a versatile, mature scientific so

Re: [Numpy-discussion] how to get the corresponding eigenvector for a specific eigen value?

2009-02-09 Thread Charles R Harris
On Mon, Feb 9, 2009 at 12:32 AM, Xiaoyu Chu wrote: > Hey all, > I am currently working on a large matrix, and I already have a > specific eigen value that I want to use in order to find out its > corresponding eigen vector. Is there an easy way to do so? > > I have tried with linalg.solv

Re: [Numpy-discussion] Comparison of arrays

2009-02-09 Thread Brent Pedersen
On Mon, Feb 9, 2009 at 6:02 AM, Neil wrote: > >> > I have two integer arrays of different shape, e.g. >> > >> > >>> a >> > >> > array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) >> > >> > >>> b >> > >> > array([ 3, 4, 5, 6, 7, 8, 9, 10]) >> > >> > How can I extract the values that belong to th

Re: [Numpy-discussion] Comparison of arrays

2009-02-09 Thread Francesc Alted
A Monday 09 February 2009, Neil escrigué: > > > I have two integer arrays of different shape, e.g. > > > > > > >>> a > > > > > > array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) > > > > > > >>> b > > > > > > array([ 3, 4, 5, 6, 7, 8, 9, 10]) > > > > > > How can I extract the values that belong

Re: [Numpy-discussion] Comparison of arrays

2009-02-09 Thread Neil
> > I have two integer arrays of different shape, e.g. > > > > >>> a > > > > array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) > > > > >>> b > > > > array([ 3, 4, 5, 6, 7, 8, 9, 10]) > > > > How can I extract the values that belong to the array a > > exclusively i.e. array([1,2]) ? > You cou

Re: [Numpy-discussion] linalg.norm missing an 'axis' kwarg?!

2009-02-09 Thread Stéfan van der Walt
2009/2/9 Hans Meine : >> Here's a basic implementation. docstring + tests not updated yet, also I >> wonder whether axis should be the first argument, but that could create >> compatibility problems. > > AFAICS, I never received an answer, but IMHO this should be integrated into > NumPy. Any obje

Re: [Numpy-discussion] linalg.norm missing an 'axis' kwarg?!

2009-02-09 Thread Hans Meine
On Thursday 20 November 2008 11:54:52 Hans Meine wrote: > On Thursday 20 November 2008 11:11:14 Hans Meine wrote: > > I have a 2D matrix comprising a sequence of vectors, and I want to > > compute the norm of each vector. np.linalg.norm seems to be the best > > bet, but it does not support axis.

Re: [Numpy-discussion] Comparison of arrays

2009-02-09 Thread Nils Wagner
On Mon, 9 Feb 2009 09:45:02 +0100 Francesc Alted wrote: > A Monday 09 February 2009, Nils Wagner escrigué: >> Hi all, >> >> I have two integer arrays of different shape, e.g. >> >> >>> a >> >> array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) >> >> >>> b >> >> array([ 3, 4, 5, 6, 7, 8, 9, 10

Re: [Numpy-discussion] Comparison of arrays

2009-02-09 Thread Francesc Alted
A Monday 09 February 2009, Nils Wagner escrigué: > Hi all, > > I have two integer arrays of different shape, e.g. > > >>> a > > array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) > > >>> b > > array([ 3, 4, 5, 6, 7, 8, 9, 10]) > > How can I extract the values that belong to the array a > exclusi

[Numpy-discussion] Comparison of arrays

2009-02-09 Thread Nils Wagner
Hi all, I have two integer arrays of different shape, e.g. >>> a array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]) >>> b array([ 3, 4, 5, 6, 7, 8, 9, 10]) How can I extract the values that belong to the array a exclusively i.e. array([1,2]) ? Nils __

Re: [Numpy-discussion] Selection of only a certain number of fields

2009-02-09 Thread Francesc Alted
A Sunday 08 February 2009, Neil escrigué: > > The first one (and most important IMO), is that newarr continues to > > be an structured array (BTW, when changed this name from the > > original record array?), and you can use all the features of these > > beasts with it. Other reason (albeit a bit s