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 r2600395 = 1.0;
        double r2600396 = s;
        double r2600397 = -r2600396;
        double r2600398 = exp(r2600397);
        double r2600399 = r2600395 + r2600398;
        double r2600400 = r2600395 / r2600399;
        double r2600401 = c_p;
        double r2600402 = pow(r2600400, r2600401);
        double r2600403 = r2600395 - r2600400;
        double r2600404 = c_n;
        double r2600405 = pow(r2600403, r2600404);
        double r2600406 = r2600402 * r2600405;
        double r2600407 = t;
        double r2600408 = -r2600407;
        double r2600409 = exp(r2600408);
        double r2600410 = r2600395 + r2600409;
        double r2600411 = r2600395 / r2600410;
        double r2600412 = pow(r2600411, r2600401);
        double r2600413 = r2600395 - r2600411;
        double r2600414 = pow(r2600413, r2600404);
        double r2600415 = r2600412 * r2600414;
        double r2600416 = r2600406 / r2600415;
        return r2600416;
}