A fourth order accurate matching operator is constructed on a hexagonal grid, for the interpolation of the mixed boundary value problem of Laplace’s equation, by using the harmonic properties of the solution. With the ...