On 14 July 2011 18:33, Ethan Coon <ecoon at lanl.gov> wrote: > On Thu, 2011-07-14 at 17:47 -0300, Lisandro Dalcin wrote: > >> >> 1) readMatDense() is wrong. The I,J,V arrays are the CSR structure, >> not the COO (coordinate) format, Then you cannot simply: >> > > Ok, thanks, I wasn't sure of which format you were using,
Well, PETSc's AIJ matrices are CSR in memory (but not in disk, binary files store row counts instead of row pointers) >> >> 2) It would be nice to provie a 'scipy' Mat format returning a >> 'scipy.sparse.csr_matrix' instance. A possible implementation would >> be: >> > > Agreed, I considered doing that as the default for sparse, but wasn't > sure what others would want. ?Providing both makes sense. > Note however that this would require a SciPy installation. PS: Sorry for not being more helpful and contributing actual code, but I'm really busy. -- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169