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 r4815401 = 1.0;
        double r4815402 = s;
        double r4815403 = -r4815402;
        double r4815404 = exp(r4815403);
        double r4815405 = r4815401 + r4815404;
        double r4815406 = r4815401 / r4815405;
        double r4815407 = c_p;
        double r4815408 = pow(r4815406, r4815407);
        double r4815409 = r4815401 - r4815406;
        double r4815410 = c_n;
        double r4815411 = pow(r4815409, r4815410);
        double r4815412 = r4815408 * r4815411;
        double r4815413 = t;
        double r4815414 = -r4815413;
        double r4815415 = exp(r4815414);
        double r4815416 = r4815401 + r4815415;
        double r4815417 = r4815401 / r4815416;
        double r4815418 = pow(r4815417, r4815407);
        double r4815419 = r4815401 - r4815417;
        double r4815420 = pow(r4815419, r4815410);
        double r4815421 = r4815418 * r4815420;
        double r4815422 = r4815412 / r4815421;
        return r4815422;
}