[petsc-users] petsc4py, plotting DA and writing to file

Barry Smith bsmith at mcs.anl.gov
Mon Sep 27 12:07:14 CDT 2010


On Sep 27, 2010, at 11:54 AM, Jed Brown wrote:

> On Mon, Sep 27, 2010 at 18:47, Lisandro Dalcin <dalcinl at gmail.com> wrote:
>> viewer = PETSc.Viewer.DRAW(globalvec.comm)
>> globalvec.view(viewer) # or viewer(globalvec)
> 
> What do you think about having a top-level matplotlib viewer?  Ideally
> it would have an option to drop you into an interactive python session
> after the initial view.
> 
> I know this brings up the multiple dispatch problem, since the viewer
> would naturally be distributed with petsc4p, but then you don't get to
> modify the case statement in VecView_XX.

   This could be done by putting into the C code the usual "case" handling for the new viewer and then have
it dispatch back to the python code. 

   Barry

> 
> Jed



More information about the petsc-users mailing list