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 r171813 = 1.0;
        double r171814 = s;
        double r171815 = -r171814;
        double r171816 = exp(r171815);
        double r171817 = r171813 + r171816;
        double r171818 = r171813 / r171817;
        double r171819 = c_p;
        double r171820 = pow(r171818, r171819);
        double r171821 = r171813 - r171818;
        double r171822 = c_n;
        double r171823 = pow(r171821, r171822);
        double r171824 = r171820 * r171823;
        double r171825 = t;
        double r171826 = -r171825;
        double r171827 = exp(r171826);
        double r171828 = r171813 + r171827;
        double r171829 = r171813 / r171828;
        double r171830 = pow(r171829, r171819);
        double r171831 = r171813 - r171829;
        double r171832 = pow(r171831, r171822);
        double r171833 = r171830 * r171832;
        double r171834 = r171824 / r171833;
        return r171834;
}