Re: [Numpy-discussion] Tuning sparse stuff in NumPy

2007-03-27 Thread David Koch
On 3/27/07, Robert Cimrman [EMAIL PROTECTED] wrote: ok. now which version of scipy (scipy.__version__) do you use (you may have posted it, but I missed it)? Not so long ago, there was an effort by Nathan Bell and others reimplementing sparsetools + scipy.sparse to get better usability and

Re: [Numpy-discussion] Tuning sparse stuff in NumPy

2007-03-27 Thread Robert Cimrman
David Koch wrote: On 3/27/07, Robert Cimrman [EMAIL PROTECTED] wrote: ok. now which version of scipy (scipy.__version__) do you use (you may have posted it, but I missed it)? Not so long ago, there was an effort by Nathan Bell and others reimplementing sparsetools + scipy.sparse to get

Re: [Numpy-discussion] Tuning sparse stuff in NumPy

2007-03-26 Thread David Koch
On 3/26/07, Robert Cimrman [EMAIL PROTECTED] wrote: Could you be more specific on which type of the sparse matrix storage did you use? Hi Robert, I used csc_matrix. /David ___ Numpy-discussion mailing list Numpy-discussion@scipy.org

Re: [Numpy-discussion] Tuning sparse stuff in NumPy

2007-03-26 Thread Robert Cimrman
David Koch wrote: On 3/26/07, Robert Cimrman [EMAIL PROTECTED] wrote: Could you be more specific on which type of the sparse matrix storage did you use? Hi Robert, I used csc_matrix. OK, good. Would you mind measuring csc * csr, csc * csc, csr * csc and csr * csr? I am curious how

Re: [Numpy-discussion] Tuning sparse stuff in NumPy

2007-03-26 Thread David Koch
Ok, I did and the results are: csc * csc: 372.601957083 csc * csc: 3.90811300278 csr * csc: 15.3202679157 csr * csr: 3.84498214722 Mhm, quite insightful. Note, that in an operation X.transpose() * X, where X is csc_matrix, then X.tranpose() is automatically cast to csr_matrix. A re-cast to csc