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 r111768 = 1.0;
        double r111769 = s;
        double r111770 = -r111769;
        double r111771 = exp(r111770);
        double r111772 = r111768 + r111771;
        double r111773 = r111768 / r111772;
        double r111774 = c_p;
        double r111775 = pow(r111773, r111774);
        double r111776 = r111768 - r111773;
        double r111777 = c_n;
        double r111778 = pow(r111776, r111777);
        double r111779 = r111775 * r111778;
        double r111780 = t;
        double r111781 = -r111780;
        double r111782 = exp(r111781);
        double r111783 = r111768 + r111782;
        double r111784 = r111768 / r111783;
        double r111785 = pow(r111784, r111774);
        double r111786 = r111768 - r111784;
        double r111787 = pow(r111786, r111777);
        double r111788 = r111785 * r111787;
        double r111789 = r111779 / r111788;
        return r111789;
}