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 r167465 = 1.0;
        double r167466 = s;
        double r167467 = -r167466;
        double r167468 = exp(r167467);
        double r167469 = r167465 + r167468;
        double r167470 = r167465 / r167469;
        double r167471 = c_p;
        double r167472 = pow(r167470, r167471);
        double r167473 = r167465 - r167470;
        double r167474 = c_n;
        double r167475 = pow(r167473, r167474);
        double r167476 = r167472 * r167475;
        double r167477 = t;
        double r167478 = -r167477;
        double r167479 = exp(r167478);
        double r167480 = r167465 + r167479;
        double r167481 = r167465 / r167480;
        double r167482 = pow(r167481, r167471);
        double r167483 = r167465 - r167481;
        double r167484 = pow(r167483, r167474);
        double r167485 = r167482 * r167484;
        double r167486 = r167476 / r167485;
        return r167486;
}