I know that we can get the stiffness matrix and external forces through ``assemble_system'' and then solve the augmented equation with a supplemented constraint.
Obviously, it is not a fenics approach. Can we model it as a variational problem with constraints?
Any hint is appreciated!
