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 r127542 = 1.0;
        double r127543 = s;
        double r127544 = -r127543;
        double r127545 = exp(r127544);
        double r127546 = r127542 + r127545;
        double r127547 = r127542 / r127546;
        double r127548 = c_p;
        double r127549 = pow(r127547, r127548);
        double r127550 = r127542 - r127547;
        double r127551 = c_n;
        double r127552 = pow(r127550, r127551);
        double r127553 = r127549 * r127552;
        double r127554 = t;
        double r127555 = -r127554;
        double r127556 = exp(r127555);
        double r127557 = r127542 + r127556;
        double r127558 = r127542 / r127557;
        double r127559 = pow(r127558, r127548);
        double r127560 = r127542 - r127558;
        double r127561 = pow(r127560, r127551);
        double r127562 = r127559 * r127561;
        double r127563 = r127553 / r127562;
        return r127563;
}