Cannot sample enough valid points. (more)

\[0.0 \lt c_p \land 0.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 r129587 = 1.0;
        double r129588 = s;
        double r129589 = -r129588;
        double r129590 = exp(r129589);
        double r129591 = r129587 + r129590;
        double r129592 = r129587 / r129591;
        double r129593 = c_p;
        double r129594 = pow(r129592, r129593);
        double r129595 = r129587 - r129592;
        double r129596 = c_n;
        double r129597 = pow(r129595, r129596);
        double r129598 = r129594 * r129597;
        double r129599 = t;
        double r129600 = -r129599;
        double r129601 = exp(r129600);
        double r129602 = r129587 + r129601;
        double r129603 = r129587 / r129602;
        double r129604 = pow(r129603, r129593);
        double r129605 = r129587 - r129603;
        double r129606 = pow(r129605, r129596);
        double r129607 = r129604 * r129606;
        double r129608 = r129598 / r129607;
        return r129608;
}