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 r7511861 = 1.0;
        double r7511862 = s;
        double r7511863 = -r7511862;
        double r7511864 = exp(r7511863);
        double r7511865 = r7511861 + r7511864;
        double r7511866 = r7511861 / r7511865;
        double r7511867 = c_p;
        double r7511868 = pow(r7511866, r7511867);
        double r7511869 = r7511861 - r7511866;
        double r7511870 = c_n;
        double r7511871 = pow(r7511869, r7511870);
        double r7511872 = r7511868 * r7511871;
        double r7511873 = t;
        double r7511874 = -r7511873;
        double r7511875 = exp(r7511874);
        double r7511876 = r7511861 + r7511875;
        double r7511877 = r7511861 / r7511876;
        double r7511878 = pow(r7511877, r7511867);
        double r7511879 = r7511861 - r7511877;
        double r7511880 = pow(r7511879, r7511870);
        double r7511881 = r7511878 * r7511880;
        double r7511882 = r7511872 / r7511881;
        return r7511882;
}