A linearized circle packing algorithm

This paper presents a geometric algorithm for approximating radii and centers for a variety of univalent circle packings, including maximal circle packings on the unit disc and the sphere and certain polygonal circle packings in the plane. This method involves an iterative process which alternates between estimates of circle radii and locations of circle centers. The algorithm employs sparse linear systems and in practice achieves a consistent linear convergence rate that is far superior to traditional packing methods. It is deployed in a MATLAB® package which is freely available. This paper gives background on circle packing, a description of the linearized algorithm, illustrations of its use, sample performance data, and remaining challenges.