Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r3977647 = 1.0;
        double r3977648 = s;
        double r3977649 = -r3977648;
        double r3977650 = exp(r3977649);
        double r3977651 = r3977647 + r3977650;
        double r3977652 = r3977647 / r3977651;
        double r3977653 = c_p;
        double r3977654 = pow(r3977652, r3977653);
        double r3977655 = r3977647 - r3977652;
        double r3977656 = c_n;
        double r3977657 = pow(r3977655, r3977656);
        double r3977658 = r3977654 * r3977657;
        double r3977659 = t;
        double r3977660 = -r3977659;
        double r3977661 = exp(r3977660);
        double r3977662 = r3977647 + r3977661;
        double r3977663 = r3977647 / r3977662;
        double r3977664 = pow(r3977663, r3977653);
        double r3977665 = r3977647 - r3977663;
        double r3977666 = pow(r3977665, r3977656);
        double r3977667 = r3977664 * r3977666;
        double r3977668 = r3977658 / r3977667;
        return r3977668;
}