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 r133766 = 1.0;
        double r133767 = s;
        double r133768 = -r133767;
        double r133769 = exp(r133768);
        double r133770 = r133766 + r133769;
        double r133771 = r133766 / r133770;
        double r133772 = c_p;
        double r133773 = pow(r133771, r133772);
        double r133774 = r133766 - r133771;
        double r133775 = c_n;
        double r133776 = pow(r133774, r133775);
        double r133777 = r133773 * r133776;
        double r133778 = t;
        double r133779 = -r133778;
        double r133780 = exp(r133779);
        double r133781 = r133766 + r133780;
        double r133782 = r133766 / r133781;
        double r133783 = pow(r133782, r133772);
        double r133784 = r133766 - r133782;
        double r133785 = pow(r133784, r133775);
        double r133786 = r133783 * r133785;
        double r133787 = r133777 / r133786;
        return r133787;
}