[Paraview] VTK parallel unstructured grid and halo elements

[Louie Cardone-Noott]

Hi Paraview list,

I use a finite element code that generates VTK output, but I'm having
trouble dealing with halo (ghost) nodes/elements correctly.

We currently partition meshes node-wise, and each process writes all its
owned nodes and one level of halo nodes to a chunk. Once all the chunks are
loaded there is a "double layer" of internal boundaries between processes.
I'd really like to be able to use the Extract Surface filter correctly but
as you would expect it also finds all the internal surfaces. From searching
I've tried:

* Clean To Grid, but this removes both sides of the double layer leaving
unsightly gaps.
* Writing a vtkGhostPoints array (on nodes) but this seems not to be used
by this filter (or Paraview at all?)
* Writing a vtkGhostLevels array (on elements). This correctly hides all
the ghost elements at the internal boundaries (great!) but still doesn't
let me do everything I want to. For example, the extracted surface looks
right, but the array limits still seem to be based on the max/min of the
external AND internal surfaces, so I suspect the internal ones are just
"hidden from view".

What's the right way to deal with this?

Thanks in advance.

Louie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20150117/3a8b8d54/attachment.html>