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 r261396 = 1.0;
        double r261397 = s;
        double r261398 = -r261397;
        double r261399 = exp(r261398);
        double r261400 = r261396 + r261399;
        double r261401 = r261396 / r261400;
        double r261402 = c_p;
        double r261403 = pow(r261401, r261402);
        double r261404 = r261396 - r261401;
        double r261405 = c_n;
        double r261406 = pow(r261404, r261405);
        double r261407 = r261403 * r261406;
        double r261408 = t;
        double r261409 = -r261408;
        double r261410 = exp(r261409);
        double r261411 = r261396 + r261410;
        double r261412 = r261396 / r261411;
        double r261413 = pow(r261412, r261402);
        double r261414 = r261396 - r261412;
        double r261415 = pow(r261414, r261405);
        double r261416 = r261413 * r261415;
        double r261417 = r261407 / r261416;
        return r261417;
}