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 r118734 = 1.0;
        double r118735 = s;
        double r118736 = -r118735;
        double r118737 = exp(r118736);
        double r118738 = r118734 + r118737;
        double r118739 = r118734 / r118738;
        double r118740 = c_p;
        double r118741 = pow(r118739, r118740);
        double r118742 = r118734 - r118739;
        double r118743 = c_n;
        double r118744 = pow(r118742, r118743);
        double r118745 = r118741 * r118744;
        double r118746 = t;
        double r118747 = -r118746;
        double r118748 = exp(r118747);
        double r118749 = r118734 + r118748;
        double r118750 = r118734 / r118749;
        double r118751 = pow(r118750, r118740);
        double r118752 = r118734 - r118750;
        double r118753 = pow(r118752, r118743);
        double r118754 = r118751 * r118753;
        double r118755 = r118745 / r118754;
        return r118755;
}