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 r165421 = 1.0;
        double r165422 = s;
        double r165423 = -r165422;
        double r165424 = exp(r165423);
        double r165425 = r165421 + r165424;
        double r165426 = r165421 / r165425;
        double r165427 = c_p;
        double r165428 = pow(r165426, r165427);
        double r165429 = r165421 - r165426;
        double r165430 = c_n;
        double r165431 = pow(r165429, r165430);
        double r165432 = r165428 * r165431;
        double r165433 = t;
        double r165434 = -r165433;
        double r165435 = exp(r165434);
        double r165436 = r165421 + r165435;
        double r165437 = r165421 / r165436;
        double r165438 = pow(r165437, r165427);
        double r165439 = r165421 - r165437;
        double r165440 = pow(r165439, r165430);
        double r165441 = r165438 * r165440;
        double r165442 = r165432 / r165441;
        return r165442;
}