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

Reply via email to