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 r231533 = 1.0;
        double r231534 = s;
        double r231535 = -r231534;
        double r231536 = exp(r231535);
        double r231537 = r231533 + r231536;
        double r231538 = r231533 / r231537;
        double r231539 = c_p;
        double r231540 = pow(r231538, r231539);
        double r231541 = r231533 - r231538;
        double r231542 = c_n;
        double r231543 = pow(r231541, r231542);
        double r231544 = r231540 * r231543;
        double r231545 = t;
        double r231546 = -r231545;
        double r231547 = exp(r231546);
        double r231548 = r231533 + r231547;
        double r231549 = r231533 / r231548;
        double r231550 = pow(r231549, r231539);
        double r231551 = r231533 - r231549;
        double r231552 = pow(r231551, r231542);
        double r231553 = r231550 * r231552;
        double r231554 = r231544 / r231553;
        return r231554;
}