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 r161294 = 1.0;
        double r161295 = s;
        double r161296 = -r161295;
        double r161297 = exp(r161296);
        double r161298 = r161294 + r161297;
        double r161299 = r161294 / r161298;
        double r161300 = c_p;
        double r161301 = pow(r161299, r161300);
        double r161302 = r161294 - r161299;
        double r161303 = c_n;
        double r161304 = pow(r161302, r161303);
        double r161305 = r161301 * r161304;
        double r161306 = t;
        double r161307 = -r161306;
        double r161308 = exp(r161307);
        double r161309 = r161294 + r161308;
        double r161310 = r161294 / r161309;
        double r161311 = pow(r161310, r161300);
        double r161312 = r161294 - r161310;
        double r161313 = pow(r161312, r161303);
        double r161314 = r161311 * r161313;
        double r161315 = r161305 / r161314;
        return r161315;
}