[petsc-dev] [Fwd: Re: [petsc-users] PETSc recommended visualization packages]

[Satish Balay]

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.
> 
> 
> 
>