dolfin team mailing list archive

["Garth N. Wells" <g.n.wells@xxxxxxxxxx>]

I'm running into a few problems trying to understand the new mesh
structure. If I do

    UnitSquare mesh(1,1);

    for (FacetIterator facet(mesh); !facet.end(); ++facet)
      cout << "Number of connected cell "
          << (*facet).numConnections(mesh.topology().dim()) << endl;

I get the error

    *** Assertion (entity < num_entities) failed

but it I create a boundary mesh first,

    UnitSquare mesh(1,1);
    MeshFunction<dolfin::uint> vertex_map;
    MeshFunction<dolfin::uint> cell_map;
    BoundaryMesh boundary(mesh, vertex_map, cell_map);

    for (FacetIterator facet(mesh); !facet.end(); ++facet)
      cout << "Number of connected cell "
          << (*facet).numConnections(mesh.topology().dim()) << endl;

everything works as expected. I guess the mesh is not building the list
of facets, and this is being done when creating a boundary mesh. Is this
a bug?

Garth