<div dir="ltr"><div dir="ltr">On Wed, Mar 27, 2019 at 10:56 PM Jed Brown via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</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">Prior to this PR, which was merged for 3.10<br>
<br>
  <a href="https://bitbucket.org/petsc/petsc/pull-requests/1029/dmda-vtk-viewing-output-multiple-dof/diff" rel="noreferrer" target="_blank">https://bitbucket.org/petsc/petsc/pull-requests/1029/dmda-vtk-viewing-output-multiple-dof/diff</a><br>
  <a href="https://bitbucket.org/petsc/petsc/commits/ea2d7708fa6" rel="noreferrer" target="_blank">https://bitbucket.org/petsc/petsc/commits/ea2d7708fa6</a><br>
<br>
we could have files that look like the following, and thus were easy to<br>
navigate in Paraview and Visit.<br></blockquote><div><br></div><div>Tell me if I am wrong (since I do not run it this way). I think the point of Patrick's change is a good one, and</div><div>necessary to view things with DMStag. What Jed dislikes, rightly, is that some metadata has gone missing.</div><div>Does Xdmf allow you to name the components? The way I do this petsc_gen_xdmf.py is to write the whole</div><div>thing as Patrick does, but then create views into the data using the component names. It sound like this should</div><div>be doable here.</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">
<VTKFile type="StructuredGrid" version="0.1" byte_order="LittleEndian"><br>
  <StructuredGrid WholeExtent="0 12 0 24 0 24"><br>
    <Piece Extent="0 12 0 24 0 24"><br>
      <Points><br>
        <DataArray type="Float64" Name="Position" NumberOfComponents="3" format="appended" offset="0" /><br>
      </Points><br>
      <PointData Scalars="ScalarPointData"><br>
        <DataArray type="Float64" Name="Density" NumberOfComponents="1" format="appended" offset="195004" /><br>
        <DataArray type="Float64" Name="MomentumX" NumberOfComponents="1" format="appended" offset="260008" /><br>
        <DataArray type="Float64" Name="MomentumY" NumberOfComponents="1" format="appended" offset="325012" /><br>
        <DataArray type="Float64" Name="MomentumZ" NumberOfComponents="1" format="appended" offset="390016" /><br>
        <DataArray type="Float64" Name="Total Energy" NumberOfComponents="1" format="appended" offset="455020" /><br>
      </PointData><br>
    </Piece><br>
  </StructuredGrid><br>
  <AppendedData encoding="raw"><br>
<br>
Now the same viewer setup and VecView yields a file that is far more<br>
opaque and requires the user to remember what all the field numbers<br>
mean.<br>
<br>
<VTKFile type="StructuredGrid" version="0.1" byte_order="LittleEndian"><br>
  <StructuredGrid WholeExtent="0 12 0 24 0 24"><br>
    <Piece Extent="0 12 0 24 0 24"><br>
      <Points><br>
        <DataArray type="Float64" Name="Position" NumberOfComponents="3" format="appended" offset="0" /><br>
      </Points><br>
      <PointData Scalars="ScalarPointData"><br>
        <DataArray type="Float64" Name="Unnamed Vec data" NumberOfComponents="5" format="appended" offset="195004" /><br>
      </PointData><br>
    </Piece><br>
  </StructuredGrid><br>
  <AppendedData encoding="raw"><br>
<br>
<br>
Do others consider this a regression?  Is it feasible to fix for 3.11?<br>
</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>