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 r116630 = 1.0;
        double r116631 = s;
        double r116632 = -r116631;
        double r116633 = exp(r116632);
        double r116634 = r116630 + r116633;
        double r116635 = r116630 / r116634;
        double r116636 = c_p;
        double r116637 = pow(r116635, r116636);
        double r116638 = r116630 - r116635;
        double r116639 = c_n;
        double r116640 = pow(r116638, r116639);
        double r116641 = r116637 * r116640;
        double r116642 = t;
        double r116643 = -r116642;
        double r116644 = exp(r116643);
        double r116645 = r116630 + r116644;
        double r116646 = r116630 / r116645;
        double r116647 = pow(r116646, r116636);
        double r116648 = r116630 - r116646;
        double r116649 = pow(r116648, r116639);
        double r116650 = r116647 * r116649;
        double r116651 = r116641 / r116650;
        return r116651;
}