[petsc-users] Creating HDF5 groups

[Tim Gallagher]

Hi,

I'm a little confused about why this section of code doesn't work, so hopefully 
somebody can help me. I'm able to run the vec/vec/examples/tutorials/ex19.c 
test that creates an HDF5 file, and h5dump verifies that it is making the 
groups and putting vectors in them. However, my code does not make the group 
(so it obviously doesn't put the vector in it!). The function is:

#undef __FUNCT__
#define __FUNCT__ "Body::SaveGrid"
PetscErrorCode Body::SaveGrid(const char filename[])
{
  PetscViewer hdf5File;
  Vec coordVec;

  // Grab the coordinate vector                                                 
                                                                                
                      
  m_lastError = DMDAGetCoordinates(m_nodes,&coordVec);
  CHKERRXX(m_lastError);

  m_lastError = PetscObjectSetName((PetscObject) coordVec, "Coordinates");
  CHKERRXX(m_lastError);

  // Open the file                                                              
                                                                                
                      
  m_lastError = PetscViewerHDF5Open(m_commGroup, filename, FILE_MODE_WRITE, 
&hdf5File);
  CHKERRXX(m_lastError);

  // Set the group to dump the grid in                                          
                                                                                
                      
  m_lastError = PetscViewerHDF5PushGroup(hdf5File, "/grid");
  CHKERRXX(m_lastError);

  // Dump the coordinates                                                       
                                                                                
                      
  m_lastError = VecView(coordVec, hdf5File);
  CHKERRXX(m_lastError);

  m_lastError = PetscViewerHDF5PopGroup(hdf5File);
  CHKERRXX(m_lastError);

  // Clean up                                                                   
                                                                                
                      
  m_lastError = PetscViewerDestroy(&hdf5File);
  CHKERRXX(m_lastError);
}

The only thing I can really see different about this and ex19 is that I don't 
first write any dataset to the root group, but that's not required to be valid 
HDF5. If I call PetscViewerHDF5GetGroup after the PushGroup call, it returns 
/grid. 

Any ideas?

Thanks,

Tim