<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Jun 8, 2014 at 3:02 PM, Blaise A Bourdin <span dir="ltr"><<a href="mailto:bourdin@lsu.edu" target="_blank">bourdin@lsu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I hate to be the party pooper, especially when I don’t have code to contribute, but I/O with DMComplex is still far from release ready. Is it reasonable that having functional I/O in a stable release is going to wait until 3.6?<br>
</blockquote><div><br></div><div>I agree with Barry that since this will not mean any interface changes, it can be done with patch releases.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<br>
I still can’t figure out how to deal with cell sets of different types (say quads and tri).<br></blockquote><div><br></div><div>I have not tried this. I had been focusing getting higher order (Q2) and either quads or tris to work. This</div>
<div>all seems to be working correctly now. I have been using it to visualize/restart a magma dynamics app</div><div>with 3 fields all with different discretizations.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

All cell and vertex sets seem to be contained in the hdf5 file, but not in a way that is usable by post processing tools (visit, paraview, ensight).<br></blockquote><div><br></div><div>For restarting, mixed meshes should work fine. They are stored in</div>
<div><br></div><div>  /topology/cells</div><div>                   /cones</div><div>                   /order</div><div>                   /orientation</div><div><br></div><div>and the field values are in</div><div><br></div>
<div>  /fields/<name></div><div><br></div><div>For visualization and most post-processing, there are separate arrays</div><div><br></div><div>  /viz/topology</div><div>        /geometry</div><div><br></div><div>which are exactly what I needed to make Paraview understand the xdmf. The fields</div>
<div>sampled down to cells and vertices are in</div><div><br></div><div>  /vertex_fields/<name></div><div>  /cell_fields/<name></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

The xdmf generation script bin/pythonscripts/petsc_gen_xdmf.py is quite fragile.<br></blockquote><div><br></div><div>I have not had it fail for me, but would be happy to look at the failure you are getting.</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">
I understand that it is probably too late to get this fixed and tested in 3.5. What is the best course of action?<br>
<br>
Blaise<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
>   We’ll be making a PETSc release in a few days so please get anything you need into next for testing and movement over to master.<br>
><br>
>  Barry<br>
><br>
> Tentative plan for release Wed June 11<br>
><br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Department of Mathematics and Center for Computation & Technology<br>
Louisiana State University, Baton Rouge, LA 70803, USA<br>
Tel. <a href="tel:%2B1%20%28225%29%20578%201612" value="+12255781612">+1 (225) 578 1612</a>, Fax  <a href="tel:%2B1%20%28225%29%20578%204276" value="+12255784276">+1 (225) 578 4276</a> <a href="http://www.math.lsu.edu/~bourdin" target="_blank">http://www.math.lsu.edu/~bourdin</a><br>

<br>
<br>
<br>
<br>
<br>
<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <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
</div></div>