[petsc-dev] getting data out of a PetscFE simulation

[Matthew Knepley]

On Tue, Jan 28, 2014 at 6:04 PM, Geoffrey Irving <irving at naml.us> wrote:

> What's the cleanest way to save data from a petsc finite element
> simulation including boundary conditions?  First, is it correct that
> the local sections contains dofs for boundary conditions, while the
> global section does not?  If so, is the following the right plan?
>
> 1. At start, save the DMPlex
> 2. At start, save the local section
> 3. Every frame, sync the global vector a local vector and write out in
> petsc binary format.
>
> Are there standard binary formats for DMPlex and sections?  All I see
> is the following sad commented out call to DMPlexView_Binary.


We are having a big discussion about this on petsc-dev. I will write a
longer
thing, but now the only thing that works all the way is VTK. You should just
be able to create the Viewer, do VecView for your fields, and then destroy
it.

I create a global section with all the boundary values in it.

I plan on having HDF5+Xdmf and Exodus.

   Matt


>
> Geoffrey
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140128/e277a9d9/attachment.html>