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 r143642 = 1.0;
        double r143643 = s;
        double r143644 = -r143643;
        double r143645 = exp(r143644);
        double r143646 = r143642 + r143645;
        double r143647 = r143642 / r143646;
        double r143648 = c_p;
        double r143649 = pow(r143647, r143648);
        double r143650 = r143642 - r143647;
        double r143651 = c_n;
        double r143652 = pow(r143650, r143651);
        double r143653 = r143649 * r143652;
        double r143654 = t;
        double r143655 = -r143654;
        double r143656 = exp(r143655);
        double r143657 = r143642 + r143656;
        double r143658 = r143642 / r143657;
        double r143659 = pow(r143658, r143648);
        double r143660 = r143642 - r143658;
        double r143661 = pow(r143660, r143651);
        double r143662 = r143659 * r143661;
        double r143663 = r143653 / r143662;
        return r143663;
}