Hello Thomas,
sorry, I should have given some more details. My idea is the following:
You could create two "distinct" regions in the triangulation by doubling
the vertices on the interface:
2 3|6 7
+__++__+
| | |
| | |
+__++__+
0 1|4 5
Here the triangulation has 8 vertices, but the coordinates of vertices 1
and 4 resp. 3 and 6 coincide. Then you can use FE_Q (1) on both cells
and the DoFHandler distributes 8 degrees of freedom here.
I admit that this approach is a bit dirty, but it looks more natural to
me than using an hp::FECollection, since it reflects better the idea of
two separated regions connected by an interface. Perhabs someone else
has a better idea.
Best Regards,
Markus
Am 03.02.11 21:28, schrieb [email protected]:
Hi Markus,
how do you mean to declare the interface as a boundary?
Internal boundaries are treated as "boundary" for imposing boundary
conditions, but you don't have double dofs there, isn't it?
Best
Thomas
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
