The reference files are created with the analytical solution. For the
special setup with a dirac source term at position xi, eta the following
analytical solution in 2 dimensions is valid:

u(x,y) = ln(sqrt((x-xi)^2+(y-eta)^2))/(2 * Pi)