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 r170268 = 1.0;
        double r170269 = s;
        double r170270 = -r170269;
        double r170271 = exp(r170270);
        double r170272 = r170268 + r170271;
        double r170273 = r170268 / r170272;
        double r170274 = c_p;
        double r170275 = pow(r170273, r170274);
        double r170276 = r170268 - r170273;
        double r170277 = c_n;
        double r170278 = pow(r170276, r170277);
        double r170279 = r170275 * r170278;
        double r170280 = t;
        double r170281 = -r170280;
        double r170282 = exp(r170281);
        double r170283 = r170268 + r170282;
        double r170284 = r170268 / r170283;
        double r170285 = pow(r170284, r170274);
        double r170286 = r170268 - r170284;
        double r170287 = pow(r170286, r170277);
        double r170288 = r170285 * r170287;
        double r170289 = r170279 / r170288;
        return r170289;
}