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 r139179 = 1.0;
        double r139180 = s;
        double r139181 = -r139180;
        double r139182 = exp(r139181);
        double r139183 = r139179 + r139182;
        double r139184 = r139179 / r139183;
        double r139185 = c_p;
        double r139186 = pow(r139184, r139185);
        double r139187 = r139179 - r139184;
        double r139188 = c_n;
        double r139189 = pow(r139187, r139188);
        double r139190 = r139186 * r139189;
        double r139191 = t;
        double r139192 = -r139191;
        double r139193 = exp(r139192);
        double r139194 = r139179 + r139193;
        double r139195 = r139179 / r139194;
        double r139196 = pow(r139195, r139185);
        double r139197 = r139179 - r139195;
        double r139198 = pow(r139197, r139188);
        double r139199 = r139196 * r139198;
        double r139200 = r139190 / r139199;
        return r139200;
}