  This paper presents a fast wavefield evaluation method for two-dimensional
wave scattering problems. The proposed method is based on a modified version of
proxy-surface-accelerated interpolative decomposition, making it effective even
if the evaluation points are near the boundary. The commonly known fast
multipole method requires the use of direct evaluations near the boundaries of
scatterers because the analytical expansion of kernel functions does not
converge. On the one hand, the proposed method does not require the analytical
expansion of kernel functions. The validity and effectiveness of the proposed
method are demonstrated using numerical examples.
