Thanks. I'll try to implement this.
On 15 April 2013 18:20, Anders Logg <[email protected]> wrote: > 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
