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 r153377 = 1.0;
        double r153378 = s;
        double r153379 = -r153378;
        double r153380 = exp(r153379);
        double r153381 = r153377 + r153380;
        double r153382 = r153377 / r153381;
        double r153383 = c_p;
        double r153384 = pow(r153382, r153383);
        double r153385 = r153377 - r153382;
        double r153386 = c_n;
        double r153387 = pow(r153385, r153386);
        double r153388 = r153384 * r153387;
        double r153389 = t;
        double r153390 = -r153389;
        double r153391 = exp(r153390);
        double r153392 = r153377 + r153391;
        double r153393 = r153377 / r153392;
        double r153394 = pow(r153393, r153383);
        double r153395 = r153377 - r153393;
        double r153396 = pow(r153395, r153386);
        double r153397 = r153394 * r153396;
        double r153398 = r153388 / r153397;
        return r153398;
}