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 r147553 = 1.0;
        double r147554 = s;
        double r147555 = -r147554;
        double r147556 = exp(r147555);
        double r147557 = r147553 + r147556;
        double r147558 = r147553 / r147557;
        double r147559 = c_p;
        double r147560 = pow(r147558, r147559);
        double r147561 = r147553 - r147558;
        double r147562 = c_n;
        double r147563 = pow(r147561, r147562);
        double r147564 = r147560 * r147563;
        double r147565 = t;
        double r147566 = -r147565;
        double r147567 = exp(r147566);
        double r147568 = r147553 + r147567;
        double r147569 = r147553 / r147568;
        double r147570 = pow(r147569, r147559);
        double r147571 = r147553 - r147569;
        double r147572 = pow(r147571, r147562);
        double r147573 = r147570 * r147572;
        double r147574 = r147564 / r147573;
        return r147574;
}