+ for (k = 0, x = 0.; k < 2 * n2; ++k)
+ x += ma->phi[2*n1+k] * ma->z2[k] * expl(lgamma(2*n1 + k + 1) - lgamma(k + 1) + y);
+ pc[1] = ma->z1[2*n1] * x / sum;
+ for (k = 1, x = 0.; k <= 2 * n2; ++k)
+ x += ma->phi[k] * ma->z2[k] * expl(lgamma(ma->M - k + 1) - lgamma(2*n2 - k + 1) + y);
+ pc[1] += ma->z1[0] * x / sum;
+ // the variant is specific to group1
+ y = lgamma(2*n1 + 1) - lgamma(ma->M + 1);
+ for (k = 0, x = 0.; k < 2 * n1; ++k)
+ x += ma->phi[2*n2+k] * ma->z1[k] * expl(lgamma(2*n2 + k + 1) - lgamma(k + 1) + y);
+ pc[0] = ma->z2[2*n2] * x / sum;