Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r3481665 = 1.0;
        double r3481666 = s;
        double r3481667 = -r3481666;
        double r3481668 = exp(r3481667);
        double r3481669 = r3481665 + r3481668;
        double r3481670 = r3481665 / r3481669;
        double r3481671 = c_p;
        double r3481672 = pow(r3481670, r3481671);
        double r3481673 = r3481665 - r3481670;
        double r3481674 = c_n;
        double r3481675 = pow(r3481673, r3481674);
        double r3481676 = r3481672 * r3481675;
        double r3481677 = t;
        double r3481678 = -r3481677;
        double r3481679 = exp(r3481678);
        double r3481680 = r3481665 + r3481679;
        double r3481681 = r3481665 / r3481680;
        double r3481682 = pow(r3481681, r3481671);
        double r3481683 = r3481665 - r3481681;
        double r3481684 = pow(r3481683, r3481674);
        double r3481685 = r3481682 * r3481684;
        double r3481686 = r3481676 / r3481685;
        return r3481686;
}