[firedrake] Output for time-stepping?

[Buesing, Henrik]

> -----Ursprüngliche Nachricht-----
> Von: firedrake-bounces at imperial.ac.uk [mailto:firedrake-
> bounces at imperial.ac.uk] Im Auftrag von Lawrence Mitchell
> Gesendet: 13 November 2015 16:20
> An: firedrake at imperial.ac.uk
> Betreff: Re: [firedrake] Output for time-stepping?
> 
> 
> > On 13 Nov 2015, at 14:56, Buesing, Henrik <HBuesing at eonerc.rwth-
> aachen.de> wrote:
> >
> >> -----Ursprüngliche Nachricht-----
> >> Von: firedrake-bounces at imperial.ac.uk [mailto:firedrake-
> >> bounces at imperial.ac.uk] Im Auftrag von Lawrence Mitchell
> >> Gesendet: 13 November 2015 13:20
> >> An: firedrake at imperial.ac.uk
> >> Betreff: Re: [firedrake] Output for time-stepping?
> >>
> >>
> >>> On 13 Nov 2015, at 12:14, Justin Chang <jychang48 at gmail.com> wrote:
> >>>
> >>> Lawrence,
> >>>
> >>> So the first option is what I have been doing. Yes only one pvd is
> >> created. But what ends up also happening is that I get hundreds of
> >> .vtu files. Are those files necessary for pvd to work? Sorry I am not
> >> in front of my computer so I can't test any of this out yet.
> >>
> >>
> >> Aha, I see.  Yes, unfortunately that's how it is :(.  Eventually we
> >> intend to have nice HDF5 output for viz, but it's not there yet.
> >>
> >> Cheers,
> >>
> >> Lawrence
> >
> > [Buesing, Henrik] Since you can also access PETSc directly, shouldn't
> you be able to use PetscViewerVTKOpen or PetscViewerHDF5Open for the
> output?
> 
> I think this requires that the Vec you get knows about the topology
> through a DM of some kind.  Unfortunately for a large class of the
> firedrake meshes (particularly the extruded meshes), although we have a
> DM describing the topology it's not "correct" for visualisation.
> 
[Buesing, Henrik] Ah, I see. That's true, you would need a DM or DMPLEX. I thought one might extract this from the KSP or SNES.