Implement point source in the variational formulation

MiroK has outlined the code for a mollified Dirac delta function. See here.

For the precise formulation of a point source in your variational formulation, and FEniCS code, use the PointSource class documented here.