Context Navigation

source: Sophya/trunk/SophyaLib/NTools/nbmath.h@ 633

Visit:

Last change on this file since 633 was 244, checked in by ansari, 26 years ago
beaucoup modifs rz+cmv 22/4/99
File size: 5.1 KB

Line
1	/* quelques fonctions mathematiques et constantes utiles */
2	/* cmv 23/06/94 */
3	#ifndef NBMATH_H_SEEN
4	#define NBMATH_H_SEEN
5
6	#include "machdefs.h"
7
8	#ifdef __cplusplus
9	extern "C" {
10	#endif
11
12	/* EA conflit entre les diverses copies de nbmath.h, nbconst.h... A nettoyer... */
13	#undef Pi
14	#undef Pis2
15	#undef DeuxPi
16	#undef SPi
17	#undef S2Pi
18	#undef Rac2
19	#undef Log2
20	#undef LnPi
21	#undef LgPi
22	#undef Ln10
23	#undef DftoDm
24	#undef Hln2pi
25	#undef JourSec
26	#undef AnSec
27
28	#undef GRAND
29	#undef GRAND2
30	#undef IGRAND
31
32	#define Pi (double) ( 3.14159265358979310e+00) /* c'est Pi */
33	#define Pis2 (double) (1.57079632679489655) /* c'est Pi/2 */
34	#define DeuxPi (double) ( 6.28318530717958620e+00) /* c'est 2Pi /
35	#define SPi (double) ( 1.77245385090551590e+00) /* c'est sqrt(Pi) */
36	#define S2Pi (double) ( 2.50662827463100020e+00) /* c'est sqrt(2Pi) /
37	#define Rac2 (double) ( 1.41421356237309510e+00) /* c'est sqrt(2) */
38	#define Log2 (double) ( 3.01029995663981200e-01) /* c'est log10(2) */
39	#define LnPi (double) ( 1.14472988584940020e+00) /* c'est ln(Pi) */
40	#define LgPi (double) ( 4.97149872694133850e-01) /* c'est log10(Pi) */
41	#define Ln10 (double) ( 2.30258509299404590e+00) /* c'est ln(10) */
42	#define DftoDm (double) ( 1.08573620475812940e+00) /* c'est 2.5/ln(10) */
43	#define Hln2pi (double) ( 9.18938533204672670e-01) /* c'est Ln(2pi)/2 /
44
45	#define JourSec (float) (86400.000000) /* nombre de secondes dans 24H */
46	#define AnSec (int) (31557600) /* nombre de secondes dans 365.25 jours */
47
48	#define GRAND (float) (1.e+35)
49	#define GRAND2 (double) (1.e+35)
50	#define IGRAND (int_4) (2147483647)
51
52	/* pour le fit 1-dim */
53	#define NPAR_FITFUN 10
54	#define FITFUN_FLOAT 1
55	#define FITFUN_DOUBLE 2
56	union FloatDouble { float fx; double dx; };
57	typedef union FloatDouble FLOATDOUBLE;
58
59	double FitFun( double (FunFit) (double ,double ,double *) ,int_4 type
60	,FLOATDOUBLE temps,FLOATDOUBLE mag,FLOATDOUBLE err
61	,int_4 *ndata,int_4 I1,int_4 I2
62	,double par,double epar,double *stepar
63	,double minpar,double maxpar,int_4 npar
64	,double stochi2,int_4 NstepMX,int_4 deb);
65
66	double GausPiv(double A,int_4 nca,int_4 n,double B,int_4 ncb,int_4 m,int_4 Inv);
67
68	int_4 paramga(float sxin, float syin, float rhoin
69	, float smax, float axisrat, float *tiltdeg);
70	int_4 gaparam(float smax, float axisrat, float tiltdeg
71	, float sxin, float syin, float *rhoin);
72
73	double nberfc(double x);
74	float probnb(float sci2,int_4 inddl,int_4 *ipass);
75	double nbgammln(double xx);
76	double nbgammq(double a,double x);
77	double nbgammp(double a,double x);
78
79	void Set_Ihoq(int degre,int mIhoqN,double IhoqNX,double *IhoqNW);
80	int Get_Ihoq(int degre,double IhoqNX,double IhoqNW);
81	void nbgauleg(double x1,double x2,double x,double w,int n);
82	double Integ_Fun(double xmin,double xmax,double (*fonc)(double),int npas);
83	double Integ_Fun_2D(double (*fonc)(double x,double y)
84	,double xmin,double xmax,double ymin,double ymax
85	,int npasx,int npasy);
86
87	void Set_FitFunDPol(int DPol);
88	double Gauss1DPolF(double x,double Par,double DgDpar);
89	double GaussI1DPol(double x,double Par,double DgDpar);
90	double Polyn1D(double x,double Par,double DgDpar);
91
92	double FitProp(double x,double y,double ey,int n,double *a1);
93	double FitLin(double x,double y,double ey,int n
94	,double a0,double a1);
95	double FitPar(double x,double y,double ey,int n
96	,double a0,double a1,double *a2);
97	double FitParLin(double xx,double y,double ey,int n,double x0,int Deg_d
98	,double a0,double a1,double b1,double b2);
99	double FitPropClean(double x,double y,double ey,int n
100	,double per_clean,double *a1);
101	double FitLinClean(double x,double y,double ey,int n
102	,double per_clean,double a0,double a1);
103	double FitParClean(double x,double y,double ey,int n
104	,double per_clean,double a0,double a1,double *a2);
105
106	float interpole(float x,int n,float t1,float t2);
107
108	int_4 MeanSig(float mag,float err,int_4 ndata,float nsigma,float mean,float *sigma);
109	int Most_Probable( float *val, int nval, float binin, int nmoy
110	, float low, float high, float *most, int deb);
111	float Ajust_GaFn(int Ns,float *fcfr
112	,float haut,float ehaut,float mean,float emean
113	,float sigma,float esigma,float fond,float efond
114	,int fixfond,int NBIN_RESOL
115	,float perc_net,float frac_sgb,float dyn_sgb,int deb);
116	float HFit_GaFn(int Ns,float *fcfr
117	,float haut,float ehaut,float mean,float emean
118	,float sigma,float esigma,float fond,float efond
119	,int fixfond,int nbin,float xmin,float xmax,int deb);
120	int Ajust_MnSg(int Ns,float *fcfr
121	,float mean,float emean,float sigma,float esigma
122	,float perc_net,float n_sigma,int deb);
123
124	float Int3DCube(float(*fonction)(float, float, float), float step,
125	float precision);
126
127	int IntFaisDr(int n,float cs,float sn,float p,float x0,float *y0
128	,int npass,float perclean,int lp);
129
130	#ifdef __cplusplus
131	}
132	#endif
133
134	#endif

Note: See TracBrowser for help on using the repository browser.

Download in other formats: