A High-Accurate Fast Poisson Solver Based on Harmonic Surface Mapping Algorithm (Communications in computational physics, 2021)