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 r7907728 = 1.0;
        double r7907729 = s;
        double r7907730 = -r7907729;
        double r7907731 = exp(r7907730);
        double r7907732 = r7907728 + r7907731;
        double r7907733 = r7907728 / r7907732;
        double r7907734 = c_p;
        double r7907735 = pow(r7907733, r7907734);
        double r7907736 = r7907728 - r7907733;
        double r7907737 = c_n;
        double r7907738 = pow(r7907736, r7907737);
        double r7907739 = r7907735 * r7907738;
        double r7907740 = t;
        double r7907741 = -r7907740;
        double r7907742 = exp(r7907741);
        double r7907743 = r7907728 + r7907742;
        double r7907744 = r7907728 / r7907743;
        double r7907745 = pow(r7907744, r7907734);
        double r7907746 = r7907728 - r7907744;
        double r7907747 = pow(r7907746, r7907737);
        double r7907748 = r7907745 * r7907747;
        double r7907749 = r7907739 / r7907748;
        return r7907749;
}