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 r6970327 = 1.0;
        double r6970328 = s;
        double r6970329 = -r6970328;
        double r6970330 = exp(r6970329);
        double r6970331 = r6970327 + r6970330;
        double r6970332 = r6970327 / r6970331;
        double r6970333 = c_p;
        double r6970334 = pow(r6970332, r6970333);
        double r6970335 = r6970327 - r6970332;
        double r6970336 = c_n;
        double r6970337 = pow(r6970335, r6970336);
        double r6970338 = r6970334 * r6970337;
        double r6970339 = t;
        double r6970340 = -r6970339;
        double r6970341 = exp(r6970340);
        double r6970342 = r6970327 + r6970341;
        double r6970343 = r6970327 / r6970342;
        double r6970344 = pow(r6970343, r6970333);
        double r6970345 = r6970327 - r6970343;
        double r6970346 = pow(r6970345, r6970336);
        double r6970347 = r6970344 * r6970346;
        double r6970348 = r6970338 / r6970347;
        return r6970348;
}