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 r169669 = 1.0;
        double r169670 = s;
        double r169671 = -r169670;
        double r169672 = exp(r169671);
        double r169673 = r169669 + r169672;
        double r169674 = r169669 / r169673;
        double r169675 = c_p;
        double r169676 = pow(r169674, r169675);
        double r169677 = r169669 - r169674;
        double r169678 = c_n;
        double r169679 = pow(r169677, r169678);
        double r169680 = r169676 * r169679;
        double r169681 = t;
        double r169682 = -r169681;
        double r169683 = exp(r169682);
        double r169684 = r169669 + r169683;
        double r169685 = r169669 / r169684;
        double r169686 = pow(r169685, r169675);
        double r169687 = r169669 - r169685;
        double r169688 = pow(r169687, r169678);
        double r169689 = r169686 * r169688;
        double r169690 = r169680 / r169689;
        return r169690;
}