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 r127692 = 1.0;
        double r127693 = s;
        double r127694 = -r127693;
        double r127695 = exp(r127694);
        double r127696 = r127692 + r127695;
        double r127697 = r127692 / r127696;
        double r127698 = c_p;
        double r127699 = pow(r127697, r127698);
        double r127700 = r127692 - r127697;
        double r127701 = c_n;
        double r127702 = pow(r127700, r127701);
        double r127703 = r127699 * r127702;
        double r127704 = t;
        double r127705 = -r127704;
        double r127706 = exp(r127705);
        double r127707 = r127692 + r127706;
        double r127708 = r127692 / r127707;
        double r127709 = pow(r127708, r127698);
        double r127710 = r127692 - r127708;
        double r127711 = pow(r127710, r127701);
        double r127712 = r127709 * r127711;
        double r127713 = r127703 / r127712;
        return r127713;
}