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 r6708753 = 1.0;
        double r6708754 = s;
        double r6708755 = -r6708754;
        double r6708756 = exp(r6708755);
        double r6708757 = r6708753 + r6708756;
        double r6708758 = r6708753 / r6708757;
        double r6708759 = c_p;
        double r6708760 = pow(r6708758, r6708759);
        double r6708761 = r6708753 - r6708758;
        double r6708762 = c_n;
        double r6708763 = pow(r6708761, r6708762);
        double r6708764 = r6708760 * r6708763;
        double r6708765 = t;
        double r6708766 = -r6708765;
        double r6708767 = exp(r6708766);
        double r6708768 = r6708753 + r6708767;
        double r6708769 = r6708753 / r6708768;
        double r6708770 = pow(r6708769, r6708759);
        double r6708771 = r6708753 - r6708769;
        double r6708772 = pow(r6708771, r6708762);
        double r6708773 = r6708770 * r6708772;
        double r6708774 = r6708764 / r6708773;
        return r6708774;
}