I also posted recently an example of circle packing with adaptive radius, so the sizing is adjusted inside the goal calculation itself, which I think is a lot faster than what you’d get by putting the whole solver in some external loop.
That was for equal radius, but I think could probably be adapted to adjust a global multiplier for a range of radii from some image or 3d field.