Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r7897294 = 1.0;
        double r7897295 = s;
        double r7897296 = -r7897295;
        double r7897297 = exp(r7897296);
        double r7897298 = r7897294 + r7897297;
        double r7897299 = r7897294 / r7897298;
        double r7897300 = c_p;
        double r7897301 = pow(r7897299, r7897300);
        double r7897302 = r7897294 - r7897299;
        double r7897303 = c_n;
        double r7897304 = pow(r7897302, r7897303);
        double r7897305 = r7897301 * r7897304;
        double r7897306 = t;
        double r7897307 = -r7897306;
        double r7897308 = exp(r7897307);
        double r7897309 = r7897294 + r7897308;
        double r7897310 = r7897294 / r7897309;
        double r7897311 = pow(r7897310, r7897300);
        double r7897312 = r7897294 - r7897310;
        double r7897313 = pow(r7897312, r7897303);
        double r7897314 = r7897311 * r7897313;
        double r7897315 = r7897305 / r7897314;
        return r7897315;
}