The R Project SVN – R – Blame – /branches/R-exp-API/src/include/Mathlib.h

Rev	Author	Line No.	Line
571	ihaka	1	`#ifndef MATHLIB_H`
		2	`#define MATHLIB_H`
2	r	3
		4	`#include "Arith.h"`
		5
571	ihaka	6	`#include <errno.h>`
		7	`#include <float.h>`
		8	`#include <limits.h>`
2	r	9	`#include <math.h>`
		10	`#include <stdlib.h>`
		11
777	maechler	12	`/* (>=) 30 Decimal-place constants computed with bc (scale=32; proper round) */`
2	r	13
777	maechler	14	`#ifndef M_SQRT_2`
		15	`#define M_SQRT_2 1.4142135623730950488016887242097`
		16	`#define M_1_SQRT_2 0.707106781186547524400844362105 /* 1/sqrt(2) */`
		17	`#define M_SQRT_32 5.656854249492380195206754896838 /* sqrt(32) */`
2	r	18	`#endif`
		19
777	maechler	20	`#ifndef M_LN_2`
		21	`#define M_LN_2 0.693147180559945309417232121458176568`
		22	`#define M_LOG10_2 0.301029995663981195213738894724493027`
		23	`#endif`
		24
2	r	25	`#ifndef M_PI`
777	maechler	26	`#define M_PI 3.141592653589793238462643383279502884197169399375`
2	r	27	`#endif`
		28	`#ifndef M_PI_half`
777	maechler	29	`#define M_PI_half 1.570796326794896619231321691640`
2	r	30	`#endif`
		31
		32	`#ifndef M_SQRT_PI`
777	maechler	33	`/* sqrt(pi), 1/sqrt(2pi), sqrt(2/pi) : */`
		34	`#define M_SQRT_PI 1.772453850905516027298167483341`
		35	`#define M_1_SQRT_2PI 0.398942280401432677939946059934`
		36	`#define M_SQRT_2dPI 0.79788456080286535587989211986876`
2	r	37	`#endif`
		38
		39
777	maechler	40	`#ifndef M_LN_SQRT_PI`
		41	`/* log(sqrt(pi)) = log(pi)/2 : */`
		42	`#define M_LN_SQRT_PI 0.5723649429247000870717136756765293558`
2	r	43	`/* log(sqrt(2pi)) = log(2pi)/2 : */`
777	maechler	44	`#define M_LN_SQRT_2PI 0.91893853320467274178032973640562`
		45	`/* log(sqrt(pi/2)) = log(pi/2)/2 : */`
		46	`#define M_LN_SQRT_PId2 0.225791352644727432363097614947441`
2	r	47	`#endif`
		48
777	maechler	49
		50
571	ihaka	51	`#define MATHLIB_ERROR(x) { printf("%s\n",x); exit(1); }`
		52	`#define MATHLIB_WARNING(x) { printf("%s\n",x); }`
2	r	53
571	ihaka	54	`#define ME_NONE 0`
		55	`#define ME_DOMAIN 1`
		56	`#define ME_RANGE 2`
		57	`#define ME_NOCONV 3`
		58	`#define ME_PRECISION 4`
		59	`#define ME_UNDERFLOW 5`
		60
		61	`#undef ML_PRECISION_WARNINGS`
		62
		63	`#ifdef IEEE_754`
		64	`extern double m_zero;`
		65	`extern double m_one;`
		66	`extern double m_tiny;`
		67	`#define ML_ERROR(x) /* nothing */`
		68	`#define ML_POSINF (m_one / m_zero)`
		69	`#define ML_NEGINF ((-m_one) / m_zero)`
		70	`#define ML_NAN (m_zero / m_zero)`
		71	`#define ML_UNDERFLOW (m_tiny * m_tiny)`
777	maechler	72	`#define ML_VALID(x) (!isnan(x))`
571	ihaka	73	`#else`
		74	`#define ML_ERROR(x) ml_error(x)`
		75	`#define ML_POSINF DBL_MAX`
		76	`#define ML_NEGINF (-DBL_MAX)`
		77	`#define ML_NAN (-DBL_MAX)`
		78	`#define ML_UNDERFLOW 0`
777	maechler	79	`#define ML_VALID(x) (errno == 0)`
2	r	80	`#endif`
		81
571	ihaka	82	`/* Splus Compatibility */`
2	r	83
571	ihaka	84	`#define snorm norm_rand`
		85	`#define sunif unif_rand`
		86	`#define sexp exp_rand`
2	r	87
571	ihaka	88	`/* Name Hiding to Avoid Clashes with Fortran */`
2	r	89
571	ihaka	90	`#ifdef HIDE_NAMES`
		91	`#define d1mach c_d1mach`
		92	`#define i1mach c_i1mach`
		93	`#endif`
2	r	94
571	ihaka	95	`#define rround fround`
		96	`#define prec fprec`
		97	`#define trunc ftrunc`
829	maechler	98	`/* NO! fsign(.) has 2 arguments; sign(.) has 1..`
		99	`#define sign fsign`
		100	`*/`
2	r	101
571	ihaka	102	`/* Machine Characteristics */`
2	r	103
571	ihaka	104	`double d1mach(int);`
		105	`double d1mach_(int*);`
		106	`int i1mach(int);`
		107	`int i1mach_(int*);`
2	r	108
571	ihaka	109	`/* General Support Functions */`
2	r	110
571	ihaka	111	`int imax2(int, int);`
		112	`int imin2(int, int);`
829	maechler	113	`double sign(double);`
571	ihaka	114	`double fmax2(double, double);`
		115	`double fmin2(double, double);`
		116	`double fmod(double, double);`
		117	`double fprec(double, double);`
		118	`double fround(double, double);`
		119	`double ftrunc(double);`
		120	`double fsign(double, double);`
		121	`double fsquare(double);`
		122	`double fcube(double);`
2	r	123
571	ihaka	124	`/* Random Number Generators */`
2	r	125
571	ihaka	126	`double snorm(void);`
		127	`double sunif(void);`
		128	`double sexp(void);`
2	r	129
571	ihaka	130	`/* Chebyshev Series */`
2	r	131
571	ihaka	132	`int chebyshev_init(double*, int, double);`
		133	`double chebyshev_eval(double, double *, int);`
2	r	134
571	ihaka	135	`/* Gamma and Related Functions */`
2	r	136
571	ihaka	137	`double logrelerr(double);`
		138	`void gammalims(double, double);`
		139	`double lgammacor(double);`
		140	`double gamma(double);`
		141	`double lgamma(double);`
		142	`void dpsifn(double, int, int, int, double, int, int*);`
		143	`double digamma(double);`
		144	`double trigamma(double);`
		145	`double tetragamma(double);`
		146	`double pentagamma(double);`
2	r	147
571	ihaka	148	`double choose(double, double);`
		149	`double lchoose(double, double);`
		150	`double fastchoose(double, double);`
		151	`double lfastchoose(double, double);`
2	r	152
571	ihaka	153	`/* Beta and Related Functions */`
2	r	154
571	ihaka	155	`double beta(double, double);`
		156	`double lbeta(double, double);`
2	r	157
571	ihaka	158	`/* Normal Distribution */`
2	r	159
571	ihaka	160	`double dnorm(double, double, double);`
		161	`double pnorm(double, double, double);`
		162	`double qnorm(double, double, double);`
		163	`double rnorm(double, double);`
2	r	164
571	ihaka	165	`/* Uniform Distribution */`
2	r	166
571	ihaka	167	`double dunif(double, double, double);`
		168	`double punif(double, double, double);`
		169	`double qunif(double, double, double);`
		170	`double runif(double, double);`
2	r	171
571	ihaka	172	`/* Gamma Distribution */`
2	r	173
571	ihaka	174	`double dgamma(double, double, double);`
		175	`double pgamma(double, double, double);`
		176	`double qgamma(double, double, double);`
		177	`double rgamma(double, double);`
2	r	178
571	ihaka	179	`/* Beta Distribution */`
2	r	180
571	ihaka	181	`double dbeta(double, double, double);`
		182	`double pbeta(double, double, double);`
		183	`double pbeta_raw(double, double, double);`
		184	`double qbeta(double, double, double);`
		185	`double rbeta(double, double);`
2	r	186
571	ihaka	187	`/* Lognormal Distribution */`
2	r	188
571	ihaka	189	`double dlnorm(double, double, double);`
		190	`double plnorm(double, double, double);`
		191	`double qlnorm(double, double, double);`
		192	`double rlnorm(double, double);`
2	r	193
571	ihaka	194	`/* Chi-squared Distribution */`
2	r	195
571	ihaka	196	`double dchisq(double, double);`
		197	`double pchisq(double, double);`
		198	`double qchisq(double, double);`
		199	`double rchisq(double);`
2	r	200
571	ihaka	201	`/* Non-central Chi-squared Distribution */`
2	r	202
605	ihaka	203	`double dnchisq(double, double, double);`
571	ihaka	204	`double pnchisq(double, double, double);`
		205	`double qnchisq(double, double, double);`
605	ihaka	206	`double rnchisq(double, double);`
571	ihaka	207
		208	`/* F Distibution */`
		209
		210	`double df(double, double, double);`
		211	`double pf(double, double, double);`
		212	`double qf(double, double, double);`
		213	`double rf(double, double);`
		214
		215	`/* Student t Distibution */`
		216
		217	`double dt(double, double);`
		218	`double pt(double, double);`
		219	`double qt(double, double);`
		220	`double rt(double);`
		221
		222	`/* Binomial Distribution */`
		223
		224	`double dbinom(double, double, double);`
		225	`double pbinom(double, double, double);`
		226	`double qbinom(double, double, double);`
		227	`double rbinom(double, double);`
		228
		229	`/* Cauchy Distribution */`
		230
		231	`double dcauchy(double, double, double);`
		232	`double pcauchy(double, double, double);`
		233	`double qcauchy(double, double, double);`
		234	`double rcauchy(double, double);`
		235
		236	`/* Exponential Distribution */`
		237
		238	`double dexp(double, double);`
		239	`double pexp(double, double);`
		240	`double qexp(double, double);`
		241	`double rexp(double);`
		242
		243	`/* Geometric Distribution */`
		244
		245	`double dgeom(double, double);`
		246	`double pgeom(double, double);`
		247	`double qgeom(double, double);`
		248	`double rgeom(double);`
		249
		250	`/* Hypergeometric Distibution */`
		251
		252	`double dhyper(double, double, double, double);`
		253	`double phyper(double, double, double, double);`
		254	`double qhyper(double, double, double, double);`
		255	`double rhyper(double, double, double);`
		256
		257	`/* Negative Binomial Distribution */`
		258
		259	`double dnbinom(double, double, double);`
		260	`double pnbinom(double, double, double);`
		261	`double qnbinom(double, double, double);`
		262	`double rnbinom(double, double);`
		263
		264	`/* Poisson Distribution */`
		265
		266	`double dpois(double, double);`
		267	`double ppois(double, double);`
		268	`double qpois(double, double);`
		269	`double rpois(double);`
		270
		271	`/* Weibull Distribution */`
		272
		273	`double dweibull(double, double, double);`
		274	`double pweibull(double, double, double);`
		275	`double qweibull(double, double, double);`
		276	`double rweibull(double, double);`
		277
		278	`/* Logistic Distribution */`
		279
		280	`double dlogis(double, double, double);`
		281	`double plogis(double, double, double);`
		282	`double qlogis(double, double, double);`
		283	`double rlogis(double, double);`
		284
602	ihaka	285	`/* Non-central Beta Distribution */`
		286
		287	`double dnbeta(double, double, double, double);`
		288	`double pnbeta(double, double, double, double);`
		289	`double qnbeta(double, double, double, double);`
		290	`double rnbeta(double, double, double);`
		291
		292	`/* Non-central F Distribution */`
		293
		294	`double dnf(double, double, double, double);`
		295	`double pnf(double, double, double, double);`
		296	`double qnf(double, double, double, double);`
		297	`double rnf(double, double, double);`
		298
		299	`/* Non-central Student t Distribution */`
		300
		301	`double dnt(double, double, double);`
		302	`double pnt(double, double, double);`
		303	`double qnt(double, double, double);`
		304	`double rnt(double, double);`
		305
624	ihaka	306	`/* Studentized Range Distribution */`
		307
		308	`double dtukey(double, double, double, double);`
		309	`double ptukey(double, double, double, double);`
		310	`double qtukey(double, double, double, double);`
		311	`double rtukey(double, double, double);`
		312
2	r	313	`#endif`

The R Project SVN R

(root)/branches/R-exp-API/src/include/Mathlib.h @ 28362 – Rev 829