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 r6737225 = 1.0;
        double r6737226 = s;
        double r6737227 = -r6737226;
        double r6737228 = exp(r6737227);
        double r6737229 = r6737225 + r6737228;
        double r6737230 = r6737225 / r6737229;
        double r6737231 = c_p;
        double r6737232 = pow(r6737230, r6737231);
        double r6737233 = r6737225 - r6737230;
        double r6737234 = c_n;
        double r6737235 = pow(r6737233, r6737234);
        double r6737236 = r6737232 * r6737235;
        double r6737237 = t;
        double r6737238 = -r6737237;
        double r6737239 = exp(r6737238);
        double r6737240 = r6737225 + r6737239;
        double r6737241 = r6737225 / r6737240;
        double r6737242 = pow(r6737241, r6737231);
        double r6737243 = r6737225 - r6737241;
        double r6737244 = pow(r6737243, r6737234);
        double r6737245 = r6737242 * r6737244;
        double r6737246 = r6737236 / r6737245;
        return r6737246;
}