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.
PointSource