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`
995	ihaka	64
		65	`#ifdef HAVE_IEEEFP_H`
989	maechler	66	`#include <ieeefp.h>`
995	ihaka	67	`#endif`
		68
571	ihaka	69	`extern double m_zero;`
		70	`extern double m_one;`
		71	`extern double m_tiny;`
		72	`#define ML_ERROR(x) /* nothing */`
		73	`#define ML_POSINF (m_one / m_zero)`
		74	`#define ML_NEGINF ((-m_one) / m_zero)`
		75	`#define ML_NAN (m_zero / m_zero)`
		76	`#define ML_UNDERFLOW (m_tiny * m_tiny)`
777	maechler	77	`#define ML_VALID(x) (!isnan(x))`
571	ihaka	78	`#else`
		79	`#define ML_ERROR(x) ml_error(x)`
		80	`#define ML_POSINF DBL_MAX`
		81	`#define ML_NEGINF (-DBL_MAX)`
		82	`#define ML_NAN (-DBL_MAX)`
		83	`#define ML_UNDERFLOW 0`
777	maechler	84	`#define ML_VALID(x) (errno == 0)`
2	r	85	`#endif`
		86
571	ihaka	87	`/* Splus Compatibility */`
2	r	88
571	ihaka	89	`#define snorm norm_rand`
		90	`#define sunif unif_rand`
		91	`#define sexp exp_rand`
2	r	92
571	ihaka	93	`/* Name Hiding to Avoid Clashes with Fortran */`
2	r	94
571	ihaka	95	`#ifdef HIDE_NAMES`
		96	`#define d1mach c_d1mach`
		97	`#define i1mach c_i1mach`
		98	`#endif`
2	r	99
571	ihaka	100	`#define rround fround`
		101	`#define prec fprec`
		102	`#define trunc ftrunc`
829	maechler	103	`/* NO! fsign(.) has 2 arguments; sign(.) has 1..`
		104	`#define sign fsign`
		105	`*/`
2	r	106
571	ihaka	107	`/* Machine Characteristics */`
2	r	108
571	ihaka	109	`double d1mach(int);`
		110	`double d1mach_(int*);`
		111	`int i1mach(int);`
		112	`int i1mach_(int*);`
2	r	113
571	ihaka	114	`/* General Support Functions */`
2	r	115
571	ihaka	116	`int imax2(int, int);`
		117	`int imin2(int, int);`
829	maechler	118	`double sign(double);`
571	ihaka	119	`double fmax2(double, double);`
		120	`double fmin2(double, double);`
		121	`double fmod(double, double);`
		122	`double fprec(double, double);`
		123	`double fround(double, double);`
		124	`double ftrunc(double);`
		125	`double fsign(double, double);`
		126	`double fsquare(double);`
		127	`double fcube(double);`
2	r	128
571	ihaka	129	`/* Random Number Generators */`
2	r	130
571	ihaka	131	`double snorm(void);`
		132	`double sunif(void);`
		133	`double sexp(void);`
2	r	134
571	ihaka	135	`/* Chebyshev Series */`
2	r	136
571	ihaka	137	`int chebyshev_init(double*, int, double);`
		138	`double chebyshev_eval(double, double *, int);`
2	r	139
571	ihaka	140	`/* Gamma and Related Functions */`
2	r	141
571	ihaka	142	`double logrelerr(double);`
		143	`void gammalims(double, double);`
		144	`double lgammacor(double);`
		145	`double gamma(double);`
		146	`double lgamma(double);`
		147	`void dpsifn(double, int, int, int, double, int, int*);`
		148	`double digamma(double);`
		149	`double trigamma(double);`
		150	`double tetragamma(double);`
		151	`double pentagamma(double);`
2	r	152
571	ihaka	153	`double choose(double, double);`
		154	`double lchoose(double, double);`
		155	`double fastchoose(double, double);`
		156	`double lfastchoose(double, double);`
2	r	157
571	ihaka	158	`/* Beta and Related Functions */`
2	r	159
571	ihaka	160	`double beta(double, double);`
		161	`double lbeta(double, double);`
2	r	162
571	ihaka	163	`/* Normal Distribution */`
2	r	164
571	ihaka	165	`double dnorm(double, double, double);`
		166	`double pnorm(double, double, double);`
		167	`double qnorm(double, double, double);`
		168	`double rnorm(double, double);`
2	r	169
571	ihaka	170	`/* Uniform Distribution */`
2	r	171
571	ihaka	172	`double dunif(double, double, double);`
		173	`double punif(double, double, double);`
		174	`double qunif(double, double, double);`
		175	`double runif(double, double);`
2	r	176
571	ihaka	177	`/* Gamma Distribution */`
2	r	178
571	ihaka	179	`double dgamma(double, double, double);`
		180	`double pgamma(double, double, double);`
		181	`double qgamma(double, double, double);`
		182	`double rgamma(double, double);`
2	r	183
571	ihaka	184	`/* Beta Distribution */`
2	r	185
571	ihaka	186	`double dbeta(double, double, double);`
		187	`double pbeta(double, double, double);`
		188	`double pbeta_raw(double, double, double);`
		189	`double qbeta(double, double, double);`
		190	`double rbeta(double, double);`
2	r	191
571	ihaka	192	`/* Lognormal Distribution */`
2	r	193
571	ihaka	194	`double dlnorm(double, double, double);`
		195	`double plnorm(double, double, double);`
		196	`double qlnorm(double, double, double);`
		197	`double rlnorm(double, double);`
2	r	198
571	ihaka	199	`/* Chi-squared Distribution */`
2	r	200
571	ihaka	201	`double dchisq(double, double);`
		202	`double pchisq(double, double);`
		203	`double qchisq(double, double);`
		204	`double rchisq(double);`
2	r	205
571	ihaka	206	`/* Non-central Chi-squared Distribution */`
2	r	207
605	ihaka	208	`double dnchisq(double, double, double);`
571	ihaka	209	`double pnchisq(double, double, double);`
		210	`double qnchisq(double, double, double);`
605	ihaka	211	`double rnchisq(double, double);`
571	ihaka	212
		213	`/* F Distibution */`
		214
		215	`double df(double, double, double);`
		216	`double pf(double, double, double);`
		217	`double qf(double, double, double);`
		218	`double rf(double, double);`
		219
		220	`/* Student t Distibution */`
		221
		222	`double dt(double, double);`
		223	`double pt(double, double);`
		224	`double qt(double, double);`
		225	`double rt(double);`
		226
		227	`/* Binomial Distribution */`
		228
		229	`double dbinom(double, double, double);`
		230	`double pbinom(double, double, double);`
		231	`double qbinom(double, double, double);`
		232	`double rbinom(double, double);`
		233
		234	`/* Cauchy Distribution */`
		235
		236	`double dcauchy(double, double, double);`
		237	`double pcauchy(double, double, double);`
		238	`double qcauchy(double, double, double);`
		239	`double rcauchy(double, double);`
		240
		241	`/* Exponential Distribution */`
		242
		243	`double dexp(double, double);`
		244	`double pexp(double, double);`
		245	`double qexp(double, double);`
		246	`double rexp(double);`
		247
		248	`/* Geometric Distribution */`
		249
		250	`double dgeom(double, double);`
		251	`double pgeom(double, double);`
		252	`double qgeom(double, double);`
		253	`double rgeom(double);`
		254
		255	`/* Hypergeometric Distibution */`
		256
		257	`double dhyper(double, double, double, double);`
		258	`double phyper(double, double, double, double);`
		259	`double qhyper(double, double, double, double);`
		260	`double rhyper(double, double, double);`
		261
		262	`/* Negative Binomial Distribution */`
		263
		264	`double dnbinom(double, double, double);`
		265	`double pnbinom(double, double, double);`
		266	`double qnbinom(double, double, double);`
		267	`double rnbinom(double, double);`
		268
		269	`/* Poisson Distribution */`
		270
		271	`double dpois(double, double);`
		272	`double ppois(double, double);`
		273	`double qpois(double, double);`
		274	`double rpois(double);`
		275
		276	`/* Weibull Distribution */`
		277
		278	`double dweibull(double, double, double);`
		279	`double pweibull(double, double, double);`
		280	`double qweibull(double, double, double);`
		281	`double rweibull(double, double);`
		282
		283	`/* Logistic Distribution */`
		284
		285	`double dlogis(double, double, double);`
		286	`double plogis(double, double, double);`
		287	`double qlogis(double, double, double);`
		288	`double rlogis(double, double);`
		289
602	ihaka	290	`/* Non-central Beta Distribution */`
		291
		292	`double dnbeta(double, double, double, double);`
		293	`double pnbeta(double, double, double, double);`
		294	`double qnbeta(double, double, double, double);`
		295	`double rnbeta(double, double, double);`
		296
		297	`/* Non-central F Distribution */`
		298
		299	`double dnf(double, double, double, double);`
		300	`double pnf(double, double, double, double);`
		301	`double qnf(double, double, double, double);`
		302	`double rnf(double, double, double);`
		303
		304	`/* Non-central Student t Distribution */`
		305
		306	`double dnt(double, double, double);`
		307	`double pnt(double, double, double);`
		308	`double qnt(double, double, double);`
		309	`double rnt(double, double);`
		310
624	ihaka	311	`/* Studentized Range Distribution */`
		312
		313	`double dtukey(double, double, double, double);`
		314	`double ptukey(double, double, double, double);`
		315	`double qtukey(double, double, double, double);`
		316	`double rtukey(double, double, double);`
		317
2	r	318	`#endif`

The R Project SVN R

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