def test_gaussian_fit_length_default(self): t0 = time.perf_counter() rms_gauss = gaussian_fit(self.time_array, self.gaussian_dist)[2] t1 = time.perf_counter() runtime = t1 - t0 result = (rms_gauss - self.length_gauss) / self.length_gauss return runtime, result
def test_gaussian_fit_length_minimize_trust_constr(self): fitOpt = FitOptions(fittingRoutine='minimize', method='trust-constr') t0 = time.perf_counter() rms_gauss = gaussian_fit(self.time_array, self.gaussian_dist, fitOpt=fitOpt)[2] t1 = time.perf_counter() runtime = t1 - t0 result = (rms_gauss - self.length_gauss) / self.length_gauss return runtime, result
binom_params_gauss = binomial_from_width_ratio(time_array, gaussian_dist, plotOpt=plotOpt, ratio_LUT=new_LUT) print('Gauss: Initial ->', initial_params_gauss, '/ Final ->', _binomial_full_to_rms(*binom_params_gauss[-2:])) # ## 3. Distribution fitting routines # In[11]: from blond_common.fitting.profile import gaussian_fit, parabolic_amplitude_fit, binomial_amplitudeN_fit, PlotOptions fitparams_gauss = gaussian_fit(time_array, gaussian_dist) fitparams_parabamp = parabolic_amplitude_fit(time_array, parabamp_dist) fitparams_binom = binomial_amplitudeN_fit(time_array, binom_dist) print('Gauss: Initial ->', initial_params_gauss, '/ Final ->', fitparams_gauss) print('Parab. amp.: Initial ->', initial_params_parabamp, '/ Final ->', fitparams_parabamp) print('Binomial: Initial ->', initial_params_binom, '/ Final ->', fitparams_binom) # In[12]: plotOpt = PlotOptions(figname='Fit-1', clf=False)