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 r196491 = 1.0;
        double r196492 = s;
        double r196493 = -r196492;
        double r196494 = exp(r196493);
        double r196495 = r196491 + r196494;
        double r196496 = r196491 / r196495;
        double r196497 = c_p;
        double r196498 = pow(r196496, r196497);
        double r196499 = r196491 - r196496;
        double r196500 = c_n;
        double r196501 = pow(r196499, r196500);
        double r196502 = r196498 * r196501;
        double r196503 = t;
        double r196504 = -r196503;
        double r196505 = exp(r196504);
        double r196506 = r196491 + r196505;
        double r196507 = r196491 / r196506;
        double r196508 = pow(r196507, r196497);
        double r196509 = r196491 - r196507;
        double r196510 = pow(r196509, r196500);
        double r196511 = r196508 * r196510;
        double r196512 = r196502 / r196511;
        return r196512;
}