1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
#ifndef MISC_H
#define MISC_H
#include <stdlib.h> /* for size_t */
#define LN_MAX 100
#define LNFACT_MAX 100
double ln(unsigned int n);
double lnfact(unsigned int n);
double kahan_sum(double *x, size_t n);
double interp1d(double x, double *xp, double *yp, size_t n);
int isclose(double a, double b, double rel_tol, double abs_tol);
int allclose(double *a, double *b, size_t n, double rel_tol, double abs_tol);
double logsumexp(double *a, size_t n);
double norm(double x, double mu, double sigma);
double norm_cdf(double x, double mu, double sigma);
#endif
|