Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-27 Thread Matthew Knepley
On Mon, Oct 26, 2015 at 10:03 PM, Adrian Croucher wrote: > hi > > When I read the HDF5 file back in again using PetscViewerHDF5Open() and > VecLoad(), is there any easy way to navigate to the last time step in the > results? > > Or, equivalently, to find how many timesteps there are in there, so

Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-26 Thread Adrian Croucher
hi When I read the HDF5 file back in again using PetscViewerHDF5Open() and VecLoad(), is there any easy way to navigate to the last time step in the results? Or, equivalently, to find how many timesteps there are in there, so I could then use PetscViewerHDF5SetTimestep() before doing VecLoad

Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-18 Thread Adrian Croucher
hi, On 17/10/15 03:01, Matthew Knepley wrote: Now I remember. I did not want the output to depend on the viewer. Does your example work if you replace PetscViewerHDF5SetTimestep() with DMSetOutputSequenceNumber()? Aha, yes, that appears to do the trick. It gives me the time values in the o

Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-16 Thread Matthew Knepley
On Thu, Oct 15, 2015 at 11:14 PM, Adrian Croucher wrote: > hi > > On 16/10/15 16:54, Matthew Knepley wrote: > > > 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? > > > No, it doesn't

Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-15 Thread Adrian Croucher
hi On 16/10/15 16:54, Matthew Knepley wrote: 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? No, it doesn't. Is there any extra setup I might need to do to get it to work? I've a

Re: [petsc-users] HDF5 vec output with DMPlex

2015-10-15 Thread Matthew Knepley
On Thu, Oct 15, 2015 at 10:51 PM, Adrian Croucher 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

[petsc-users] HDF5 vec output with DMPlex

2015-10-15 Thread Adrian Croucher
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 PetscViewerH