Added Dirichlet constant source+doublet method
Tested and working for 2D