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 r6814260 = 1.0;
        double r6814261 = s;
        double r6814262 = -r6814261;
        double r6814263 = exp(r6814262);
        double r6814264 = r6814260 + r6814263;
        double r6814265 = r6814260 / r6814264;
        double r6814266 = c_p;
        double r6814267 = pow(r6814265, r6814266);
        double r6814268 = r6814260 - r6814265;
        double r6814269 = c_n;
        double r6814270 = pow(r6814268, r6814269);
        double r6814271 = r6814267 * r6814270;
        double r6814272 = t;
        double r6814273 = -r6814272;
        double r6814274 = exp(r6814273);
        double r6814275 = r6814260 + r6814274;
        double r6814276 = r6814260 / r6814275;
        double r6814277 = pow(r6814276, r6814266);
        double r6814278 = r6814260 - r6814276;
        double r6814279 = pow(r6814278, r6814269);
        double r6814280 = r6814277 * r6814279;
        double r6814281 = r6814271 / r6814280;
        return r6814281;
}