[aspect-devel] Duplicate Vertices

[Walter Landry]

Wolfgang Bangerth <bangerth at math.tamu.edu> wrote:
> 
>> For the XDMF format, you can store the values at the integration
>> points in HDF5 files and then compute the values at vertices.  I did
>> this for Gale for the pressure.  Doing this kind of computation in the
>> XDMF format is rather painful, but it is possible.
> 
> Interesting. Can you do this for only a subset of the fields that you 
> store in one file?

Yes.  I do it for the P-1 pressure elements, but not the Q2 velocity
elements.  This is because I was reusing the checkpoint data, which
stores the values at the integration points for each cell.  So there
was lots duplicate data for the Q2 elements, but no duplicate data for
the P-1 elements.  Presumably you could be more clever and not
duplicate any Q2 data.

Cheers,
Walter Landry