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 r6869612 = 1.0;
        double r6869613 = s;
        double r6869614 = -r6869613;
        double r6869615 = exp(r6869614);
        double r6869616 = r6869612 + r6869615;
        double r6869617 = r6869612 / r6869616;
        double r6869618 = c_p;
        double r6869619 = pow(r6869617, r6869618);
        double r6869620 = r6869612 - r6869617;
        double r6869621 = c_n;
        double r6869622 = pow(r6869620, r6869621);
        double r6869623 = r6869619 * r6869622;
        double r6869624 = t;
        double r6869625 = -r6869624;
        double r6869626 = exp(r6869625);
        double r6869627 = r6869612 + r6869626;
        double r6869628 = r6869612 / r6869627;
        double r6869629 = pow(r6869628, r6869618);
        double r6869630 = r6869612 - r6869628;
        double r6869631 = pow(r6869630, r6869621);
        double r6869632 = r6869629 * r6869631;
        double r6869633 = r6869623 / r6869632;
        return r6869633;
}