It seems I had a misunderstanding about int2d
suppose that I have a mixed fespace
like
fespace -> [RT0Orth,P1]
I like to have RT0Orth
to be defined all over the domain and P1
only inside the coil
.
In fact I like to make a matrix like K=[A B; B' C]
which
A=(RT0Orth,RT0Orth);
B=(RT0Orth,P1);
C=(P1,P1);