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 r128686 = 1.0;
        double r128687 = s;
        double r128688 = -r128687;
        double r128689 = exp(r128688);
        double r128690 = r128686 + r128689;
        double r128691 = r128686 / r128690;
        double r128692 = c_p;
        double r128693 = pow(r128691, r128692);
        double r128694 = r128686 - r128691;
        double r128695 = c_n;
        double r128696 = pow(r128694, r128695);
        double r128697 = r128693 * r128696;
        double r128698 = t;
        double r128699 = -r128698;
        double r128700 = exp(r128699);
        double r128701 = r128686 + r128700;
        double r128702 = r128686 / r128701;
        double r128703 = pow(r128702, r128692);
        double r128704 = r128686 - r128702;
        double r128705 = pow(r128704, r128695);
        double r128706 = r128703 * r128705;
        double r128707 = r128697 / r128706;
        return r128707;
}