A flexible derivative-free solver for (bound constrained) nonlinear least-squares minimization