← Back to team overview

dolfin team mailing list archive

Re: new mesh

 

On Mon, Oct 16, 2006 at 07:54:07PM +0200, Garth N. Wells wrote:
> 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

I'll see if I can fix this.

/Anders


References