Cannot read MeshFunctions in parallel

[Garth Wells]

Reading MeshFunctions from file fails when running in parallel. Without a rewrite of the XMLMeshFunction, where the mapping between the mesh entities and and the vertices is stored as part of the file format, only vertex and cell valued mesh functions can be used in parallel. Hence, we only address the cases of cell and vertex entities for now.

MeshFunctions over vertices and cells are now working. A nice error is issued when trying to read facets. Closing this bug now. What we need to do is store MeshFunctions in term of the vertices to file in order to reproduce the global numbering in parallel.

[Anders Logg (logg)]

Testing notification from launchpad to dolfin-dev

[Anders Logg (logg)]

On Thu, Sep 03, 2009 at 08:52:50PM -0000, Garth Wells wrote:
> ** Changed in: dolfin
> Assignee: logg (logg) => Skavhaug (skavhaug)

Nice. :-)

--
Anders

[Garth Wells (garth-wells)]

Skavhaug wrote:
> ** Description changed:
>
> - Reading MeshFunctions from file fails when running in parallel
> + Reading MeshFunctions from file fails when running in parallel. Without
> + a rewrite of the XMLMeshFunction, where the mapping between the mesh
> + entities and and the vertices is stored as part of the file format, only
> + vertex and cell valued mesh functions can be used in parallel. Hence, we
> + only address the cases of cell and vertex entities for now.
>

The key to getting a number of demos running is to read facet
MeshFunctions. How much effort will a rewrite require?

Garth

[Anders Logg (logg)]

On Tue, Sep 08, 2009 at 06:38:40PM -0000, Garth Wells wrote:
>
> Skavhaug wrote:
> > ** Description changed:
> >
> > - Reading MeshFunctions from file fails when running in parallel
> > + Reading MeshFunctions from file fails when running in parallel. Without
> > + a rewrite of the XMLMeshFunction, where the mapping between the mesh
> > + entities and and the vertices is stored as part of the file format, only
> > + vertex and cell valued mesh functions can be used in parallel. Hence, we
> > + only address the cases of cell and vertex entities for now.
> >
>
> The key to getting a number of demos running is to read facet
> MeshFunctions. How much effort will a rewrite require?

Quite a bit I think. It involves changing the XML format. I suggest we
move this to 0.9.4 and add a function not_implemented_in_parallel that
we can put in some key places, like in parsing of MeshFunctions. Then
all demos that fail in parallel will give a sensible message about why
they are failing.

--
Anders