Cannot sample enough valid points. (more)

\[0 \lt c_p \land 0 \lt c_n\]
\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}
double f(double c_p, double c_n, double t, double s) {
        double r42251512 = 1.0;
        double r42251513 = s;
        double r42251514 = -r42251513;
        double r42251515 = exp(r42251514);
        double r42251516 = r42251512 + r42251515;
        double r42251517 = r42251512 / r42251516;
        double r42251518 = c_p;
        double r42251519 = pow(r42251517, r42251518);
        double r42251520 = r42251512 - r42251517;
        double r42251521 = c_n;
        double r42251522 = pow(r42251520, r42251521);
        double r42251523 = r42251519 * r42251522;
        double r42251524 = t;
        double r42251525 = -r42251524;
        double r42251526 = exp(r42251525);
        double r42251527 = r42251512 + r42251526;
        double r42251528 = r42251512 / r42251527;
        double r42251529 = pow(r42251528, r42251518);
        double r42251530 = r42251512 - r42251528;
        double r42251531 = pow(r42251530, r42251521);
        double r42251532 = r42251529 * r42251531;
        double r42251533 = r42251523 / r42251532;
        return r42251533;
}