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 r153308 = 1.0;
        double r153309 = s;
        double r153310 = -r153309;
        double r153311 = exp(r153310);
        double r153312 = r153308 + r153311;
        double r153313 = r153308 / r153312;
        double r153314 = c_p;
        double r153315 = pow(r153313, r153314);
        double r153316 = r153308 - r153313;
        double r153317 = c_n;
        double r153318 = pow(r153316, r153317);
        double r153319 = r153315 * r153318;
        double r153320 = t;
        double r153321 = -r153320;
        double r153322 = exp(r153321);
        double r153323 = r153308 + r153322;
        double r153324 = r153308 / r153323;
        double r153325 = pow(r153324, r153314);
        double r153326 = r153308 - r153324;
        double r153327 = pow(r153326, r153317);
        double r153328 = r153325 * r153327;
        double r153329 = r153319 / r153328;
        return r153329;
}