fit_sphere.hpp is working, though a printf is needed to avoid NUMERICAL failure