A Kernel-Free Boundary Integral Method for Variable Coefficients Elliptic PDEs. Commun. Comput. Phys. 2020;15(4):1108-1140. doi:10.4208/cicp.170313.071113s