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 r53675199 = 1.0;
        double r53675200 = s;
        double r53675201 = -r53675200;
        double r53675202 = exp(r53675201);
        double r53675203 = r53675199 + r53675202;
        double r53675204 = r53675199 / r53675203;
        double r53675205 = c_p;
        double r53675206 = pow(r53675204, r53675205);
        double r53675207 = r53675199 - r53675204;
        double r53675208 = c_n;
        double r53675209 = pow(r53675207, r53675208);
        double r53675210 = r53675206 * r53675209;
        double r53675211 = t;
        double r53675212 = -r53675211;
        double r53675213 = exp(r53675212);
        double r53675214 = r53675199 + r53675213;
        double r53675215 = r53675199 / r53675214;
        double r53675216 = pow(r53675215, r53675205);
        double r53675217 = r53675199 - r53675215;
        double r53675218 = pow(r53675217, r53675208);
        double r53675219 = r53675216 * r53675218;
        double r53675220 = r53675210 / r53675219;
        return r53675220;
}