#define mxCopy(a,n,e) memcpy(mxCalloc(n,e), a, (n)*(e))

double Rand(void);
double GammaRand(double a);
double BetaRand(double a, double b);
int BinoRand(double p, int n);

double logSum(double a, double b);
double pochhammer(double x, int n);
double di_pochhammer(double x, int n);
double tri_pochhammer(double x, int n);
double gammaln2(double x, double d);
double gammaln(double x);
double digamma(double x);
double trigamma(double x);

unsigned *ismember_sorted(double *a, unsigned a_len, double *s, unsigned s_len);