示例#1
0
文件: LMfit.cpp 项目: eugeneai/fityk
void LMfit::autoiter()
{
    wssr_before_ = chi2;
    realt prev_chi2 = chi2;
    const SettingsMgr *sm = F_->settings_mgr();
    if (F_->get_verbosity() >= 1) {
        F_->ui()->mesg("\t === Levenberg-Marquardt method ===");
        F_->ui()->mesg("Initial values:  lambda=" + S(lambda) +
                       "  WSSR=" + sm->format_double(chi2));
        F_->ui()->mesg("Max. number of iterations: " + max_iterations_);
    }
    realt stop_rel = F_->get_settings()->lm_stop_rel_change;
    realt max_lambda = F_->get_settings()->lm_max_lambda;
    if (stop_rel > 0 && F_->get_verbosity() >= 1)
        F_->ui()->mesg("Will stop when relative change of WSSR is "
                       "twice in row below " + S (stop_rel * 100.) + "%");
    int small_change_counter = 0;
    for (int iter = 0; !common_termination_criteria(iter); iter++) {
        bool better_fit = do_iteration();
        if (better_fit) {
            realt d = prev_chi2 - chi2;
            if (F_->get_verbosity() >= 1)
                F_->ui()->mesg("#" + S(iter_nr_) + ":"
                               "  WSSR=" + sm->format_double(chi2) +
                               "  lambda=" + S(lambda) + "  d(WSSR)=" + S(-d) +
                               "  (" + S (d / prev_chi2 * 100) + "%)");
            // another termination criterium: negligible change of chi2
            if (d / prev_chi2 < stop_rel || chi2 == 0) {
                small_change_counter++;
                if (small_change_counter >= 2 || chi2 == 0) {
                    F_->msg("... converged.");
                    break;
                }
            }
            else
                small_change_counter = 0;
            prev_chi2 = chi2;
        }
        else { // no better fit
            if (F_->get_verbosity() >= 1)
                F_->ui()->mesg("#" + S(iter_nr_) + ":"
                               " (WSSR=" + sm->format_double(chi2_) + ")"
                               "  lambda=" + S(lambda));
            if (lambda > max_lambda) { // another termination criterium
                F_->msg("In L-M method: lambda=" + S(lambda) + " > "
                        + S(max_lambda) + ", stopped.");
                break;
            }
        }
        iteration_plot(a, chi2);
    }
    post_fit (a, chi2);
}
示例#2
0
文件: NMfit.cpp 项目: darckense/fityk
void NMfit::autoiter()
{
    realt convergence = F_->get_settings()->nm_convergence;
    wssr_before_ = compute_wssr(a_orig_, dmdm_);
    F_->msg("WSSR before starting simplex fit: " + S(wssr_before_));
    for (int iter = 0; !termination_criteria(iter, convergence); ++iter) {
        ++iter_nr_;
        change_simplex();
        find_best_worst();
        iteration_plot(best->a, best->wssr);
    }
    post_fit (best->a, best->wssr);
}