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 r5455234 = 1.0;
        double r5455235 = s;
        double r5455236 = -r5455235;
        double r5455237 = exp(r5455236);
        double r5455238 = r5455234 + r5455237;
        double r5455239 = r5455234 / r5455238;
        double r5455240 = c_p;
        double r5455241 = pow(r5455239, r5455240);
        double r5455242 = r5455234 - r5455239;
        double r5455243 = c_n;
        double r5455244 = pow(r5455242, r5455243);
        double r5455245 = r5455241 * r5455244;
        double r5455246 = t;
        double r5455247 = -r5455246;
        double r5455248 = exp(r5455247);
        double r5455249 = r5455234 + r5455248;
        double r5455250 = r5455234 / r5455249;
        double r5455251 = pow(r5455250, r5455240);
        double r5455252 = r5455234 - r5455250;
        double r5455253 = pow(r5455252, r5455243);
        double r5455254 = r5455251 * r5455253;
        double r5455255 = r5455245 / r5455254;
        return r5455255;
}