The terms related to the Lagrange multipliers are badly handled in the ImposedDrivingVariable:
For the second point, one must do this
- r(this->c) = a*u1(pos);
+ r(this->c) -= a*u1(pos);
Those mistakes don't have any impact on the results, but changes the values of the residual.
Fixed in current revisions of
rliv-2.0,rliv-3.0andtrunk.