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 r174425 = 1.0;
        double r174426 = s;
        double r174427 = -r174426;
        double r174428 = exp(r174427);
        double r174429 = r174425 + r174428;
        double r174430 = r174425 / r174429;
        double r174431 = c_p;
        double r174432 = pow(r174430, r174431);
        double r174433 = r174425 - r174430;
        double r174434 = c_n;
        double r174435 = pow(r174433, r174434);
        double r174436 = r174432 * r174435;
        double r174437 = t;
        double r174438 = -r174437;
        double r174439 = exp(r174438);
        double r174440 = r174425 + r174439;
        double r174441 = r174425 / r174440;
        double r174442 = pow(r174441, r174431);
        double r174443 = r174425 - r174441;
        double r174444 = pow(r174443, r174434);
        double r174445 = r174442 * r174444;
        double r174446 = r174436 / r174445;
        return r174446;
}