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 r162587 = 1.0;
        double r162588 = s;
        double r162589 = -r162588;
        double r162590 = exp(r162589);
        double r162591 = r162587 + r162590;
        double r162592 = r162587 / r162591;
        double r162593 = c_p;
        double r162594 = pow(r162592, r162593);
        double r162595 = r162587 - r162592;
        double r162596 = c_n;
        double r162597 = pow(r162595, r162596);
        double r162598 = r162594 * r162597;
        double r162599 = t;
        double r162600 = -r162599;
        double r162601 = exp(r162600);
        double r162602 = r162587 + r162601;
        double r162603 = r162587 / r162602;
        double r162604 = pow(r162603, r162593);
        double r162605 = r162587 - r162603;
        double r162606 = pow(r162605, r162596);
        double r162607 = r162604 * r162606;
        double r162608 = r162598 / r162607;
        return r162608;
}