[Numpy-discussion] Optimization suggestion sought

2010-12-26 Thread Enzo Michelangeli
For a pivoted algorithm, I have to perform an operation that in fully vectorized form can be expressed as: pivot = tableau[locat,:]/tableau[locat,cand] tableau -= tableau[:,cand:cand+1]*pivot tableau[locat,:] = pivot tableau is a rather large bidimensional array, and I'd like to avoid

Re: [Numpy-discussion] Optimization suggestion sought

2010-12-27 Thread Enzo Michelangeli
4b433fbf94fa137eaa5ee69a06987eda\Release\docume~1\admin\locals~1\temp\admin\python26_compiled\sc_c9159c98b571d3181d8848337bf1e50a1.o" failed with exit status 1 \-- So, for the time being, no speed breakthrough... Enzo - Original Mes

Re: [Numpy-discussion] Optimization suggestion sought

2011-01-01 Thread Enzo Michelangeli
- Original Message - From: "Robert Bradshaw" Sent: Wednesday, December 29, 2010 4:47 PM [...] >> Regarding Justin's suggestion, before trying Cython (which, according to >> http://wiki.cython.org/tutorials/numpy , seems to require a bit of work >> to >> handle numpy arrays properly) > >

[Numpy-discussion] Arrays with aliased elements?

2011-01-01 Thread Enzo Michelangeli
Is there any way, not involving compilation of C code, to define ndarrays where some rows or columns share the same data buffers? For example, something built by a hypothetical variant of the np.repeat() function, such that, if a = array([2,3]), calling: b = np.aliasedrepeat(x, [1, 2], axis=0

Re: [Numpy-discussion] Arrays with aliased elements?

2011-01-01 Thread Enzo Michelangeli
s the intended savings in memory and CPU cycles... Enzo - Original Message - From: "Robert Kern" To: "Discussion of Numerical Python" Sent: Sunday, January 02, 2011 10:08 AM Subject: Re: [Numpy-discussion] Arrays with aliased elements? > On Sat, Jan 1, 2011 at