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 r3028581 = 1.0;
        double r3028582 = s;
        double r3028583 = -r3028582;
        double r3028584 = exp(r3028583);
        double r3028585 = r3028581 + r3028584;
        double r3028586 = r3028581 / r3028585;
        double r3028587 = c_p;
        double r3028588 = pow(r3028586, r3028587);
        double r3028589 = r3028581 - r3028586;
        double r3028590 = c_n;
        double r3028591 = pow(r3028589, r3028590);
        double r3028592 = r3028588 * r3028591;
        double r3028593 = t;
        double r3028594 = -r3028593;
        double r3028595 = exp(r3028594);
        double r3028596 = r3028581 + r3028595;
        double r3028597 = r3028581 / r3028596;
        double r3028598 = pow(r3028597, r3028587);
        double r3028599 = r3028581 - r3028597;
        double r3028600 = pow(r3028599, r3028590);
        double r3028601 = r3028598 * r3028600;
        double r3028602 = r3028592 / r3028601;
        return r3028602;
}