Cannot sample enough valid points. (more)

\[0.0 \lt c_p \land 0.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 r11853253 = 1.0;
        double r11853254 = s;
        double r11853255 = -r11853254;
        double r11853256 = exp(r11853255);
        double r11853257 = r11853253 + r11853256;
        double r11853258 = r11853253 / r11853257;
        double r11853259 = c_p;
        double r11853260 = pow(r11853258, r11853259);
        double r11853261 = r11853253 - r11853258;
        double r11853262 = c_n;
        double r11853263 = pow(r11853261, r11853262);
        double r11853264 = r11853260 * r11853263;
        double r11853265 = t;
        double r11853266 = -r11853265;
        double r11853267 = exp(r11853266);
        double r11853268 = r11853253 + r11853267;
        double r11853269 = r11853253 / r11853268;
        double r11853270 = pow(r11853269, r11853259);
        double r11853271 = r11853253 - r11853269;
        double r11853272 = pow(r11853271, r11853262);
        double r11853273 = r11853270 * r11853272;
        double r11853274 = r11853264 / r11853273;
        return r11853274;
}