1*5ddc57e5SXin Li /*
2*5ddc57e5SXin Li * Library: lmfit (Levenberg-Marquardt least squares fitting)
3*5ddc57e5SXin Li *
4*5ddc57e5SXin Li * File: lmcurve.c
5*5ddc57e5SXin Li *
6*5ddc57e5SXin Li * Contents: Implements lmcurve, a simplified API for curve fitting
7*5ddc57e5SXin Li * using the generic Levenberg-Marquardt routine lmmin.
8*5ddc57e5SXin Li *
9*5ddc57e5SXin Li * Copyright: Joachim Wuttke, Forschungszentrum Juelich GmbH (2004-2013)
10*5ddc57e5SXin Li *
11*5ddc57e5SXin Li * License: see ../COPYING (FreeBSD)
12*5ddc57e5SXin Li *
13*5ddc57e5SXin Li * Homepage: apps.jcns.fz-juelich.de/lmfit
14*5ddc57e5SXin Li *
15*5ddc57e5SXin Li * Note to programmers: Don't patch and fork, but copy and variate!
16*5ddc57e5SXin Li * If you need to compute residues differently, then please do not patch
17*5ddc57e5SXin Li * lmcurve.c, but copy it to a differently named file, and rename lmcurve(),
18*5ddc57e5SXin Li * lmcurve_evaluate() and lmcurve_data_struct before adapting them to your
19*5ddc57e5SXin Li * needs, like we have done in lmcurve_tyd.c.
20*5ddc57e5SXin Li */
21*5ddc57e5SXin Li
22*5ddc57e5SXin Li #include "lmmin.h"
23*5ddc57e5SXin Li
24*5ddc57e5SXin Li typedef struct {
25*5ddc57e5SXin Li const double* t;
26*5ddc57e5SXin Li const double* y;
27*5ddc57e5SXin Li double (*f)(const double t, const double* par);
28*5ddc57e5SXin Li } lmcurve_data_struct;
29*5ddc57e5SXin Li
lmcurve_evaluate(const double * par,const int m_dat,const void * data,double * fvec,int * info)30*5ddc57e5SXin Li void lmcurve_evaluate(
31*5ddc57e5SXin Li const double* par, const int m_dat, const void* data, double* fvec,
32*5ddc57e5SXin Li int* info)
33*5ddc57e5SXin Li {
34*5ddc57e5SXin Li lmcurve_data_struct* D = (lmcurve_data_struct*)data;
35*5ddc57e5SXin Li int i;
36*5ddc57e5SXin Li for (i = 0; i < m_dat; i++)
37*5ddc57e5SXin Li fvec[i] = D->y[i] - D->f(D->t[i], par);
38*5ddc57e5SXin Li }
39*5ddc57e5SXin Li
lmcurve(const int n_par,double * par,const int m_dat,const double * t,const double * y,double (* f)(const double t,const double * par),const lm_control_struct * control,lm_status_struct * status)40*5ddc57e5SXin Li void lmcurve(
41*5ddc57e5SXin Li const int n_par, double* par, const int m_dat,
42*5ddc57e5SXin Li const double* t, const double* y,
43*5ddc57e5SXin Li double (*f)(const double t, const double* par),
44*5ddc57e5SXin Li const lm_control_struct* control, lm_status_struct* status)
45*5ddc57e5SXin Li {
46*5ddc57e5SXin Li lmcurve_data_struct data = { t, y, f };
47*5ddc57e5SXin Li
48*5ddc57e5SXin Li lmmin(n_par, par, m_dat, (const void*)&data, lmcurve_evaluate,
49*5ddc57e5SXin Li control, status);
50*5ddc57e5SXin Li }
51