<div dir="ltr"><div dir="ltr">On Mon, Jul 13, 2020 at 2:51 AM Lisandro Dalcin <<a href="mailto:dalcinl@gmail.com">dalcinl@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Fri, 10 Jul 2020 at 15:37, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><br></div><div class="gmail_quote"><div>That is true. Do they also get rid of the single mesh and single timstep requirements? HDF5+XDMF makes it much easier since</div><div>we can put multiple meshes and timesteps in one file.</div></div></div></blockquote><div><br></div><div>Well, I use ParaView's *.pvd files for that, and I dump each timestep to its own *.vtu file in a folder to pack the files. But you still have a point, the format is indeed restricting.</div><div>Why do you consider it so important to put multiple meshes and timesteps in one file? It is just that you hate to have so many files scattered around? Or something deeper?</div><div>I do hate the fact that the VTK formats (either legacy or XML) do not allow you to dump a single mesh to be reused for multiple timestep.</div></div></div></blockquote><div><br></div><div>For now its complexity of moving simulation data around and scripting for it. However, now I have at least two meshes in my problem, and</div><div>I anticipate having several more. I believe this will be the long term trend.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div>I would probably move out of VTK files in favor of something else if I had a way to encode VTK's (the library, not the file format) high-order Lagrange elements.<br></div><div>Actually, I'm toying with dumping files with PETSc's raw binary I/O with MPI, and writing a proper ParaView plugin in Python to read the data.</div></div></div></blockquote><div><br></div><div>I have again discussed higher order with the Firedrake people. They are using the Paraview mechanism, but it is so fragile and baroque that I</div><div>refuse. Currently, the cleanest way is still to refine the mesh and project to P_1. There have been so many attempts at high order in VTK, I can</div><div>only conclude that Kitware does not give a crap about it.</div><div><br></div><div>  Thanks,</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div>-- <br></div></div><div dir="ltr"><div dir="ltr"><div>Lisandro Dalcin<br>============<br>Research Scientist<br>Extreme Computing Research Center (ECRC)<br>King Abdullah University of Science and Technology (KAUST)<br><a href="http://ecrc.kaust.edu.sa/" target="_blank">http://ecrc.kaust.edu.sa/</a><br></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>