pushed the patch to petsc-dev satish
On Tue, 19 Jul 2011, Ethan Coon wrote: > Taking this to dev. > > Attached is a patch for PetscBinaryRead, which now has been tested on > Mat, IS, and Vec. Supports returning dense numpy.matrix, sparse CSR, > and scipy.sparse matrices (all have been tested). > > Also cc'ed Ata directly since I'm not sure he's on dev. Ata -- could > you please test the attached PetscBinaryRead.py for me on your system to > make sure I've correctly handled your version of numpy? Note that the > return format has changed from what you have now. Thanks. > > PetscBinaryWrite.py when I get to it... > > Thanks, > > Ethan > > -------- Forwarded Message -------- > From: Lisandro Dalcin <dalcinl at gmail.com> > Reply-to: PETSc users list <petsc-users at mcs.anl.gov> > To: PETSc users list <petsc-users at mcs.anl.gov> > Subject: Re: [petsc-users] PETSc recommended visualization packages > Date: Thu, 14 Jul 2011 19:34:42 -0300 > > 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. > > > >
