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 r189310 = 1.0;
        double r189311 = s;
        double r189312 = -r189311;
        double r189313 = exp(r189312);
        double r189314 = r189310 + r189313;
        double r189315 = r189310 / r189314;
        double r189316 = c_p;
        double r189317 = pow(r189315, r189316);
        double r189318 = r189310 - r189315;
        double r189319 = c_n;
        double r189320 = pow(r189318, r189319);
        double r189321 = r189317 * r189320;
        double r189322 = t;
        double r189323 = -r189322;
        double r189324 = exp(r189323);
        double r189325 = r189310 + r189324;
        double r189326 = r189310 / r189325;
        double r189327 = pow(r189326, r189316);
        double r189328 = r189310 - r189326;
        double r189329 = pow(r189328, r189319);
        double r189330 = r189327 * r189329;
        double r189331 = r189321 / r189330;
        return r189331;
}