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 r2814097 = 1.0;
        double r2814098 = s;
        double r2814099 = -r2814098;
        double r2814100 = exp(r2814099);
        double r2814101 = r2814097 + r2814100;
        double r2814102 = r2814097 / r2814101;
        double r2814103 = c_p;
        double r2814104 = pow(r2814102, r2814103);
        double r2814105 = r2814097 - r2814102;
        double r2814106 = c_n;
        double r2814107 = pow(r2814105, r2814106);
        double r2814108 = r2814104 * r2814107;
        double r2814109 = t;
        double r2814110 = -r2814109;
        double r2814111 = exp(r2814110);
        double r2814112 = r2814097 + r2814111;
        double r2814113 = r2814097 / r2814112;
        double r2814114 = pow(r2814113, r2814103);
        double r2814115 = r2814097 - r2814113;
        double r2814116 = pow(r2814115, r2814106);
        double r2814117 = r2814114 * r2814116;
        double r2814118 = r2814108 / r2814117;
        return r2814118;
}