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 r149280 = 1.0;
        double r149281 = s;
        double r149282 = -r149281;
        double r149283 = exp(r149282);
        double r149284 = r149280 + r149283;
        double r149285 = r149280 / r149284;
        double r149286 = c_p;
        double r149287 = pow(r149285, r149286);
        double r149288 = r149280 - r149285;
        double r149289 = c_n;
        double r149290 = pow(r149288, r149289);
        double r149291 = r149287 * r149290;
        double r149292 = t;
        double r149293 = -r149292;
        double r149294 = exp(r149293);
        double r149295 = r149280 + r149294;
        double r149296 = r149280 / r149295;
        double r149297 = pow(r149296, r149286);
        double r149298 = r149280 - r149296;
        double r149299 = pow(r149298, r149289);
        double r149300 = r149297 * r149299;
        double r149301 = r149291 / r149300;
        return r149301;
}