}
}
-#define char2int(s) (((int)s[0])<<8|s[1])
-
static int cal_pdg(const bcf1_t *b, bcf_p1aux_t *ma)
{
int i, j, k;
{
int k, n1 = ma->n1, n2 = ma->n - ma->n1;
long double sum1, sum2;
- pc[0] = pc[1] = pc[2] = pc[3] = 0.;
+ pc[0] = pc[1] = pc[2] = pc[3] = -1.;
if (n1 <= 0 || n2 <= 0) return;
for (k = 0, sum1 = 0.; k <= 2*n1; ++k) sum1 += ma->phi1[k] * ma->z1[k];
for (k = 0, sum2 = 0.; k <= 2*n2; ++k) sum2 += ma->phi2[k] * ma->z2[k];
long double sum = 0.;
// set PL and PL_len
for (i = 0; i < b->n_gi; ++i) {
- if (b->gi[i].fmt == char2int("PL")) {
+ if (b->gi[i].fmt == bcf_str2int("PL", 2)) {
ma->PL = (uint8_t*)b->gi[i].data;
ma->PL_len = b->gi[i].len;
break;