int bcf_p1_read_prior(bcf_p1aux_t *ma, const char *fn);
long double bcf_p1_cal_g3(bcf_p1aux_t *p1a, double g[3]);
int bcf_p1_set_n1(bcf_p1aux_t *b, int n1);
+ void bcf_p1_set_folded(bcf_p1aux_t *p1a); // only effective when set_n1() is not called
#ifdef __cplusplus
}