A fast and accurate algorithm for the computation of effective electric and mathematically equivalent properties of composites with nonsmooth interfaces is reported. The algorithm is based on an integral equation reformulation of the electrostatic partial differential equation and a fast hierarchical technique for potential field evaluation. In a numerical example, 200 large and strongly inhomogeneous aggregates of randomly overlapping disks are solved with a relative error of 0.0005.