- ma->theta = theta;
- for (i = 0, sum = 0.; i < 2 * ma->n; ++i)
- sum += (ma->alpha[i] = ma->theta / (2 * ma->n - i));
- ma->alpha[2 * ma->n] = 1. - sum;
+ if (type == MC_PTYPE_COND2) {
+ for (i = 0; i <= 2 * ma->n; ++i)
+ ma->alpha[i] = 2. * (i + 1) / (2 * ma->n + 1) / (2 * ma->n + 2);
+ } else {
+ double sum;
+ for (i = 0, sum = 0.; i < 2 * ma->n; ++i)
+ sum += (ma->alpha[i] = theta / (2 * ma->n - i));
+ ma->alpha[2 * ma->n] = 1. - sum;
+ }