1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
// parameters for calling the high-level interface lmfit
// ( lmfit.c provides lm_initialize_control which sets default values ):
typedef struct {
double ftol; // relative error desired in the sum of squares.
double xtol; // relative error between last two approximations.
double gtol; // orthogonality desired between fvec and its derivs.
double epsilon; // step used to calculate the jacobian.
double stepbound; // initial bound to steps in the outer loop.
double fnorm; // norm of the residue vector fvec.
int maxcall; // maximum number of iterations.
int nfev; // actual number of iterations.
int info; // status of minimization.
} lm_control_type;
// the subroutine that calculates fvec:
typedef void (lm_evaluate_ftype) (
double* par, int m_dat, double* fvec, void *data, int *info );
// default implementation therof, provided by lm_eval.c:
void lm_evaluate_default (
double* par, int m_dat, double* fvec, void *data, int *info );
// the subroutine that informs about fit progress:
typedef void (lm_print_ftype) (
int n_par, double* par, int m_dat, double* fvec, void *data,
int iflag, int iter, int nfev );
// default implementation therof, provided by lm_eval.c:
void lm_print_default (
int n_par, double* par, int m_dat, double* fvec, void *data,
int iflag, int iter, int nfev );
// compact high-level interface:
void lm_initialize_control( lm_control_type *control );
void lm_minimize ( int m_dat, int n_par, double* par,
lm_evaluate_ftype *evaluate, lm_print_ftype *printout,
void *data, lm_control_type *control );
double lm_enorm( int, double* );
// low-level interface for full control:
void lm_lmdif( int m, int n, double* x, double* fvec, double ftol, double xtol,
double gtol, int maxfev, double epsfcn, double* diag, int mode,
double factor, int *info, int *nfev,
double* fjac, int* ipvt, double* qtf,
double* wa1, double* wa2, double* wa3, double* wa4,
lm_evaluate_ftype *evaluate, lm_print_ftype *printout,
void *data );
#ifndef _LMDIF
extern char *lm_infmsg[];
extern char *lm_shortmsg[];
#endif
|
