On Thu, Oct 15, 2015 at 10:51 PM, Adrian Croucher <[email protected]
> wrote:
> hi,
>
> I am wanting to output time-stepped simulation results to an HDF5 file,
> from a PETSc Vec created on a DMPlex using DMCreateGlobalVector().
>
> I tried doing it the way it's done in
> src/vec/vec/examples/tutorials/ex19.c, using PetscViewerHDF5Open(),
> PetscViewerHDF5PushGroup() and PetscViewerHDF5SetTimestep(), but the
> resulting HDF5 file only ever seems to have results for the last time step
> in it.
>
> It works OK for a Vec created simply using VecCreate(), as in the ex19
> example. Is there any reason it shouldn't work with DMPlex?
>
No, that should work. This is how we do it in PyLith. There must be some
setup problem. When you look at the HDF5 file. Does
it have a dimension for timestep?
Thanks,
Matt
> Cheers, Adrian
>
> --
> Dr Adrian Croucher
> Senior Research Fellow
> Department of Engineering Science
> University of Auckland, New Zealand
> email: [email protected]
> tel: +64 (0)9 923 84611
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
