Solution of non‐linear boundary integral equations in complex geometries with auxiliary integral subtraction