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 r145311 = 1.0;
        double r145312 = s;
        double r145313 = -r145312;
        double r145314 = exp(r145313);
        double r145315 = r145311 + r145314;
        double r145316 = r145311 / r145315;
        double r145317 = c_p;
        double r145318 = pow(r145316, r145317);
        double r145319 = r145311 - r145316;
        double r145320 = c_n;
        double r145321 = pow(r145319, r145320);
        double r145322 = r145318 * r145321;
        double r145323 = t;
        double r145324 = -r145323;
        double r145325 = exp(r145324);
        double r145326 = r145311 + r145325;
        double r145327 = r145311 / r145326;
        double r145328 = pow(r145327, r145317);
        double r145329 = r145311 - r145327;
        double r145330 = pow(r145329, r145320);
        double r145331 = r145328 * r145330;
        double r145332 = r145322 / r145331;
        return r145332;
}