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 r7446635 = 1.0;
        double r7446636 = s;
        double r7446637 = -r7446636;
        double r7446638 = exp(r7446637);
        double r7446639 = r7446635 + r7446638;
        double r7446640 = r7446635 / r7446639;
        double r7446641 = c_p;
        double r7446642 = pow(r7446640, r7446641);
        double r7446643 = r7446635 - r7446640;
        double r7446644 = c_n;
        double r7446645 = pow(r7446643, r7446644);
        double r7446646 = r7446642 * r7446645;
        double r7446647 = t;
        double r7446648 = -r7446647;
        double r7446649 = exp(r7446648);
        double r7446650 = r7446635 + r7446649;
        double r7446651 = r7446635 / r7446650;
        double r7446652 = pow(r7446651, r7446641);
        double r7446653 = r7446635 - r7446651;
        double r7446654 = pow(r7446653, r7446644);
        double r7446655 = r7446652 * r7446654;
        double r7446656 = r7446646 / r7446655;
        return r7446656;
}