Re: [Numpy-discussion] confusion about eigenvector

2008-03-01 Thread [EMAIL PROTECTED]
> I dont know if this made anything any clearer. However, a simple > example may be clearer: thanks Arnar for the kind response,now things are a lot clearer...will try out in code .. D ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Damian Eads
Lisandro Dalcin wrote: > On 3/1/08, Charles R Harris <[EMAIL PROTECTED]> wrote: >> So they differ in the least significant bit. Not surprising, I expect the >> Fortran compiler might well perform operations in different order, >> accumulate in different places, etc. It might also accumulate in high

Re: [Numpy-discussion] how to pronounce numpy?

2008-03-01 Thread Damian Eads
Robert Kern wrote: > On Sat, Mar 1, 2008 at 6:45 PM, Lisandro Dalcin <[EMAIL PROTECTED]> wrote: >> Sorry for the stupid question, but my English knowledge just covers >> reading and writting (the last, not so good) >> >> At the very begining, http://scipy.org/ says >> >> SciPy (pronounced "Sigh

Re: [Numpy-discussion] A little help please?

2008-03-01 Thread Travis E. Oliphant
Travis E. Oliphant wrote: > Neal Becker wrote: > >> Travis E. Oliphant wrote: >> >> >> >> >> The code for this is a bit hard to understand. It does appear that it only >> searches for a conversion on the 2nd argument. I don't think that's >> desirable behavior. >> >> What I'm wondering is,

Re: [Numpy-discussion] Rename record array fields (with object arrays)

2008-03-01 Thread Travis E. Oliphant
Travis E. Oliphant wrote: > Sameer DCosta wrote: > >> Hi, >> >> I'm having trouble renaming record array fields if they contain object >> arrays in them. I followed the solutions posted by Robert Kern and >> Stefan van der Walt (Thanks again) but it doesn't look like this >> method works in all

Re: [Numpy-discussion] UFUNC_CHECK_STATUS cpu hog

2008-03-01 Thread Thomas Grill
Am 02.03.2008 um 04:24 schrieb Travis E. Oliphant: if(__builtin_expect(fpstatus,0)) \ Why the use of __builtin_expect here instead of fpstatus == 0? Oops, nevertheless it should rather be something like if(__builtin_expect(fpstatus == 0,1)) or if(__builtin_expect(fpstatus,0) == 0) sorry f

Re: [Numpy-discussion] UFUNC_CHECK_STATUS cpu hog

2008-03-01 Thread Thomas Grill
Am 02.03.2008 um 04:24 schrieb Travis E. Oliphant: Thomas Grill wrote: Hi all, i did some profiling on OS X/Intel 10.5 (numpy 1.0.4) and was surprised to find calls to the system function feclearexcept to be by far the biggest cpu hog, taking away about 30% of the cpu in my case. Would it be p

Re: [Numpy-discussion] UFUNC_CHECK_STATUS cpu hog

2008-03-01 Thread Travis E. Oliphant
Thomas Grill wrote: > Hi all, > i did some profiling on OS X/Intel 10.5 (numpy 1.0.4) and was > surprised to find calls to the system function feclearexcept to be by > far the biggest cpu hog, taking away about 30% of the cpu in my case. > Would it be possible to change UFUNC_CHECK_STATUS in ufu

[Numpy-discussion] UFUNC_CHECK_STATUS cpu hog

2008-03-01 Thread Thomas Grill
Hi all, i did some profiling on OS X/Intel 10.5 (numpy 1.0.4) and was surprised to find calls to the system function feclearexcept to be by far the biggest cpu hog, taking away about 30% of the cpu in my case. Would it be possible to change UFUNC_CHECK_STATUS in ufuncobject.h in a way that

Re: [Numpy-discussion] how to pronounce numpy?

2008-03-01 Thread Robert Kern
On Sat, Mar 1, 2008 at 6:45 PM, Lisandro Dalcin <[EMAIL PROTECTED]> wrote: > Sorry for the stupid question, but my English knowledge just covers > reading and writting (the last, not so good) > > At the very begining, http://scipy.org/ says > > SciPy (pronounced "Sigh Pie") ... > > Then, for th

[Numpy-discussion] how to pronounce numpy?

2008-03-01 Thread Lisandro Dalcin
Sorry for the stupid question, but my English knowledge just covers reading and writting (the last, not so good) At the very begining, http://scipy.org/ says SciPy (pronounced "Sigh Pie") ... Then, for the other guy, this assertion NumPy (pronounced "Num Pie", "Num" as in "Number") ... whould

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Lisandro Dalcin
On 3/1/08, Charles R Harris <[EMAIL PROTECTED]> wrote: > So they differ in the least significant bit. Not surprising, I expect the > Fortran compiler might well perform operations in different order, > accumulate in different places, etc. It might also accumulate in higher > precision registers or

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Charles R Harris
2008/3/1 Lisandro Dalcin <[EMAIL PROTECTED]>: > Dear Charles, > > As I said, I have no time to code the pure Python+numpy nonlinear and > linear loops, and the matrix-free stuff to mimic the PETSc > implementation. However, I post the F90 code and the numpy code, and a > small script for testing w

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Lisandro Dalcin
Dear Charles, As I said, I have no time to code the pure Python+numpy nonlinear and linear loops, and the matrix-free stuff to mimic the PETSc implementation. However, I post the F90 code and the numpy code, and a small script for testing with random input. When I have some spare time, I'll try to

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Lisandro Dalcin
On 3/1/08, Pauli Virtanen <[EMAIL PROTECTED]> wrote: > A silly question: did you check directly that the pure-numpy code and > the F90 code give the same results for the Jacobian-vector product > J(z0) z for some randomly chosen vectors z0, z? No, I did not do that. However, I've checked the out

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Charles R Harris
On Sat, Mar 1, 2008 at 12:43 PM, Lisandro Dalcin <[EMAIL PROTECTED]> wrote: > Dear all, > > I want to comment some extrange stuff I'm experiencing with numpy. > Please, let me know if this is expected and known. > > I'm trying to solve a model nonlinear PDE, 2D Bratu problem (-Lapacian > u - alpha

Re: [Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Pauli Virtanen
Hi, la, 2008-03-01 kello 16:43 -0300, Lisandro Dalcin kirjoitti: > I want to comment some extrange stuff I'm experiencing with numpy. > Please, let me know if this is expected and known. > > I'm trying to solve a model nonlinear PDE, 2D Bratu problem (-Lapacian > u - alpha * exp(u), homogeneus bo

[Numpy-discussion] numpy and roundoff(?)

2008-03-01 Thread Lisandro Dalcin
Dear all, I want to comment some extrange stuff I'm experiencing with numpy. Please, let me know if this is expected and known. I'm trying to solve a model nonlinear PDE, 2D Bratu problem (-Lapacian u - alpha * exp(u), homogeneus bondary conditions), using the simple finite differences with a 5-p

Re: [Numpy-discussion] svd() and eigh()

2008-03-01 Thread Arnar Flatberg
On Sat, Mar 1, 2008 at 2:43 PM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > hi > i have a set of images of faces which i make into a 2d array using > numpy.ndarray > each row represents a face image > faces= > [[ 173. 87. ... 88. 165.] > [ 158. 103. .. 73. 143.] > [ 180. 87

Re: [Numpy-discussion] confusion about eigenvector

2008-03-01 Thread Arnar Flatberg
On Sat, Mar 1, 2008 at 8:27 AM, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > > This example assumes that facearray is an ndarray.(like you described > > in original post ;-) ) It looks like you are using a matrix. > > hi Arnar > thanks .. > a few doubts however > > 1.when i use say 10 ima

[Numpy-discussion] svd() and eigh()

2008-03-01 Thread [EMAIL PROTECTED]
hi i have a set of images of faces which i make into a 2d array using numpy.ndarray each row represents a face image faces= [[ 173. 87. ... 88. 165.] [ 158. 103. .. 73. 143.] [ 180. 87. .. 55. 143.] [ 155. 117. .. 93. 155.]] from which i can get the mean image => avgface=a

Re: [Numpy-discussion] failure building numpy using icc

2008-03-01 Thread Robert Kern
On Thu, Feb 28, 2008 at 1:21 PM, Glen W. Mabey <[EMAIL PROTECTED]> wrote: > Hello, > > I'm using svn numpy and get the following error upon executing > > /usr/local/bin/python2.5 setup.py config --noisy > --cc=/opt/intel/cce/10.0.025/bin/icc --compiler=intel --fcompiler=intel > build_clib buil