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 r673634 = 1.0;
        double r673635 = s;
        double r673636 = -r673635;
        double r673637 = exp(r673636);
        double r673638 = r673634 + r673637;
        double r673639 = r673634 / r673638;
        double r673640 = c_p;
        double r673641 = pow(r673639, r673640);
        double r673642 = r673634 - r673639;
        double r673643 = c_n;
        double r673644 = pow(r673642, r673643);
        double r673645 = r673641 * r673644;
        double r673646 = t;
        double r673647 = -r673646;
        double r673648 = exp(r673647);
        double r673649 = r673634 + r673648;
        double r673650 = r673634 / r673649;
        double r673651 = pow(r673650, r673640);
        double r673652 = r673634 - r673650;
        double r673653 = pow(r673652, r673643);
        double r673654 = r673651 * r673653;
        double r673655 = r673645 / r673654;
        return r673655;
}