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 r185589 = 1.0;
        double r185590 = s;
        double r185591 = -r185590;
        double r185592 = exp(r185591);
        double r185593 = r185589 + r185592;
        double r185594 = r185589 / r185593;
        double r185595 = c_p;
        double r185596 = pow(r185594, r185595);
        double r185597 = r185589 - r185594;
        double r185598 = c_n;
        double r185599 = pow(r185597, r185598);
        double r185600 = r185596 * r185599;
        double r185601 = t;
        double r185602 = -r185601;
        double r185603 = exp(r185602);
        double r185604 = r185589 + r185603;
        double r185605 = r185589 / r185604;
        double r185606 = pow(r185605, r185595);
        double r185607 = r185589 - r185605;
        double r185608 = pow(r185607, r185598);
        double r185609 = r185606 * r185608;
        double r185610 = r185600 / r185609;
        return r185610;
}