Add a local operator for the Laplace problem

This operator only assembles the Laplace stiffness matrix, and nothing
else.  No volume terms, no Neumann terms, no constraints.

Supported is:
- explicitly assembling the matrix (aka "jacobian_volume")
- matrix-vectors products without explicit matrix assembly
  (aka "alpha_volume")