-
Notifications
You must be signed in to change notification settings - Fork 20
/
utils.h
executable file
·110 lines (73 loc) · 2.27 KB
/
utils.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#ifndef UTILSH
#define UTILSH
#include <gsl/gsl_vector.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_sf.h>
#include <gsl/gsl_permutation.h>
#include <math.h>
#include <time.h>
#include <assert.h>
#include <sys/stat.h>
#include <sys/types.h>
#include <stdlib.h>
#include <stdarg.h>
time_t TIME;
char line1[80];
void outlog(char* fmt, ...);
typedef struct int_vector
{
int size;
int* val;
} int_vector;
int_vector* new_int_vector(int size);
void delete_int_vector(int_vector* iv);
void ivappend(int_vector* v, int val);
static inline int ivget(int_vector* v, int n)
{ return(v->val[n]); };
static inline void ivset(int_vector* v, int n, int val)
{ v->val[n] = val; };
static inline void vset(gsl_vector* v, int i, double x)
{ gsl_vector_set(v, i, x); };
static inline double vget(const gsl_vector* v, int i)
{ return(gsl_vector_get(v, i)); };
static inline void vinc(gsl_vector* v, int i, double x)
{ vset(v, i, vget(v, i) + x); };
static inline double log_sum(double log_a, double log_b)
{
if (log_a < log_b)
return(log_b+log(1 + exp(log_a-log_b)));
else
return(log_a+log(1 + exp(log_b-log_a)));
};
static inline double lgam(double x)
{ return(gsl_sf_lngamma(x)); }
double runif();
double rgauss(double mean, double stdev);
int sample_from_log(gsl_vector* log_prob);
void init_random_number_generator();
void vct_fscanf(const char* filename, gsl_vector* v);
void iv_permute(int_vector* iv);
void iv_permute_from_perm(int_vector* iv, gsl_permutation* p);
int_vector* iv_copy(int_vector* iv);
double sum(gsl_vector* v);
void print_vector(gsl_vector* v);
gsl_permutation* rpermutation(int size);
void write_vect(gsl_vector* vect, char* name, FILE* file);
gsl_vector* read_vect(char* name, int size, FILE* file);
int read_int(char* name, FILE* file);
void write_int(int x, char* name, FILE* file);
void write_double(double x, char* name, FILE* file);
double read_double(char* name, FILE* file);
int directory_exist(const char *dname);
void make_directory(char* name);
int rbernoulli(double p);
double rbeta(double a,
double b);
double rgamma(double shape,
double scale);
double log_dgamma(double x,
double shape,
double scale);
void resize(gsl_vector * vec, size_t newsize);
#endif