Re: [Numpy-discussion] Zero row in SVD's unitary matrix on some Mac's

2011-04-25 Thread Rob Beezer
Thanks for the reply, Pauli. I suspected this might be the case, but was hoping maybe this was something that had been seen before. I've included your suggestions on the bug report for Sage. Rob On 04/25/2011 10:57 AM, Pauli Virtanen wrote: > The Numpy routine is a very thin wrapper of LAPAC

Re: [Numpy-discussion] Zero row in SVD's unitary matrix on some Mac's

2011-04-25 Thread Paul Anton Letnes
On 25. apr. 2011, at 19.57, Pauli Virtanen wrote: > On Mon, 25 Apr 2011 10:16:13 -0700, Rob Beezer wrote: > [clip] >> Many more details and complete transcripts are at: >> http://trac.sagemath.org/sage_trac/ticket/11248 >> >> Any thoughts or advice to help us understand this would be greatly >> a

Re: [Numpy-discussion] [Scikit-learn-general] random number generator, entropy and pickling

2011-04-25 Thread Gael Varoquaux
On Mon, Apr 25, 2011 at 01:23:12PM -0500, Robert Kern wrote: > > Yes, that's exactly why we want the different objects to able to recieve > > their own PRNG. > But seriously, they are running A+B, the combination of A and B. If > A+B changes to A+B', then the results may be different. That's to be

Re: [Numpy-discussion] [Scikit-learn-general] random number generator, entropy and pickling

2011-04-25 Thread Robert Kern
On Mon, Apr 25, 2011 at 13:15, Gael Varoquaux wrote: > On Mon, Apr 25, 2011 at 11:05:05AM -0700, T J wrote: >> If code A relies on code B (eg, some numpy function) and code B >> changes, then the stream of random numbers will no longer be the same. >>  The point here is that the user wrote code A

Re: [Numpy-discussion] random number generator, entropy and pickling

2011-04-25 Thread Gael Varoquaux
On Mon, Apr 25, 2011 at 11:05:05AM -0700, T J wrote: > If code A relies on code B (eg, some numpy function) and code B > changes, then the stream of random numbers will no longer be the same. > The point here is that the user wrote code A but depended on code B, > and even though code A was unchan

Re: [Numpy-discussion] random number generator, entropy and pickling

2011-04-25 Thread T J
On Mon, Apr 25, 2011 at 9:57 AM, Gael Varoquaux wrote: > > We thought that we could simply have a PRNG per object, as in: > >    def __init__(self, prng=None): >        if prng is None: >            prng = np.random.RandomState() >        self.prng = prng > > I don't like this option, because it m

Re: [Numpy-discussion] Zero row in SVD's unitary matrix on some Mac's

2011-04-25 Thread Pauli Virtanen
On Mon, 25 Apr 2011 10:16:13 -0700, Rob Beezer wrote: [clip] > Many more details and complete transcripts are at: > http://trac.sagemath.org/sage_trac/ticket/11248 > > Any thoughts or advice to help us understand this would be greatly > appreciated. The Numpy routine is a very thin wrapper of LAP

Re: [Numpy-discussion] [Scikit-learn-general] random number generator, entropy and pickling

2011-04-25 Thread Robert Kern
On Mon, Apr 25, 2011 at 11:57, Gael Varoquaux wrote: > Hi there, > > We are courrently having a discussion on the scikits learn mailing list > about which patterns to adopt for random number generation. One thing > that is absolutely clear is that making the stream of random numbers > reproducible

Re: [Numpy-discussion] [ANN]: OpenPiv. A python package for PIV image analysis

2011-04-25 Thread Alex Liberzon
Dear Fabrice, Let's take this discussion off-line, contact us directly through http://github.com/alexlib or to e-mail. It's a coincidence but I just finished wrapping the same library using SWIG :). Such a pity I didn't know about your work. I also believe that we talk about the same thing as

[Numpy-discussion] Zero row in SVD's unitary matrix on some Mac's

2011-04-25 Thread Rob Beezer
I am working to make many of NumPy's matrix decomposition routines available in Sage. As part of testing a new routine, we have found some odd behavior with the singular value decomposition. On certain Mac's, the numpy built in Sage will return the second of the unitary matrices with a row of

[Numpy-discussion] random number generator, entropy and pickling

2011-04-25 Thread Gael Varoquaux
Hi there, We are courrently having a discussion on the scikits learn mailing list about which patterns to adopt for random number generation. One thing that is absolutely clear is that making the stream of random numbers reproducible is critical. We have several objects that can serve as random va

Re: [Numpy-discussion] using numpy to save arrays in a loop

2011-04-25 Thread Isaiah Norton
There are various numpy and scipy input/output functions described here: http://www.scipy.org/Cookbook/InputOutput http://docs.scipy.org/doc/scipy/reference/io.html On Mon, Apr 25, 2011 at 10:32 AM, Zhaoru Zhang wrote: > Hi guys, > > I have a time loop and at each time step I get a 1-d array and

[Numpy-discussion] using numpy to save arrays in a loop

2011-04-25 Thread Zhaoru Zhang
Hi guys, I have a time loop and at each time step I get a 1-d array and want to save it in the data file 'data.dat '(or whatever suffix), so finally i got . I made this by using the python 'write' function, but I' thinking about if there is any numpy function that can realize this? Thanks. Zha

Re: [Numpy-discussion] Array as Variable using "from cdms2 import MV2 as MV"

2011-04-25 Thread Keith Goodman
On Mon, Apr 25, 2011 at 12:17 AM, wrote: > On Mon, Apr 25, 2011 at 2:50 AM, dileep kunjaai > wrote: >> Dear sir, >> >> I am have 2  mxn numpy array say "obs" & "fcst". I have to >> calculate sum of squre of (obs[i, j]-fcst[i, j]) using from cdms2 >> import MV2 as MV   in CDAT witho

Re: [Numpy-discussion] [ANN]: OpenPiv. A python package for PIV image analysis

2011-04-25 Thread Fabrice Silva
Le dimanche 24 avril 2011 à 16:53 +0200, Davide a écrit : > Hi all, > I am pleased to announce the re-birth of OpenPiv, a python package for > the analysis of PIV images. PIV is a well established optical > measurement technique for fluid flows. The package has the goal of > provide an implement

Re: [Numpy-discussion] Memmap with multiprocessing

2011-04-25 Thread Thiago Franco Moraes
Hi, Has anyone confirmed if this is a bug? Should I post this in the bug tracker? Thanks! On Tue, Apr 19, 2011 at 9:01 PM, Thiago Franco de Moraes wrote: > Hi all, > > I'm having a error using memmap objects shared among processes created > by the multprocessing module. This error only happen i

Re: [Numpy-discussion] cdms2.tvariable.TransientVariable to float variable

2011-04-25 Thread dileep kunjaai
SORRRY FRIENDS THAT I WRONGLY SENT On Mon, Apr 25, 2011 at 3:20 PM, dileep kunjaai wrote: > Dear sir, > How to change cdms2.tvariable.TransientVariable to float. Actually > this "cdms2.tvariable.TransientVariable" obtained from using from cdms2 > import MV2 > > T

[Numpy-discussion] cdms2.tvariable.TransientVariable to float variable

2011-04-25 Thread dileep kunjaai
Dear sir, How to change cdms2.tvariable.TransientVariable to float. Actually this "cdms2.tvariable.TransientVariable" obtained from using from cdms2 import MV2 Thanking you -- DILEEPKUMAR. R J R F, IIT DELHI ___ NumPy-Discussion mailing list Nu

Re: [Numpy-discussion] Array as Variable using "from cdms2 import MV2 as MV"

2011-04-25 Thread dileep kunjaai
Thank you sir ...thanks a lot Keep in touch On Mon, Apr 25, 2011 at 12:47 PM, wrote: > On Mon, Apr 25, 2011 at 2:50 AM, dileep kunjaai > wrote: > > Dear sir, > > > > I am have 2 mxn numpy array say "obs" & "fcst". I have to > >

Re: [Numpy-discussion] Array as Variable using "from cdms2 import MV2 as MV"

2011-04-25 Thread josef . pktd
On Mon, Apr 25, 2011 at 2:50 AM, dileep kunjaai wrote: > Dear sir, > > I am have 2  mxn numpy array say "obs" & "fcst". I have to > calculate sum of squre of (obs[i, j]-fcst[i, j]) using from cdms2 > import MV2 as MV   in CDAT without using "for"  loop. > > For example: > obs= > [0.6