-
Notifications
You must be signed in to change notification settings - Fork 1
/
bkgd.h
43 lines (28 loc) · 1005 Bytes
/
bkgd.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
#ifndef __BKGD_H__
#define __BKGD_H__
#include <glib.h>
#include "bkgd_param.h"
#include "bkgd_point.h"
/* amount to scale b values to before rounding to integer */
#define BKGD_SCALE 1000.0
typedef struct {
long start;
long end;
long left_ttl; /* how many conserved sites are to left of this block */
long right_ttl; /* how many conserved sites are to right of this block */
double r; /* rec-rate for block */
double r_start;
double r_end;
} ConsBlock;
double bkgd_t_dist_exp(double t, void *v);
double bkgd_t_dist_gamma(double t, void *v);
double bkgd_blk_integrand(double t, void *v);
double bkgd_site_integrand(double t, void *v);
double bkgd_drv1_blk_integrand(double t, void *parm);
double bkgd_drv1_site_integrand(double t, void *parm);
double bkgd_drv2_blk_integrand(double t, void *parm);
double bkgd_drv2_site_integrand(double t, void *parm);
void bkgd_calc_b(BkgdPoint *bpoint, GList *cons_list,
GList *next_cons, BkgdParam *parm,
double r_chr_len);
#endif