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 r165387 = 1.0;
        double r165388 = s;
        double r165389 = -r165388;
        double r165390 = exp(r165389);
        double r165391 = r165387 + r165390;
        double r165392 = r165387 / r165391;
        double r165393 = c_p;
        double r165394 = pow(r165392, r165393);
        double r165395 = r165387 - r165392;
        double r165396 = c_n;
        double r165397 = pow(r165395, r165396);
        double r165398 = r165394 * r165397;
        double r165399 = t;
        double r165400 = -r165399;
        double r165401 = exp(r165400);
        double r165402 = r165387 + r165401;
        double r165403 = r165387 / r165402;
        double r165404 = pow(r165403, r165393);
        double r165405 = r165387 - r165403;
        double r165406 = pow(r165405, r165396);
        double r165407 = r165404 * r165406;
        double r165408 = r165398 / r165407;
        return r165408;
}