In DOLFIN, when applying Dirichlet bcs to H(div) spaces, DOLFIN insists that the bc function is a vector-valued function, whereas the physically and mathematically natural function is scalar (normal component). The present state is annoying when boundaries are not axis-aligned.
Does anyone have a nice fix for this, or will it require low-level changes? Looks like the problem is ufc::finite_elemenent::restrict. Garth
_______________________________________________ fenics mailing list fenics@fenicsproject.org http://fenicsproject.org/mailman/listinfo/fenics
