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 r180599 = 1.0;
        double r180600 = s;
        double r180601 = -r180600;
        double r180602 = exp(r180601);
        double r180603 = r180599 + r180602;
        double r180604 = r180599 / r180603;
        double r180605 = c_p;
        double r180606 = pow(r180604, r180605);
        double r180607 = r180599 - r180604;
        double r180608 = c_n;
        double r180609 = pow(r180607, r180608);
        double r180610 = r180606 * r180609;
        double r180611 = t;
        double r180612 = -r180611;
        double r180613 = exp(r180612);
        double r180614 = r180599 + r180613;
        double r180615 = r180599 / r180614;
        double r180616 = pow(r180615, r180605);
        double r180617 = r180599 - r180615;
        double r180618 = pow(r180617, r180608);
        double r180619 = r180616 * r180618;
        double r180620 = r180610 / r180619;
        return r180620;
}