This paper presents a new implementation of the dual reciprocity method (DRM) in connection with the dual interpolation boundary face method (DiBFM) for the Poisson equation. In DiBFM, the nodes of an element are categorized into two groups: (i) source nodes (ii) virtual nodes. First layer interpolation is used to interpolate the physical variables, while boundary integrals are evaluated on the source nodes only. Moreover, moving least squares (MLS) interpolation is used and provides additional constraints equations to establish the relationship between source and virtual nodes. Additionally, augmented thin plate spline (ATPS) is used to better interpolate the non-homogeneous term. Finally, it is claimed that the proposed method is much superior to the DRM for Poisson type equation with different geometries, especially for complex geometry. Numerical examples are evaluated and compared with the DRM to ensure the superiority of the proposed method.