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.
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