boundary = BoundaryMesh(mesh, "exterior")
for c0 in cells(boundary):
print c0
for c1 in cells(c0):
print " ", c1
This will iterate over all the cells on the boundary which are the
boundary facets of the original mesh.
Use boundary.entity_map to get the mapping to the original mesh.
Note that "exterior" has nothing to do with your interior/exterior
surfaces.
--
Anders
On Mon, Apr 15, 2013 at 06:12:08PM +0100, Gennadiy Rishkin wrote:
> Hi Anders,
> I can iterate over all the faces on say the exterior but how do I get
> to a neighbour facet from the current facet so I can mark it?
> Thanks,
> Gennadiy
>
> On 15 April 2013 17:49, Gennadiy Rishkin
> <[1][email protected]> wrote:
>
> Thanks. I'll try this.
> Gennadiy
>
> On 15 April 2013 16:33, Anders Logg <[2][email protected]> wrote:
>
> ok, so it has a thickness?
> My suggestion would be to:
> 1. mark boundary facets with a number, say 0
> 2. manually identify *one* facet on either the interior or exterior
> boundary
> 3. recursively iterate outward from that facet (neigbors of
> neighbors
> etc) and mark by 1
> This can be done fairly easily with mesh iterators, but there's no
> built-in functionality for it.
>
> On Mon, Apr 15, 2013 at 04:28:12PM +0100, Gennadiy Rishkin wrote:
> > So, if we assume the rugby ball wasn't deformed, it would be
> shaped
> > like a hollow ellipsoid. The outer surface would be the external
> > surface that is kicked by the foot, while the inner surface would
> be
> > inner lining. Unfortunately, the mesh is not a perfect ellipsoid.
> > Gennadiy
> >
>
> > On 15 April 2013 16:23, Anders Logg <[1][3][email protected]> wrote:
> >
> > On Mon, Apr 15, 2013 at 04:08:50PM +0100, Gennadiy Rishkin wrote:
> > > Hello,
> > > I have a mesh of a slightly or somewhat deformed rugby ball
> after
> > being
> > > pummelled during a game. I do not have marked surfaces but I
> would
> > like
> > > to enumerate the inner surface and the outer surface
> separately.
> > Is
> > > there a way I can do this in FEniCS?
> > > Thanks in advance,
> > > Gennadiy
> >
> > What do you mean by inner and outer surface?
> >
>
> > Referenser
> >
> > 1. mailto:[4][email protected]
>
> Referenser
>
> 1. mailto:[email protected]
> 2. mailto:[email protected]
> 3. mailto:[email protected]
> 4. mailto:[email protected]
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics