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 r146680 = 1.0;
        double r146681 = s;
        double r146682 = -r146681;
        double r146683 = exp(r146682);
        double r146684 = r146680 + r146683;
        double r146685 = r146680 / r146684;
        double r146686 = c_p;
        double r146687 = pow(r146685, r146686);
        double r146688 = r146680 - r146685;
        double r146689 = c_n;
        double r146690 = pow(r146688, r146689);
        double r146691 = r146687 * r146690;
        double r146692 = t;
        double r146693 = -r146692;
        double r146694 = exp(r146693);
        double r146695 = r146680 + r146694;
        double r146696 = r146680 / r146695;
        double r146697 = pow(r146696, r146686);
        double r146698 = r146680 - r146696;
        double r146699 = pow(r146698, r146689);
        double r146700 = r146697 * r146699;
        double r146701 = r146691 / r146700;
        return r146701;
}