dolfin team mailing list archive

["Ostien, Jakob T" <jtostie@xxxxxxxxxx>]

Hi,

I have a MeshFunction that I am using to assemble SubDomains, but when I put it into the constructor of my NonlinearProblem, it trips the assert in MeshFunction.h line 90 that says

  dolfin_assert(&entity.mesh() == _mesh);

And this is puzzling because when I add two statements to see what is different about the locally owned mesh and entity mesh, e.g.

  (&entity.mesh())->disp();
  _mesh->disp();

The results _are_ identical.

Has anyone else seen this and/or know what to do about it?

Jake