<div dir="ltr">Good evening,<div><br></div><div>My group and I are in the process of developing an application using PETSc, and we are currently planning and figuring out the best objects to do this. I currently have some questions about visualization and I/O.</div><div><br></div><div>Some of the most important features that we want in the process of developing this code is the ability for easy checkpoint/restarts and visualization. As we are working with a time dependent PDE, I have stumbled upon the "ts" and "tstrajectory" objects that seem to have some very helpful features. However, I am having some trouble extracting information out of the the "tstrajectory" object in an HDF5 format. I initially thought something like this would work, where tj was obtained using TSGetTrajectory(ts,&tj) after a successful solve of the ts object:</div><div><br></div><div>PetscViewer viewer;<br>PetscViewerHDF5Open(PETSC_COMM_WORLD,"output.h5",FILE_MODE_WRITE,&viewer);CHKERRQ(ierr);<br>TSTrajectoryView(tj,viewer);CHKERRQ(ierr); <br>PetscViewerDestroy(&viewer);CHKERRQ(ierr);<br></div><div><div><br></div><div>However, the output file comes out empty. Do I have to extract the solutions one by one? How would I do that? Does the HDF5 viewer only work with vector objects or can it be generalized to other objects? I am trying to maintain the highest degree of abstraction possible (i.e. multiple 2d ode solutions in one hdf5 file).</div><div><br></div><div>This seems favorable to me because it would be easy to use the HDF5 file for postprocessing/visualization and restarts, but let me know if you think this is a bad idea and it would be better to just use the HDF5 for output/postprocess and use native PETSc binary files for restarts. </div><div><br></div><div>Additionally, as I move towards large visualizations, I would like to use Paraview. I saw that the best way of doing this is to use the <span style="color:rgb(0,0,0);white-space:pre-wrap">petsc_gen_xdmf.py script on an HDF5 solution file. Is this still the best way for PETSc visualization with Paraview?</span></div><div><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></div><div><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></div><div>Thank you,</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><font face="arial, sans-serif">Alfredo Duarte</font><div><font face="arial, sans-serif">Graduate Research Assistant</font></div><div><font face="arial, sans-serif">The University of Texas at Austin</font></div></div></div></div></div></div></div>