Fast, High-Order Accurate Integral Equation Methods And Application To Pde-Constrained Optimization