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 r161517 = 1.0;
        double r161518 = s;
        double r161519 = -r161518;
        double r161520 = exp(r161519);
        double r161521 = r161517 + r161520;
        double r161522 = r161517 / r161521;
        double r161523 = c_p;
        double r161524 = pow(r161522, r161523);
        double r161525 = r161517 - r161522;
        double r161526 = c_n;
        double r161527 = pow(r161525, r161526);
        double r161528 = r161524 * r161527;
        double r161529 = t;
        double r161530 = -r161529;
        double r161531 = exp(r161530);
        double r161532 = r161517 + r161531;
        double r161533 = r161517 / r161532;
        double r161534 = pow(r161533, r161523);
        double r161535 = r161517 - r161533;
        double r161536 = pow(r161535, r161526);
        double r161537 = r161534 * r161536;
        double r161538 = r161528 / r161537;
        return r161538;
}