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 r153600 = 1.0;
        double r153601 = s;
        double r153602 = -r153601;
        double r153603 = exp(r153602);
        double r153604 = r153600 + r153603;
        double r153605 = r153600 / r153604;
        double r153606 = c_p;
        double r153607 = pow(r153605, r153606);
        double r153608 = r153600 - r153605;
        double r153609 = c_n;
        double r153610 = pow(r153608, r153609);
        double r153611 = r153607 * r153610;
        double r153612 = t;
        double r153613 = -r153612;
        double r153614 = exp(r153613);
        double r153615 = r153600 + r153614;
        double r153616 = r153600 / r153615;
        double r153617 = pow(r153616, r153606);
        double r153618 = r153600 - r153616;
        double r153619 = pow(r153618, r153609);
        double r153620 = r153617 * r153619;
        double r153621 = r153611 / r153620;
        return r153621;
}