dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

What does the mesh connectivity (d, d) really mean? This is related to the bug

    https://bugs.launchpad.net/dolfin/+bug/1063013

The code

    #include <dolfin.h>
    int main()
    {
      // Create mesh
      UnitSquare mesh(1, 1);
      mesh.init(2, 2);
      Cell cell(mesh2, 0);
      for (CellIterator c(cell); !c.end(); ++c)
        cout << "In cell loop" << endl;
    }

prints

    In cell loop
    In cell loop
    In cell loop
    In cell loop
    In cell loop
    In cell loop

In what sense are entities of the same dimension 'connected'? The
present behaviour is causing a problem when computing a BoundaryMesh
in 1D because it picks up the end vertices, but then iterates along
the lines of

        for (VertexIterator v(vertex); !v.end(); ++v)
        {

        }

which then yields the vertex that is one place in from the boundary.

Garth