This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

Subdomainwise discontinuous P1 space?

+1 vote

Given a non-overlapping partition $\Omega = \overline{\Omega_1 \cup \Omega_2}$, is there any chance to obtain the following function space in FEniCS?
$$ Q := \{ v \in L^2(\Omega)~~:~~v|_T \in P_1(T)~~\forall~T~~:~~v|_{\Omega_i} \in C(\Omega_i),~~i=1,2\}$$

In other words, I want to have the piecewise $H^1$-conforming linear finite element space, but I want it to be discontinuous at the interface $\Gamma := \partial\Omega_1\cap\partial\Omega_2$.

btw, I am using the Python interface.

asked Mar 20, 2014 by Christian Waluga FEniCS Expert (12,310 points)
edited Mar 20, 2014 by Christian Waluga

1 Answer

+2 votes
 
Best answer

Simple answer: no, in future, with multi-domain functionality.

Workaround: define

V = VectorFunctionSpace(mesh, 'CG', 1, dim=2)

and restrict redundant DOFs using DirichletBC.

answered Mar 20, 2014 by Jan Blechta FEniCS Expert (51,420 points)
selected Mar 21, 2014 by Christian Waluga

Okay, I think I got it now. Far from convenient, but at least it works. Can't wait for the multi-domain stuff :-). Thanks!

...