<div dir="ltr">You have to add code to get dm_view to work apparently.<div><br></div><div>I have done this before. This is some hocus pocus that Toby gave me. Is this really the way?</div><div><br></div><div><div>    PetscViewer       viewer = NULL;</div><div>    PetscBool         flg;</div><div>    PetscViewerFormat fmt;</div><div>    ierr = DMViewFromOptions(dm,NULL,"-dm_view");CHKERRQ(ierr);</div><div>    ierr = PetscOptionsGetViewer(ctx->wComm,NULL,"-x2_vec_view",&viewer,&fmt,&flg);CHKERRQ(ierr);</div><div>    if (flg) {</div><div>      ierr = PetscViewerPushFormat(viewer,fmt);CHKERRQ(ierr);</div><div>      ierr = VecView(dmpi->phi,viewer);CHKERRQ(ierr);</div><div>      ierr = VecView(dmpi->rho,viewer);CHKERRQ(ierr);</div><div>      ierr = PetscViewerPopFormat(viewer);CHKERRQ(ierr);</div></div><div><div>   } </div><div>   ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);</div></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 21, 2016 at 2:16 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  How about adding a least some documentation to bin/petsc_gen_xdmf.py<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
> On Dec 21, 2016, at 1:14 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
> On Wed, Dec 21, 2016 at 12:56 PM, Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br>
> I ran with:<br>
><br>
> -snes_monitor_residual hdf5:res.h5<br>
> -snes_monitor_solution hdf5:sol.h5<br>
><br>
> and got h5 files that look OK:<br>
><br>
> That will stash vectors in the HDF5 file, but you still need to put in the DM as well<br>
><br>
>   -dm_view hdf5:ex1.h5<br>
>  -snes_monitor_residual hdf5:ex1.h5::append<br>
>  -snes_monitor_solution hdf5:ex1.h5::append<br>
><br>
>     Matt<br>
><br>
> -rw-r--r--  1 markadams  staff   19760 Dec 21 13:46 res.h5<br>
> -rw-r--r--  1 markadams  staff   19760 Dec 21 13:46 sol.h5<br>
><br>
> but then got this. Am I doing this wrong?<br>
><br>
> 13:46 multi-species *+> ~/Codes/landaufem/Plex$ ~/Codes/petsc/bin/petsc_gen_<wbr>xdmf.py res.h5<br>
> Traceback (most recent call last):<br>
>   File "/Users/markadams/Codes/petsc/<wbr>bin/petsc_gen_xdmf.py", line 241, in <module><br>
>     generateXdmf(f)<br>
>   File "/Users/markadams/Codes/petsc/<wbr>bin/petsc_gen_xdmf.py", line 211, in generateXdmf<br>
>     geom      = h5['geometry']<br>
>   File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.<wbr>wrapper (/Users/travis/build/<wbr>MacPython/h5py-wheels/h5py/<wbr>h5py/_objects.c:2405)<br>
>   File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.<wbr>wrapper (/Users/travis/build/<wbr>MacPython/h5py-wheels/h5py/<wbr>h5py/_objects.c:2362)<br>
>   File "/Users/markadams/homebrew/<wbr>lib/python2.7/site-packages/<wbr>h5py/_hl/group.py", line 160, in __getitem__<br>
>     oid = h5o.open(<a href="http://self.id" rel="noreferrer" target="_blank">self.id</a>, self._e(name), lapl=self._lapl)<br>
>   File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.<wbr>wrapper (/Users/travis/build/<wbr>MacPython/h5py-wheels/h5py/<wbr>h5py/_objects.c:2405)<br>
>   File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.<wbr>wrapper (/Users/travis/build/<wbr>MacPython/h5py-wheels/h5py/<wbr>h5py/_objects.c:2362)<br>
>   File "h5py/h5o.pyx", line 190, in h5py.h5o.open (/Users/travis/build/<wbr>MacPython/h5py-wheels/h5py/<wbr>h5py/h5o.c:3456)<br>
> KeyError: "Unable to open object (Object 'geometry' doesn't exist)"<br>
><br>
><br>
><br>
> --<br>
> 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<br>
<br>
</div></div></blockquote></div><br></div>