#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <float.h>
#include <new>
#include <numeric/nls/lmmin.hh>

Namespaces
	numeric
	Unit headers.

	numeric::nls

Macros
#define	LM_MACHEP 1.2e-16

#define	LM_DWARF 1.0e-38

#define	LM_SQRT_DWARF 3.834e-20

#define	LM_SQRT_GIANT 1.304e19

#define	LM_USERTOL 1.e-14

#define	MIN(a, b) (((a)<=(b)) ? (a) : (b))

#define	MAX(a, b) (((a)>=(b)) ? (a) : (b))

#define	SQR(x) (x)*(x)

Functions
void	numeric::nls::lm_printout_std (int n_par, const double par, int m_dat, const void , const double *fvec, int printflags, int iflag, int iter, int nfev)

void	numeric::nls::lmmin (int n_par, double par, int m_dat, const void data, void(evaluate)(const double par, int m_dat, const void data, double fvec, int info), lm_status_struct status, void(printout)(int n_par, const double par, int m_dat, const void data, const double fvec, int printflags, int iflag, int iter, int nfev))

void	numeric::nls::lm_lmpar (int n, double r, int ldr, int ipvt, double diag, double qtb, double delta, double par, double x, double sdiag, double aux, double *xdi)

void	numeric::nls::lm_qrfac (int m, int n, double a, int pivot, int ipvt, double rdiag, double acnorm, double *wa)

void	numeric::nls::lm_qrsolv (int n, double r, int ldr, int ipvt, double diag, double qtb, double x, double sdiag, double *wa)

void	numeric::nls::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, void(evaluate)(const double par, int m_dat, const void data, double fvec, int info), void(printout)(int n_par, const double par, int m_dat, const void data, const double fvec, int printflags, int iflag, int iter, int nfev), int printflags, const void data)

double	numeric::nls::lm_enorm (int n, const double *x)