dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Sun, Apr 26, 2009 at 09:29:14PM +0200, kent-and@xxxxxxxxx wrote:
> 
> There are at least two ways to do this.
> 
> 1) Create a submesh. Anders has implemented support for
>     submeshes but I don't know if it works on meshes that
>     have lower dim than the original. Anders has to tell you how
>     to proceed.

No, a SubMesh is just a part of a Mesh (a new Mesh on a subset of the
cells) with a mapping back to the original mesh. This can be used for
creating matching submeshes when doing for example fluid-structure
interaction.

-- 
Anders


> 2) Change DofMap. As it is now DofMap checks if it has a precomputed
>     dof_map (array of integers). If it has this array it will not compute
>     it. Hence you need to create such an array. The array should contains
>     -1 for the values that should not be used. You create this array by
> running
>     over the part of the mesh that you want to use.
> 
>     I suggest you start experimenting by adding a function to FunctionSpace
>     which again calls DofMap.  This function needs as input the structure
> that
>     you want to use to specify the subdomain.
> 
> OK?
> 
> Kent
> 
> 
> > I had asked a question a while back about having a finite element function
> > that was ONLY defined on the boundary.  If I wanted to put this in, what
> > parts of the code should I look at?  Or is it too complicated?
> >
> > - Shawn
> > _______________________________________________
> > DOLFIN-dev mailing list
> > DOLFIN-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/dolfin-dev
> >
> 
>

Attachment: signature.asc
Description: Digital signature