constraints.add_line(u_i[i]);
constraints.add_line(v_i[i]);
You can only set a constraint for one of the two. Given the next line, I think
you want to set the constraint for u_i[i], but not v_i[i].
constraints.add_entry(u_i[i],v_i[i],0.5);
constraints.set_inhomogeneity(u_i[i],-0.5);
Hello all,
I designed a coupled Laplace problem to play with ConstraintMatrix for
coupled components. On face \Gamma_4 and \Gamma_2, I let u=0.5v-0.5,
v=2u+1, respectively .
[image: Screen Shot 2020-04-05 at 4.14.41 PM.png]
Both u and v use the same finite element: FE_Q<*dim*>(1),