Now:
* Constructor takes no argument.
* Assembly additionally takes FeSpace.

This way, the [Bi]linearForm can be recycled and assembled for multiple
FeSpaces, without need to set the coefficients every time (e.g., for UQ,
where the same problem is solved on multiple meshes, or assembly of
auxiliary matrices in multigrid solvers).