A short overview on the direct multi-elliptic interpolation and the related meshless methods for solving partial differential equations is given. A new technique is proposed which produces a biharmonic interpolation along the boundary and solves the original problem inside the domain. An error estimation is also derived. To implement the method, quadtree-based multi-level methods are used. The approach avoids the use of large, dense and ill-conditioned matrices and significantly reduces the computational cost. © 2008 Elsevier B.V. All rights reserved.