dolfin team mailing list archive

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

I'm struggling a bit with the new mesh structure. For a triangular mesh,
I expected that

  BoundaryMesh boundary(mesh);
  for (CellIterator facet(boundary); !facet.end(); ++facet)
  {
      cout << (*facet).numConnections(0) << endl;
      cout << (*facet).numConnections(2) << endl;
  }

would print 2 and 2 (or 1), but it prints 2 and 8. Also,

  cout << (*facet).numConnections(1) << endl;

leads to a segmentation fault. Should

  foo.numConnections(N);

return the number of mesh entities of topological dimension N which are
attached to foo?

Garth