[petsc-users] Data saving / viewer

[Gaetano Esposito]

Thanks, I was aware of such functions. I think I may have found where
my confusion is. If I have a regular Vec, HDF5 SetTimestep and
IncrementTimestep work as expected.
Instead the following does not work for me:

[...]
DMCreateGlobalVector(da, &solution);
[...]  // fill the vector solution and restore
HDF5Open
HDF5SetTimestep(H5_viewer, 0);

ierr = VecView(solution,H5_viewer); ERRCHKQ(ierr);
HDF5IncrementTimestep(H5_viewer);
ierr = VecView(solution,H5_viewer); ERRCHKQ(ierr);
[...]

I get the following external library error from HDF at the line of the
second VecView:

HDF5-DIAG: Error detected in HDF5 (1.8.6) MPI-process 0:
 #000: H5D.c line 170 in H5Dcreate2(): unable to create dataset
   major: Dataset
   minor: Unable to initialize object
 #001: H5Dint.c line 431 in H5D_create_named(): unable to create and
link to dataset
   major: Dataset
   minor: Unable to initialize object
 #002: H5L.c line 1640 in H5L_link_object(): unable to create new
link to object
   major: Links
   minor: Unable to initialize object
 #003: H5L.c line 1866 in H5L_create_real(): can't insert link
   major: Symbol table
   minor: Unable to insert object
 #004: H5Gtraverse.c line 929 in H5G_traverse(): internal path traversal failed
   major: Symbol table
   minor: Object not found
 #005: H5Gtraverse.c line 718 in H5G_traverse_real(): traversal operator failed
   major: Symbol table
   minor: Callback failed
 #006: H5L.c line 1675 in H5L_link_cb(): name already exists
   major: Symbol table
   minor: Object already exists