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 r55564454 = 1.0;
        double r55564455 = s;
        double r55564456 = -r55564455;
        double r55564457 = exp(r55564456);
        double r55564458 = r55564454 + r55564457;
        double r55564459 = r55564454 / r55564458;
        double r55564460 = c_p;
        double r55564461 = pow(r55564459, r55564460);
        double r55564462 = r55564454 - r55564459;
        double r55564463 = c_n;
        double r55564464 = pow(r55564462, r55564463);
        double r55564465 = r55564461 * r55564464;
        double r55564466 = t;
        double r55564467 = -r55564466;
        double r55564468 = exp(r55564467);
        double r55564469 = r55564454 + r55564468;
        double r55564470 = r55564454 / r55564469;
        double r55564471 = pow(r55564470, r55564460);
        double r55564472 = r55564454 - r55564470;
        double r55564473 = pow(r55564472, r55564463);
        double r55564474 = r55564471 * r55564473;
        double r55564475 = r55564465 / r55564474;
        return r55564475;
}