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
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : dolfin@lists.launchpad.net
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
