Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

--- On Mon, 6/8/09, Robert Kern wrote: > The OpenFormula standard is probably better: > > http://www.oasis-open.org/committees/documents.php?wg_abbrev=office-formula > > -- > Robert Kern OK, thanks Robert (as always); I'll go ahead and use this until/unless someone provide a printed referen

[Numpy-discussion] More on Summer NumPy Doc Marathon

Hi again, folks. I have a special request. Part of the vision for my job is that I'll focus my writing efforts on the docs no one else is gung-ho to work on. So, even if you're not quite ready to commit, if you're leaning toward volunteering to be a team lead for one (or more) categories, pl

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Tue, Jun 09, 2009 at 01:10:24PM +1200, Jochen Schroeder wrote: > On 09/06/09 00:16, Gael Varoquaux wrote: > > On Mon, Jun 08, 2009 at 05:14:27PM -0500, Gökhan SEVER wrote: > > >IPython's edit command works in a similar fashion, too. > > >edit test.py > > The cool thing is that you can

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On Tue, Jun 9, 2009 at 12:51 AM, wrote: > > --- On Mon, 6/8/09, Skipper Seabold wrote: > >> I forgot the last payment (which doesn't earn any >> interest), so one more 100. > > So in fact they're not in agreement? > >> pretty soon.  I don't have a more permanent reference >> for fv offhand, >> bu

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

--- On Mon, 6/8/09, Skipper Seabold wrote: > I forgot the last payment (which doesn't earn any > interest), so one more 100. So in fact they're not in agreement? > pretty soon.  I don't have a more permanent reference > for fv offhand, > but it should be in any corporate finance text etc.  > M

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On Tue, Jun 9, 2009 at 12:18 AM, Robert Kern wrote: > On Mon, Jun 8, 2009 at 22:54, David Goldsmith wrote: >> >> So would we regard a hard-copy of the users guide or reference manual for >> such a spreadsheet as sufficiently "permanent" to pass muster for use as a >> reference? > > The OpenFormul

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On Mon, Jun 8, 2009 at 22:54, David Goldsmith wrote: > > So would we regard a hard-copy of the users guide or reference manual for > such a spreadsheet as sufficiently "permanent" to pass muster for use as a > reference? The OpenFormula standard is probably better: http://www.oasis-open.org/com

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

So would we regard a hard-copy of the users guide or reference manual for such a spreadsheet as sufficiently "permanent" to pass muster for use as a reference? DG --- On Mon, 6/8/09, Alan G Isaac wrote: > From: Alan G Isaac > Subject: Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt,

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On 6/8/2009 11:18 PM David Goldsmith apparently wrote: > I formally "move" that numpy.financial (or at least that > subset of it consisting of functions which are commonly > subject to multiple definitions) be moved out of numpy. My recollection is that Travis O. added this with the explicit

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On Mon, Jun 8, 2009 at 9:01 PM, Skipper Seabold wrote: > Just quickly comparing > > In [3]: np.lib.financial.fv(.1,10,-100,-350) > Out[3]: 2501.5523211350032 > > With OO Calc > =fv(.1,10,-100,-350) > =2501.55 > > Both return the value of 350*1.1**10 + 100*1.1**9 + ... + 100*1.1 > which is what I w

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

--- On Mon, 6/8/09, Skipper Seabold wrote: > There was a recent discussion about numpy.financial in this > thread > . > > Skipper Thanks, Skipper. Having now read that thread (but not the arguments, provided elsewhere, for

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

--- On Mon, 6/8/09, Robert Kern wrote: > Goldsmith > wrote: > > > The advantage of my implementation is that I didn't have to > think too > hard about it. > -- > Robert Kern > Agreed. :-) DG ___ Numpy-discussion mailing list Numpy-discussi

[Numpy-discussion] From CorePy: New ExtBuffer object

Hi, Just a heads-up on something they're talking about over at CorePy. Regards Stéfan -- Forwarded message -- From: Andrew Friedley Date: 2009/6/8 Subject: [Corepy-devel] New ExtBuffer object To: CorePy Development I wrote a new buffer object today, called ExtBuffer, that can

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

Hi Juan 2009/6/8 Juanjo Gomez Navarro : > I'm new in numpy. Actually, I'm new in Python. In order to learn a bit, I > want to create a program to plot the Mandelbrot set. This program is quite > simple, and I have already programmed it. The problem is that I come from > fortran, so I use to think

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On 09/06/09 00:16, Gael Varoquaux wrote: > On Mon, Jun 08, 2009 at 05:14:27PM -0500, Gökhan SEVER wrote: > >IPython's edit command works in a similar fashion, too. > > >edit test.py > > The cool thing is that you can select text in the editor and execute in > EPDLab. On the other hand, I

Re: [Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

On Mon, Jun 8, 2009 at 5:44 PM, wrote: > > Hi, folks.  Unable to find a printed reference for the definition we use to > compute the functions in the Subject line of this email, I posted a couple > queries for help in this regard in the Discussion for fv > (http://docs.scipy.org/numpy/docs/numpy

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

On Mon, Jun 8, 2009 at 18:01, David Goldsmith wrote: > > --- On Mon, 6/8/09, Anne Archibald wrote: > >> > You can't, really. What you can do is just keep >> iterating with the >> > whole data set and ignore the parts that have already >> converged. Here >> > is an example: >> >> Well, yes and no.

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

--- On Mon, 6/8/09, Anne Archibald wrote: > > You can't, really. What you can do is just keep > iterating with the > > whole data set and ignore the parts that have already > converged. Here > > is an example: > > Well, yes and no. This is only worth doing if the number of > problem > points th

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

Thanks, Robert! DG --- On Mon, 6/8/09, Robert Kern wrote: > I haven't been able to figure out (though I haven't yet > really buckled down and tried real hard) how to broadcast a > conditionally-terminated iteration where the number of > iterations will vary among the array elements.  Hopefully

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

2009/6/8 Robert Kern : > On Mon, Jun 8, 2009 at 17:04, David Goldsmith wrote: >> >> I look forward to an instructive reply: the "Pythonic" way to do it would be >> to take advantage of the facts that Numpy is "pre-vectorized" and uses >> broadcasting, but so far I haven't been able to figure out

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

On Mon, Jun 8, 2009 at 17:04, David Goldsmith wrote: > > I look forward to an instructive reply: the "Pythonic" way to do it would be > to take advantage of the facts that Numpy is "pre-vectorized" and uses > broadcasting, but so far I haven't been able to figure out (though I haven't > yet real

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 08, 2009 at 05:14:27PM -0500, Gökhan SEVER wrote: >IPython's edit command works in a similar fashion, too. >edit test.py The cool thing is that you can select text in the editor and execute in EPDLab. On the other hand, I know that IPython has hooks to grow this in the code ba

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 4:47 PM, Christopher Barker wrote: > Gael Varoquaux wrote: > > Click in the menu: 'new file in remote browser', or something like this. > > If you have editra installed, it will launch it, with a special plugin > > allowing you to execute selected code in EPDLab. > > very co

Re: [Numpy-discussion] How to remove fortran-like loops with numpy?

I look forward to an instructive reply: the "Pythonic" way to do it would be to take advantage of the facts that Numpy is "pre-vectorized" and uses broadcasting, but so far I haven't been able to figure out (though I haven't yet really buckled down and tried real hard) how to broadcast a condi

[Numpy-discussion] Summer NumPy Doc Marathon (Reply-to: scipy-...@scipy.org)

Dear SciPy Community Members: Hi! My name is David Goldsmith. I've been hired for the summer by Joe Harrington to further progress on NumPy documentation and ultimately, pending funding, SciPy documentation. Joe and I are reviving last summer’s enthusiasm in the community for this miss

[Numpy-discussion] How to remove fortran-like loops with numpy?

Hi all, I'm new in numpy. Actually, I'm new in Python. In order to learn a bit, I want to create a program to plot the Mandelbrot set. This program is quite simple, and I have already programmed it. The problem is that I come from fortran, so I use to think in "for" loops. I know that it is not th

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 08, 2009 at 01:34:39PM -0700, Christopher Barker wrote: > Travis Oliphant wrote: > > You might take a look at EPDLab as well. Thanks to Gael Varoquaux, It > > integrates IPython into an Envisage application and has a crude > > name-space browser > I was wondering when you guys woul

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

Gael Varoquaux wrote: > Click in the menu: 'new file in remote browser', or something like this. > If you have editra installed, it will launch it, with a special plugin > allowing you to execute selected code in EPDLab. very cool, thanks! -Chris -- Christopher Barker, Ph.D. Oceanographer Eme

[Numpy-discussion] Definitions of pv, fv, nper, pmt, and rate

Hi, folks.  Unable to find a printed reference for the definition we use to compute the functions in the Subject line of this email, I posted a couple queries for help in this regard in the Discussion for fv (http://docs.scipy.org/numpy/docs/numpy.lib.financial.fv/#discussion-sec).  josef Pktd'

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 08, 2009 at 05:19:11PM -0400, David Warde-Farley wrote: > On 8-Jun-09, at 12:58 PM, Jonno wrote: > > Thanks Josef, > > I shouldn't have included Matplotlib since Pydee does work well with > > its plots. I had forgotten that. It really is just the Mayavi plots > > (or scenes I guess) t

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On 8-Jun-09, at 12:58 PM, Jonno wrote: > Thanks Josef, > > I shouldn't have included Matplotlib since Pydee does work well with > its plots. I had forgotten that. It really is just the Mayavi plots > (or scenes I guess) that don't play well. I don't know how exactly matplotlib integration issues

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 08, 2009 at 12:54:25PM -0500, Travis Oliphant wrote: >You might take a look at EPDLab as well. Thanks to Gael Varoquaux, It >integrates IPython into an Envisage application and has a crude name-space >browser And it integrates with editra to have an editor where you can s

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 15:34, Robert Kern wrote: > On Mon, Jun 8, 2009 at 15:34, Christopher Barker wrote: >> Travis Oliphant wrote: >>> You might take a look at EPDLab as well.   Thanks to Gael Varoquaux, It >>> integrates IPython into an Envisage application and has a crude >>> name-space browser

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 15:34, Christopher Barker wrote: > Travis Oliphant wrote: >> You might take a look at EPDLab as well.   Thanks to Gael Varoquaux, It >> integrates IPython into an Envisage application and has a crude >> name-space browser > > I was wondering when you guys would get around to

Re: [Numpy-discussion] matrix default to column vector?

On Mon, Jun 8, 2009 at 15:21, wrote: > 2009/6/8 Stéfan van der Walt : >> 2009/6/8 Robert Kern : Remember, the example is a **teaching** example. >>> >>> I know. Honestly, I would prefer that teachers skip over the normal >>> equations entirely and move directly to decomposition approaches. If

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

Travis Oliphant wrote: > You might take a look at EPDLab as well. Thanks to Gael Varoquaux, It > integrates IPython into an Envisage application and has a crude > name-space browser I was wondering when you guys would get around to making one of those. Nice start, the iPython shell is nice, t

Re: [Numpy-discussion] matrix default to column vector?

2009/6/8 Stéfan van der Walt : > 2009/6/8 Robert Kern : >>> Remember, the example is a **teaching** example. >> >> I know. Honestly, I would prefer that teachers skip over the normal >> equations entirely and move directly to decomposition approaches. If >> you are going to make them implement leas

Re: [Numpy-discussion] matrix default to column vector?

2009/6/8 Robert Kern : >> Remember, the example is a **teaching** example. > > I know. Honestly, I would prefer that teachers skip over the normal > equations entirely and move directly to decomposition approaches. If > you are going to make them implement least-squares from more basic > tools, I t

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

Gökhan SEVER wrote: > So far, IPython-Scite is the fastest that I can build my programs. > Experiment in Ipython and build pieces in Scite. I would like to know > what others use in this respect? Peppy (http://peppy.flipturn.org/) + iPython It would be nice to have those two integrated, though

[Numpy-discussion] The SciPy Doc Marathon continues

Let's Finish Documenting SciPy! Last year, we began the SciPy Documentation Marathon to write reference pages ("docstrings") for NumPy and SciPy. It was a huge job, bigger than we first imagined, with NumPy alone having over 2,000 functions. We created the doc wiki (now at doc

Re: [Numpy-discussion] matrix default to column vector?

On Mon, Jun 8, 2009 at 3:33 PM, Robert Kern wrote: > On Mon, Jun 8, 2009 at 14:10, Alan G Isaac wrote: Going back to Alan Isaac's example: 1)  beta = (X.T*X).I * X.T * Y 2)  beta = np.dot(np.dot(la.inv(np.dot(X.T,X)),X.T),Y) >> >> >> Robert Kern wrote: >>> 4) beta = la.lstsq(X, Y)[0]

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On 8-Jun-09, at 8:33 AM, Jason Rennie wrote:Note that EM can be very slow to converge:That's absolutely true, but EM for PCA can be a life saver in cases where diagonalizing (or even computing) the full covariance matrix is not a realistic option. Diagonalization can be a lot of wasted effort if al

[Numpy-discussion] New datetime dtypes

Hello, I'm working on the new datetime64 and timedelta64 dtypes (as proposed here: http://projects.scipy.org/numpy/browser/trunk/doc/neps/datetime-proposal3.rst). I'm looking through the C code in numpy core, and can't seem to find much in the way of dtypes. Pierre suggested looking through the m

Re: [Numpy-discussion] matrix default to column vector?

On Mon, Jun 8, 2009 at 14:10, Alan G Isaac wrote: >>> Going back to Alan Isaac's example: >>> 1)  beta = (X.T*X).I * X.T * Y >>> 2)  beta = np.dot(np.dot(la.inv(np.dot(X.T,X)),X.T),Y) > > > Robert Kern wrote: >> 4) beta = la.lstsq(X, Y)[0] >> >> I really hate that example. > > > Remember, the examp

[Numpy-discussion] Fwd: Re: is my numpy installation using custom blas/lapack?

Changing the site.cfg as you suggested did the trick! For what its worth, setup.py build no longer fails as before at compilation step (line 95), (I'm still puzzled whether this earlier 'failure' was caused by some error in my build process but I should probably let it go.) and numpy.show_c

Re: [Numpy-discussion] matrix default to column vector?

>> Going back to Alan Isaac's example: >> 1) beta = (X.T*X).I * X.T * Y >> 2) beta = np.dot(np.dot(la.inv(np.dot(X.T,X)),X.T),Y) Robert Kern wrote: > 4) beta = la.lstsq(X, Y)[0] > > I really hate that example. Remember, the example is a **teaching** example. I actually use NumPy in a Master'

Re: [Numpy-discussion] matrix default to column vector?

Olivier Verdier wrote: > Well, allowing dot(A,B,C) does not remove any other possibility does it? > I won't fight for this though. I personally don't care but I think that > it would remove the last argument for matrices against arrays, namely > the fact that A*B*C is easier to write than dot(do

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Jun 8, 2009, at 12:26 PM, Gökhan SEVER wrote: On Mon, Jun 8, 2009 at 12:11 PM, Jonno wrote: On Mon, Jun 8, 2009 at 11:35 AM, Gökhan SEVER wrote: > Hello, > > To me, IPython is the right way to follow. Try "whos" to see what's in your > namespace. > > You may want see this instructiona

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 12:11 PM, Jonno wrote: > On Mon, Jun 8, 2009 at 11:35 AM, Gökhan SEVER > wrote: > > Hello, > > > > To me, IPython is the right way to follow. Try "whos" to see what's in > your > > namespace. > > > > You may want see this instructional video (A Demonstration of the > 'IPyth

Re: [Numpy-discussion] matrix default to column vector?

2009/6/8 Christopher Barker > Olivier Verdier wrote: > > One > > should realize that allowing dot(A,B,C) is just *better* than the > > present situation where the user is forced into writing dot(dot(A,B),C) > > or dot(A,dot(B,C)). > > I'm lost now -- how is this better in any significant way?

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 11:35 AM, Gökhan SEVER wrote: > Hello, > > To me, IPython is the right way to follow. Try "whos" to see what's in your > namespace. > > You may want see this instructional video (A Demonstration of the 'IPython' > Interactive Shell) to learn more about IPython's functionality

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 11:39 AM, wrote: > > I'm using now pydee as my main shell to try out new scripts and I > don't have any problems with the plots. I'm creating plots the > standard way > from matplotlib import pyplot as plt > plt.plot(x,y) > > and I can close the poping up plot windows. > if

Re: [Numpy-discussion] matrix default to column vector?

Olivier Verdier wrote: > One > should realize that allowing dot(A,B,C) is just *better* than the > present situation where the user is forced into writing dot(dot(A,B),C) > or dot(A,dot(B,C)). I'm lost now -- how is this better in any significant way? Tom K. wrote: > But, > almost all experie

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

On Mon, Jun 8, 2009 at 12:29 PM, Jonno wrote: > Hi All, > > I'm new to python and tools like matplotlib and Mayavi so I may be > missing something basic. I've been looking for a fairly lightweight > editor/interactive shell combo that allows me to create plots and > figures from a shell and play w

Re: [Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

Hello, To me, IPython is the right way to follow. Try "whos" to see what's in your namespace. You may want see this instructional video (A Demonstration of the 'IPython' Interactive Shell) to learn more about IPython's functionality or

[Numpy-discussion] Interactive Shell/Editor/Workspace(variables)View/Plots

Hi All, I'm new to python and tools like matplotlib and Mayavi so I may be missing something basic. I've been looking for a fairly lightweight editor/interactive shell combo that allows me to create plots and figures from a shell and play with them and kill them gracefully. The Mayavi documentatio

Re: [Numpy-discussion] scipy 0.7.1rc2 released

Matthieu Brucher wrote: >> Concerning the other errors: did you compile with intel compilers or GNU >> ones ? >> > > Only Intel compilers. Maybe I should check the rc branch instead of the trunk? > I just wanted to confirm - I am actually rather surprised there are not more errors :) chee

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

Jason Rennie wrote: > > I hung-out in the machine learning community appx. 1999-2007 and > thought the Salakhutdinov work was extremely refreshing to see after > listening to no end of papers applying EM to whatever was the hot > topic at the time. :) Isn't it true for any general framework who en

Re: [Numpy-discussion] scipy 0.7.1rc2 released

2009/6/8 David Cournapeau : > Matthieu Brucher wrote: >> David, >> >> I've checked out the trunk, and the segmentation fault isn't there >> anymore (the trunk is labeled 0.8.0 though) >> > > Yes, the upcoming 0.7.1 release has its code in the 0.7.x svn branch. > But the fix for #946 is a backport o

Re: [Numpy-discussion] scipy 0.7.1rc2 released

Matthieu Brucher wrote: > David, > > I've checked out the trunk, and the segmentation fault isn't there > anymore (the trunk is labeled 0.8.0 though) > Yes, the upcoming 0.7.1 release has its code in the 0.7.x svn branch. But the fix for #946 is a backport of 0.8.0, so in theory, it should be f

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 8, 2009 at 8:55 AM, David Cournapeau < da...@ar.media.kyoto-u.ac.jp> wrote: > I think it depends on what you are doing - EM is used for 'real' work > too, after all :) Certainly, but EM is really just a mediocre gradient descent/hill climbing algorithm that is relatively easy to impl

Re: [Numpy-discussion] scipy 0.7.1rc2 released

David, I've checked out the trunk, and the segmentation fault isn't there anymore (the trunk is labeled 0.8.0 though) Here is the log from the remaining errors with the MKL: == ERROR: Failure: ImportError (/data/pau112/INNO/loca

Re: [Numpy-discussion] matrix default to column vector?

Is this lack of associativity really *always* such a huge issue? I can imagine many situations where it is not. One just want to compute A*B*C, without any particular knowing of whether A*(B*C) or (A*B)*C is best. If the user is allowed to blindly use A*B*C, I don't really see why he wouldn't be al

Re: [Numpy-discussion] scipy 0.7.1rc2 released

OK, I'm stuck with #946 with the MKL as well (finally managed to compile and use it with only the static library safe for libguide). I'm trying to download the trunk at the moment to check if the segmentation fault is still there. Matthieu 2009/6/8 Matthieu Brucher : > Good luck with fixing this

Re: [Numpy-discussion] setmember1d_nu

Robert Cimrman wrote: > Hi Neil, > > Neil Crighton wrote: >> Hi all, >> >> I posted this message couple of days ago, but gmane grouped it with an old >> thread and it hasn't shown up on the front page. So here it is again... >> >> I'd really like to see the setmember1d_nu function in ticket 1036

Re: [Numpy-discussion] extract elements of an array that are contained in another array?

Robert Cimrman wrote: > Hi Josef, > > thanks for the summary! I am responding below, later I will make an > enhancement ticket. Done, see http://projects.scipy.org/numpy/ticket/1133 r. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http:/

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 08, 2009 at 06:28:06AM -0700, Keith Goodman wrote: > On Mon, Jun 8, 2009 at 6:17 AM, Gael Varoquaux > wrote: > > On Mon, Jun 08, 2009 at 09:02:12AM -0400, josef.p...@gmail.com wrote: > >> whats the actual shape of the array/data you run your PCA on. > > 50 000 dimensions, 820 datapoin

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 8, 2009 at 6:17 AM, Gael Varoquaux wrote: > On Mon, Jun 08, 2009 at 09:02:12AM -0400, josef.p...@gmail.com wrote: >> whats the actual shape of the array/data you run your PCA on. > > 50 000 dimensions, 820 datapoints. Have you tried shuffling each time series, performing PCA, looking

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

2009/6/8 Gael Varoquaux : > On Mon, Jun 08, 2009 at 09:02:12AM -0400, josef.p...@gmail.com wrote: >> whats the actual shape of the array/data you run your PCA on. > > 50 000 dimensions, 820 datapoints. You definitely can't expect to find 50 meaningfull PCs. It's impossible to robustly get them wit

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 08, 2009 at 09:02:12AM -0400, josef.p...@gmail.com wrote: > whats the actual shape of the array/data you run your PCA on. 50 000 dimensions, 820 datapoints. > Number of time periods, size of cross section at point in time? I am not sure what the question means. The data is sampled at

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

Jason Rennie wrote: > Note that EM can be very slow to converge: > > http://www.cs.toronto.edu/~roweis/papers/emecgicml03.pdf > > > EM is great for churning-out papers, not so great for getting real > work done. I think it depends on wha

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 8, 2009 at 3:29 AM, Gael Varoquaux wrote: > On Mon, Jun 08, 2009 at 08:58:29AM +0200, Matthieu Brucher wrote: >> Given the number of PCs, I think you may just be measuring noise. >> As said in several manifold reduction publications (as the ones by >> Torbjorn Vik who published on robu

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 08, 2009 at 08:33:11AM -0400, Jason Rennie wrote: >EM is great for churning-out papers, not so great for getting real work >done.� That's just what I thought. >Btw, have you considered how much the Gaussianity assumption is >hurting you? I have. And the answer is: no

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

Note that EM can be very slow to converge: http://www.cs.toronto.edu/~roweis/papers/emecgicml03.pdf EM is great for churning-out papers, not so great for getting real work done. Conjugate gradient is a much better tool, at least in my (and Salakhutdinov's) experience. Btw, have you considered h

Re: [Numpy-discussion] scipy 0.7.1rc2 released

2009/6/8 David Cournapeau : > Matthieu Brucher wrote: >> Good luck with fixing this then :| >> >> I've tried to build scipy with the MKL and ATLAS, and I have in both >> cases a segmentation fault. With the MKL, it is the same as in a >> previous mail, and for ATLAS it is there: >> Regression test

Re: [Numpy-discussion] scipy 0.7.1rc2 released

Matthieu Brucher wrote: > Good luck with fixing this then :| > > I've tried to build scipy with the MKL and ATLAS, and I have in both > cases a segmentation fault. With the MKL, it is the same as in a > previous mail, and for ATLAS it is there: > Regression test for #946. ... Segmentation fault >

Re: [Numpy-discussion] scipy 0.7.1rc2 released

Good luck with fixing this then :| I've tried to build scipy with the MKL and ATLAS, and I have in both cases a segmentation fault. With the MKL, it is the same as in a previous mail, and for ATLAS it is there: Regression test for #946. ... Segmentation fault A bad ATLAS compilation? Matthieu >

Re: [Numpy-discussion] scipy 0.7.1rc2 released

On Mon, Jun 8, 2009 at 8:45 PM, Matthieu Brucher wrote: > 2009/6/8 David Cournapeau : >> Matthieu Brucher wrote: >>> I'm trying to compile it with ICC 10.1.018, and it fails :| >>> >>> icc: scipy/special/cephes/const.c >>> scipy/special/cephes/const.c(94): error: floating-point operation >>> result

Re: [Numpy-discussion] extract elements of an array that are contained in another array?

Hi Josef, thanks for the summary! I am responding below, later I will make an enhancement ticket. josef.p...@gmail.com wrote: > On Sat, Jun 6, 2009 at 4:42 AM, Neil Crighton wrote: >> Robert Cimrman ntc.zcu.cz> writes: >> >>> Anne Archibald wrote: >>> 1. add a keyword argument to intersec

Re: [Numpy-discussion] scipy 0.7.1rc2 released

2009/6/8 David Cournapeau : > Matthieu Brucher wrote: >> I'm trying to compile it with ICC 10.1.018, and it fails :| >> >> icc: scipy/special/cephes/const.c >> scipy/special/cephes/const.c(94): error: floating-point operation >> result is out of range >>   double INFINITY = 1.0/0.0;  /* 99e999; */

Re: [Numpy-discussion] scipy 0.7.1rc2 released

Matthieu Brucher wrote: > I'm trying to compile it with ICC 10.1.018, and it fails :| > > icc: scipy/special/cephes/const.c > scipy/special/cephes/const.c(94): error: floating-point operation > result is out of range > double INFINITY = 1.0/0.0; /* 99e999; */ >^ > > scipy

Re: [Numpy-discussion] scipy 0.7.1rc2 released

I'm trying to compile it with ICC 10.1.018, and it fails :| icc: scipy/special/cephes/const.c scipy/special/cephes/const.c(94): error: floating-point operation result is out of range double INFINITY = 1.0/0.0; /* 99e999; */ ^ scipy/special/cephes/const.c(99): error: floa

Re: [Numpy-discussion] scipy 0.7.1rc2 released

2009/6/8 Matthieu Brucher : > I'm trying to compile it with ICC 10.1.018, and it fails :| > > icc: scipy/special/cephes/const.c > scipy/special/cephes/const.c(94): error: floating-point operation > result is out of range >  double INFINITY = 1.0/0.0;  /* 99e999; */ >                       ^ > > sci

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

2009/6/8 David Warde-Farley : > > On 8-Jun-09, at 1:17 AM, David Cournapeau wrote: > >> I would not be surprised if David had this paper in mind :) >> >> http://www.cs.toronto.edu/~roweis/papers/empca.pdf > > Right you are :) > > There is a slight trick to it, though, in that it won't produce an >

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

2009/6/8 Gael Varoquaux : > On Mon, Jun 08, 2009 at 08:58:29AM +0200, Matthieu Brucher wrote: >> Given the number of PCs, I think you may just be measuring noise. >> As said in several manifold reduction publications (as the ones by >> Torbjorn Vik who published on robust PCA for medical imaging),

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On Mon, Jun 08, 2009 at 08:58:29AM +0200, Matthieu Brucher wrote: > Given the number of PCs, I think you may just be measuring noise. > As said in several manifold reduction publications (as the ones by > Torbjorn Vik who published on robust PCA for medical imaging), you > cannot expect to have mor

Re: [Numpy-discussion] performance matrix multiplication vs. matlab

On 8-Jun-09, at 1:17 AM, David Cournapeau wrote: > I would not be surprised if David had this paper in mind :) > > http://www.cs.toronto.edu/~roweis/papers/empca.pdf Right you are :) There is a slight trick to it, though, in that it won't produce an orthogonal basis on its own, just something

Re: [Numpy-discussion] Functions for indexing into certain parts of an array (2d)

On Sun, Jun 7, 2009 at 6:19 AM, Bruce Southey wrote: > On Sun, Jun 7, 2009 at 3:37 AM, Fernando Perez wrote: >> One more question.  For these *_indices() functions, would you want an >> interface that accepts *either* >> >> diag_indices(size,ndim) > > As I indicated above, this is unacceptable for