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 r102558 = 1.0;
        double r102559 = s;
        double r102560 = -r102559;
        double r102561 = exp(r102560);
        double r102562 = r102558 + r102561;
        double r102563 = r102558 / r102562;
        double r102564 = c_p;
        double r102565 = pow(r102563, r102564);
        double r102566 = r102558 - r102563;
        double r102567 = c_n;
        double r102568 = pow(r102566, r102567);
        double r102569 = r102565 * r102568;
        double r102570 = t;
        double r102571 = -r102570;
        double r102572 = exp(r102571);
        double r102573 = r102558 + r102572;
        double r102574 = r102558 / r102573;
        double r102575 = pow(r102574, r102564);
        double r102576 = r102558 - r102574;
        double r102577 = pow(r102576, r102567);
        double r102578 = r102575 * r102577;
        double r102579 = r102569 / r102578;
        return r102579;
}