Fast Integral Equation Solver For Variable Coefficient Elliptic Pdes In Complex Geometries