On 13 Nov 2015, at 14:56, Buesing, Henrik <HBuesing@eonerc.rwth-aachen.de> wrote:
-----Ursprüngliche Nachricht----- Von: firedrake-bounces@imperial.ac.uk [mailto:firedrake- bounces@imperial.ac.uk] Im Auftrag von Lawrence Mitchell Gesendet: 13 November 2015 13:20 An: firedrake@imperial.ac.uk Betreff: Re: [firedrake] Output for time-stepping?
On 13 Nov 2015, at 12:14, Justin Chang <jychang48@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. Cheers, Lawrence