fit_hbd_pdr_on_best_grid_size function

Fit pulled diversification rates of birth-death models on a time grid with optimal size.