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 r167428 = 1.0;
        double r167429 = s;
        double r167430 = -r167429;
        double r167431 = exp(r167430);
        double r167432 = r167428 + r167431;
        double r167433 = r167428 / r167432;
        double r167434 = c_p;
        double r167435 = pow(r167433, r167434);
        double r167436 = r167428 - r167433;
        double r167437 = c_n;
        double r167438 = pow(r167436, r167437);
        double r167439 = r167435 * r167438;
        double r167440 = t;
        double r167441 = -r167440;
        double r167442 = exp(r167441);
        double r167443 = r167428 + r167442;
        double r167444 = r167428 / r167443;
        double r167445 = pow(r167444, r167434);
        double r167446 = r167428 - r167444;
        double r167447 = pow(r167446, r167437);
        double r167448 = r167445 * r167447;
        double r167449 = r167439 / r167448;
        return r167449;
}