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 r139716 = 1.0;
        double r139717 = s;
        double r139718 = -r139717;
        double r139719 = exp(r139718);
        double r139720 = r139716 + r139719;
        double r139721 = r139716 / r139720;
        double r139722 = c_p;
        double r139723 = pow(r139721, r139722);
        double r139724 = r139716 - r139721;
        double r139725 = c_n;
        double r139726 = pow(r139724, r139725);
        double r139727 = r139723 * r139726;
        double r139728 = t;
        double r139729 = -r139728;
        double r139730 = exp(r139729);
        double r139731 = r139716 + r139730;
        double r139732 = r139716 / r139731;
        double r139733 = pow(r139732, r139722);
        double r139734 = r139716 - r139732;
        double r139735 = pow(r139734, r139725);
        double r139736 = r139733 * r139735;
        double r139737 = r139727 / r139736;
        return r139737;
}