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 r6993467 = 1.0;
        double r6993468 = s;
        double r6993469 = -r6993468;
        double r6993470 = exp(r6993469);
        double r6993471 = r6993467 + r6993470;
        double r6993472 = r6993467 / r6993471;
        double r6993473 = c_p;
        double r6993474 = pow(r6993472, r6993473);
        double r6993475 = r6993467 - r6993472;
        double r6993476 = c_n;
        double r6993477 = pow(r6993475, r6993476);
        double r6993478 = r6993474 * r6993477;
        double r6993479 = t;
        double r6993480 = -r6993479;
        double r6993481 = exp(r6993480);
        double r6993482 = r6993467 + r6993481;
        double r6993483 = r6993467 / r6993482;
        double r6993484 = pow(r6993483, r6993473);
        double r6993485 = r6993467 - r6993483;
        double r6993486 = pow(r6993485, r6993476);
        double r6993487 = r6993484 * r6993486;
        double r6993488 = r6993478 / r6993487;
        return r6993488;
}