nr.h

数值算法

typedef struct FCOMPLEX {float r,i;} fcomplex;
typedef struct IMMENSE {unsigned long l,r;} immense;
typedef struct GREAT {unsigned short l,c,r;} great;

#define TRADITIONAL 1   /* but leave LINT_ARGS and ANSI undefined */
#ifdef LINT_ARGS
	void adi(double **, double **, double **, double **, 
		double **, double **, double **, double **, 
		int, int, double, double, double);
	void amoeba(float **, float *, int, float, float (*)(), int *);
	void anneal(float *, float *, int *, int);
	void avevar(float *, int, float *, float *);
	void balanc(float **, int);
	void bcucof(float *, float *, float *, float *, float, 
		float, float **);
	void bcuint(float *, float *, float *, float *, float, 
		float, float, float, float, float, float *, 
		float *, float *);
	float bessi(int, float);
	float bessi0(float);
	float bessi1(float);
	float bessj(int, float);
	float bessj0(float);
	float bessj1(float);
	float bessk(int, float);
	float bessk0(float);
	float bessk1(float);
	float bessy(int, float);
	float bessy0(float);
	float bessy1(float);
	float beta(float, float);
	float betacf(float, float, float);
	float betai(float, float, float);
	float bico(int, int);
	void bksub(int, int, int, int, int, float ***);
	float bnldev(float, int, int *);
	float brent(float, float, float, float (*)(), float, float *);
	void bsstep(float *, float *, int, float *, float,
		float, float *, float *, float *, void (*)());
	void caldat(long, int *, int *, int *);
	float cel(float, float, float, float);
	void chder(float, float, float *, float *, int);
	float chebev(float, float, float *, int, float);
	void chebft(float, float, float *, int, float (*)());
	void chebpc(float *, float *, int);
	void chint(float, float, float *, float *, int);
	void chsone(float *, float *, int, int, float *, 
		float *, float *);
	void chstwo(float *, float *, int, int, float *, 
		float *, float *);
	void cntab1(int **, int, int, float *, float *, 
		float *, float *, float *);
	void cntab2(int **, int, int, float *, float *, 
		float *, float *, float *, float *, float *, float *);
	void convlv(float *, int, float *, int, int, float *);
	void correl(float *, float *, int, float *);
	void cosft(float *, int, int);
	void covsrt(float **, int, int *, int);
	void crank(int, float *, float *);
	float dbrent(float, float, float, float (*)(), float (*)(),
		float, float *);
	void ddpoly(float *, int, float, float *, int);
	void des(immense, immense, int *, int, immense *);
	void ks(immense, int, great *);
	void cyfun(unsigned long, great, unsigned long *);
	float df1dim(float);
	void dfpmin(float *, int, float, int *, float *, float (*)(),
		void (*)());
	void difeq(int, int, int, int, int, int, int *, int, 
		float **, float **);
	void dlinmin(float *, float *, int, float *, float (*)(), void (*)());
	void eclass(int *, int, int *, int *, int);
	void eclazz(int *, int, int (*)());
	void eigsrt(float *, float **, int);
	float el2(float, float, float, float);
	void elmhes(float **, int);
	float erf(float);
	float erfc(float);
	float erfcc(float);
	void eulsum(float *, float, int, float *);
	float evlmem(float, float *, int, float);
	float expdev(int *);
	float f1dim(float);
	float factln(int);
	float factrl(int);
	void fgauss(float, float *, float *, float *, int);
	void fit(float *, float *, int, float *, int, float *, 
		float *, float *, float *, float *, float *);
	void fixrts(float *, int);
	void fleg(float, float *, int);
	void flmoon(int, int, long *, float *);
	void four1(float *, int, int);
	void fourn(float *, int *, int, int);
	void fpoly(float, float *, int);
	void frprmn(float *, int, float, int *, float *, float (*)(),
		void (*)());
	void ftest(float *, int, float *, int, float *, float *);
	float gamdev(int, int *);
	float gammln(float);
	float gammp(float, float);
	float gammq(float, float);
	float gasdev(int *);
	void gauleg(double, double, double *, double *, int);
	void gaussj(float **, int, float **, int);
	void gcf(float *, float, float, float *);
	float golden(float, float, float, float (*)(), float, float *);
	void gser(float *, float, float, float *);
	void hqr(float **, int, float *, float *);
	void hunt(float *, int, float, int *);
	void indexx(int, float *, int *);
	int irbit1(unsigned long int *);
	int irbit2(unsigned long int *);
	void jacobi(float **, int, float *, float **, int *);
	long julday(int, int, int);
	void kendl1(float *, float *, int, float *, float *, float *);
	void kendl2(float **, int, int, float *, float *, float *);
	void ksone(float *, int, float (*)(), float *, float *);
	void kstwo(float *, int, float *, int, float *, float *);
	void laguer(fcomplex *, int, fcomplex *, float, int);
	void lfit(float *, float *, float *, int, float *, int, 
		int *, int, float **, float *, void (*)());
	void linmin(float *, float *, int, float *, float (*)());
	void locate(float *, int, float, int *);
	void lubksb(float **, int, int *, float *);
	void ludcmp(float **, int, int *, float *);
	void mdian1(float *, int, float *);
	void mdian2(float *, int, float *);
	void medfit(float *, float *, int, float *, float *, float *);
	void memcof(float *, int, int, float *, float *);
	float midexp(float (*)(), float, float, int);
	float midinf(float (*)(), float, float, int);
	float midpnt(float (*)(), float, float, int);
	float midsql(float (*)(), float, float, int);
	float midsqu(float (*)(), float, float, int);
	void mmid(float *, float *, int, float, float, int, float *,
		void (*)());
	void mnbrak(float *, float *, float *, float *, 
		float *, float *, float (*)());
	void mnewt(int, float *, int, float, float);
	void moment(float *, int, float *, float *, float *, 
		float *, float *, float *);
	void mprove(float **, float **, int, int *, float *, float *);
	void mrqcof(float *, float *, float *, int, float *, int, 
		int *, int, float **, float *, float *, void (*)());
	void mrqmin(float *, float *, float *, int, float *, int, 
		int *, int, float **, float **, float *, void (*)(),
		float *);
	void odeint(float *, int, float, float, float, float, 
		float, int *, int *, void (*)(), void (*)());
	void pcshft(float, float, float *, int);
	void pearsn(float *, float *, int, float *, float *, float *);
	void piksr2(int, float *, float *);
	void piksrt(int, float *);
	void pinvs(int, int, int, int, int, int, float ***, float **);
	float plgndr(int, int, float);
	float poidev(float, int *);
	void polcoe(float *, float *, int, float *);
	void polcof(float *, float *, int, float *);
	void poldiv(float *, int, float *, int, float *, float *);
	void polin2(float *, float *, float **, int, int, float, 
		float, float *, float *);
	void polint(float *, float *, int, float, float *, float *);
	void powell(float *, float **, int, float, int *, float *, float (*)());
	void predic(float *, int, float *, int, float *, int);
	float probks(float);
	void pzextr(int, float, float *, float *, float *, int, int);
	void qcksrt(int, float *);
	float qgaus(float (*)(float), float, float);
	float qromb(float (*)(float), float, float);
	float qromo(float (*)(float), float, float, float (*)());
	void qroot(float *, int, float *, float *, float);
	float qsimp(float (*)(float), float, float);
	float qtrap(float (*)(float), float, float);
	float quad3d(float (*)(float), float, float);
	float ran0(int *);
	float ran1(int *);
	float ran2(long *);
	float ran3(int *);
	float ran4(int *);
	void rank(int, int *, int *);
	void ratint(float *, float *, int, float, float *, float *);
	void realft(float *, int, int);
	void red(int, int, int, int, int, int, int, int, int, int, 
		int, float ***, float **);
	void rk4(float *, float *, int, float, float, float *, void (*)());
	void rkdumb(float *, int, float, float, int, void (*)());
	void rkqc(float *, float *, int, float *, float, float, 
		float *, float *, float *, void (*)());
	float rofunc(float);
	float rtbis(float (*)(), float, float, float);
	float rtflsp(float (*)(), float, float, float);
	float rtnewt(void (*)(), float, float, float);
	float rtsafe(void (*)(), float, float, float);
	float rtsec(float (*)(), float, float, float);
	void rzextr(int, float, float *, float *, float *, int, int);
	void scrsho(float (*)());
	void shell(int, float *);
	void shoot(int, float *, float *, int, float, float, float, 
		float, float, float *, float *);
	void shootf(int, float *, float *, float *, float *, int, int, 
		float, float, float, float, float, float, float *, 
		float *, float *);
	void simp1(float **, int, int *, int, int, int *, float *);
	void simp2(float **, int, int *, int, int *, int, float *);
	void simp3(float **, int, int, int, int);
	void simplx(float **, int, int, int, int, int, 
		int *, int *, int *);
	void sinft(float *, int);
	void smooft(float *, int, float);
	void sncndn(float, float, float *, float *, float *);
	void solvde(int, float, float, float *, int *, int, int, 
		int, float **, float ***, float **);
	void sor(double **, double **, double **, double **, double **, 
		double **, double **, int, double);
	void sort(int, float *);
	void sort2(int, float *, float *);
	void sort3(int, float *, float *, float *);
	void sparse(float *, int, float *, float *);
	void spctrm(FILE *, float *, int, int, int);
	void spear(float *, float *, int, float *, 
		float *, float *, float *, float *);
	void splie2(float *, float *, float **, int, int, float **);
	void splin2(float *, float *, float **, float **, int, int, 
		float, float, float *);
	void spline(float *, float *, int, float, float, float *);
	void splint(float *, float *, float *, int, float, float *);
	void svbksb(float **, float *, float **, int, int,
		float *, float *);
	void svdcmp(float **, int, int, float *, float **);
	void svdfit(float *, float *, float *, int, float *, int, 
		float **, float **, float *, float *, void (*)());
	void svdvar(float **, int, float *, float **);
	void toeplz(float *, float *, float *, int);
	void tptest(float *, float *, int, float *, float *);
	void tqli(float *, float *, int, float **);
	float trapzd(float (*)(float), float, float, int);
	void tred2(float **, int, float *, float *);
	void tridag(float *, float *, float *, float *, float *, int);
	void ttest(float *, int, float *, int, float *, float *);
	void tutest(float *, int, float *, int, float *, float *);
	void twofft(float *, float *, float *, float *, int);
	void vander(float *, float *, float *, int);
	int zbrac(float (*)(), float *, float *);
	void zbrak(float (*)(), float, float, int, float *, float *, int *);
	float zbrent(float (*)(), float, float, float);
	void zroots(fcomplex *, int, fcomplex *, int);
#endif
#ifdef ANSI
	void  adi(double **a, double **b, double **c, double **d, double **e,
		double **f, double **g, double **u, int jmax, int k,
		double alpha, double beta, double eps);
	void  amoeba(float **p, float *y, int ndim, float ftol,
		float (*funk)(float *), int *iter);
	void  anneal(float *x, float *y, int *iorder, int ncity);
	void  avevar(float *data, int n, float *ave, float *svar);
	void  balanc(float **a, int n);
	void  bcucof(float *y, float *y1, float *y2, float *y12, float d1,
		float d2, float **c);
	void  bcuint(float *y, float *y1, float *y2, float *y12, float x1l, 
		float x1u, float x2l, float x2u, float x1, float x2, 
		float *ansy, float *ansy1, float *ansy2);
	float bessi(int n, float x);
	float bessi0(float x);
	float bessi1(float x);
	float bessj(int n, float x);
	float bessj0(float x);
	float bessj1(float x);
	float bessk(int n, float x);
	float bessk0(float x);
	float bessk1(float x);
	float bessy(int n, float x);
	float bessy0(float x);
	float bessy1(float x);
	float beta(float z, float w);
	float betacf(float a, float b, float x);
	float betai(float a, float b, float x);
	float bico(int n, int k);
	void  bksub(int ne, int nb, int jf, int k1, int k2, float ***c);
	float bnldev(float pp, int n, int *idum);
	float brent(float ax, float bx, float cx, float (*f)(float), float tol,
		float *xmin);
	void  bsstep(float *y, float *dydx, int nv, float *xx, float htry,
		float eps, float *yscal, float *hdid, float *hnext, 
		void (*derivs)(float,float *,float *));
	void  caldat(long julian, int *mm, int *id, int *iyyy);
	float cel(float qqc, float pp, float aa, float bb);
	void  chder(float a, float b, float *c, float *cder, int n);
	float chebev(float a, float b, float *c, int m, float x);
	void  chebft(float a, float b, float *c, int n, float (*func)(float));
	void  chebpc(float *c, float *d, int n);
	void  chint(float a, float b, float *c, float *cint, int n);
	void  chsone(float *bins, float *ebins, int nbins, int knstrn, 
		float *df, float *chsq, float *prob);
	void  chstwo(float *bins1, float *bins2, int nbins, int knstrn, 
		float *df, float *chsq, float *prob);
	void  cntab1(int **nn, int n1, int nj, float *chisq, float *df, 
		float *prob, float *cramrv, float *ccc);
	void  cntab2(int **nn, int ni, int nj, float *h, float *hx, float *hy, 
		float *hygx, float *hxgy, float *uygx, float *uxgy,
		float *uxy);
	void  convlv(float *data, int n, float *respns, int m, int isign, 
		float *ans);
	void  correl(float *data1, float *data2, int n, float *ans);
	void  cosft(float *y, int n, int isign);
	void  covsrt(float **covar, int ma, int *lista, int mfit);
	void  crank(int n, float *w, float *s);
	float dbrent(float ax, float bx, float cx, float (*f)(float),
		float (*df)(float), float tol, float *xmin);
	void  ddpoly(float *c, int nc, float x, float *pd, int nd);
	void  des(immense inp, immense key, int *newkey, int isw, immense *out);
	void  ks(immense key, int n, great *kn);
	void  cyfun(unsigned long ir, great k, unsigned long *iout);
	float df1dim(float x);
	void  dfpmin(float *p, int n, float ftol, int *iter, float *fret, 
		float (*func)(float *), void (*dfunc)(float *,float *));
	void  difeq(int k, int k1, int k2, int jsf, int is1, int isf, 
		int *indexv, int ne, float **s, float **y);
	void  dlinmin(float *p, float *xi, int n, float *fret,
		float (*func)(float *), void (*dfunc)(float *,float *));
	void  eclass(int *nf, int n, int *lista, int *listb, int m);
	void  eclazz(int *nf, int n, int (*equiv)(int,int));
	void  eigsrt(float *d, float **v, int n);
	float el2(float x, float qqc, float aa, float bb);
	void  elmhes(float **a, int n);
	float erf(float x);
	float erfc(float x);
	float erfcc(float x);
	void  eulsum(float *sum, float term, int jterm, float *wksp);
	float evlmem(float fdt, float *cof, int m, float pm);
	float expdev(int *idum);
	float f1dim(float x);
	float factln(int n);
	float factrl(int n);
	void  fgauss(float x, float *a, float *y, float *dyda, int na);
	void  fit(float *x, float *y, int ndata, float *sig, int mwt, float *a, 
		float *b, float *siga, float *sigb, float *chi2, float *q);
	void  fixrts(float *d, int npoles);
	void  fleg(float x, float *pl, int nl);
	void  flmoon(int n, int nph, long *jd, float *frac);
	void  four1(float *data, int nn, int isign);
	void  fourn(float *data, int *nn, int ndim, int isign);
	void  fpoly(float x, float *p, int np);
	void  frprmn(float *p, int n, float ftol, int *iter, float *fret, 
		float (*func)(float *), void (*dfunc)(float *,float *));
	void  ftest(float *data1, int n1, float *data2, int n2, float *f, 
		float *prob);
	float gamdev(int ia, int *idum);
	float gammln(float xx);
	float gammp(float a, float x);
	float gammq(float a, float x);
	float gasdev(int *idum);
	void  gauleg(double x1, double x2, double *x, double *w, int n);
	void  gaussj(float **a, int n, float **b, int m);
	void  gcf(float *gammcf, float a, float x, float *gln);
	float golden(float ax, float bx, float cx, float (*f)(float), float tol, 
		float *xmin);
	void  gser(float *gamser, float a, float x, float *gln);
	void  hqr(float **a, int n, float *wr, float *wi);
	void  hunt(float *xx, int n, float x, int *jlo);
	void  indexx(int n, float *arrin, int *indx);
	int   irbit1(unsigned long int *iseed);
	int   irbit2(unsigned long int *iseed);
	void  jacobi(float **a, int n, float *d, float **v, int *nrot);
	long  julday(int mm, int id, int iyyy);
	void  kendl1(float *data1, float *data2, int n, float *tau, float *z,
		float *prob);
	void  kendl2(float **tab, int i, int j, float *tau, float *z,
		float *prob);
	void  ksone(float *data, int n, float (*func)(float), float *d,