Matrix-vector multiplication in weak formulation

The boundary conditions have to be set in your initial guess.