<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Jan 17, 2018 at 7:35 AM, Marek Pecha <span dir="ltr"><<a href="mailto:marek.pecha@vsb.cz" target="_blank">marek.pecha@vsb.cz</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I create very simple program for writing scalar to HDF5 file (it is just test-program, because I learn how HDF5 works), see below.<br>
<br>
#include <petscviewerhdf5.h><br>
<br>
int main(int argc, char *argv[]) {<br>
PetscErrorCode ierr;<br>
PetscScalar n = 5;<br>
PetscViewer viewer;<br>
<br>
PetscInitialize(&argc,&argv,(<wbr>char*)0,NULL);<br>
<br>
PetscViewerHDF5Open(PETSC_<wbr>COMM_WORLD, "test.hdf5", FILE_MODE_WRITE, &viewer);<br>
PetscViewerHDF5PushGroup(<wbr>viewer, "/");<br>
ierr = PetscViewerHDF5WriteAttribute(<wbr>viewer, "/", "n", PETSC_SCALAR, &n);<br>
CHKERRQ(ierr);<br>
PetscViewerDestroy(&viewer);<br>
<br>
PetscFinalize();<br>
<br>
return (0);<br>
}<br>
<br>
However, I got following error:<br>
<br>
HDF5-DIAG: Error detected in HDF5 (1.8.12) MPI-process 0:<br>
#000: H5D.c line 341 in H5Dopen2(): not a dataset<br></blockquote><div><br></div><div>Sorry it took me so long to respond. I was on winter break.</div><div><br></div><div>The problem is that HDF5 attributes can only be applied to datasets. '/' is not a dataset. For example,</div><div>if you wrote a Vec into the file, you could apply an attribute to it, like 'axial vector'.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
major: Dataset<br>
minor: Inappropriate type<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>
[0]PETSC ERROR: Error in external library<br>
[0]PETSC ERROR: Error in HDF5 call H5Dopen2() Status -1<br>
[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.<br>
[0]PETSC ERROR: Petsc Release Version 3.7.5, unknown<br>
[0]PETSC ERROR: ./aijmathdf5 on a arch-llvm-opt named sweety by dev Wed Jan 17 13:33:51 2018<br>
[0]PETSC ERROR: Configure options --download-blacs --download-boost --download-fblaslapack --download-metis --download-mpich --download-mumps --download-parmetis --download-scalapack --download-superlu --download-superlu_dist --download-hdf5 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 FOPTFLAGS=-O3 PETSC_ARCH=arch-llvm-opt<br>
[0]PETSC ERROR: #1 PetscViewerHDF5WriteAttribute(<wbr>) line 653 in /opt/petsc/petsc-maint/src/<wbr>sys/classes/viewer/impls/hdf5/<wbr>hdf5v.c<br>
[0]PETSC ERROR: #2 main() line 13 in /Users/dev/Desktop/sparse_<wbr>hdf5/main.cpp<br>
[0]PETSC ERROR: No PETSc Option Table entries<br>
[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov-------<wbr>---<br>
application called MPI_Abort(MPI_COMM_WORLD, 76) - process 0<br>
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=76<br>
:<br>
system msg for write_line failure : Bad file descriptor<br>
<br>
I don’t know, what I am doing wrong.<br>
<br>
Thank you<br>
<span class="HOEnZb"><font color="#888888">Marek</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>