#ifndef MISC_H
#define MISC_H

#include <stdlib.h> /* for size_t */

#define LNFACT_MAX 100

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