The first approach is derived by means of the variational formulation of the Poisson problem, while the second one is a direct differential method.
Numerical examples on different domain geometries are implemented to verify and compare the proposed approaches; the computational efficiency of the developed methods is also studied.