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 r189850 = 1.0;
        double r189851 = s;
        double r189852 = -r189851;
        double r189853 = exp(r189852);
        double r189854 = r189850 + r189853;
        double r189855 = r189850 / r189854;
        double r189856 = c_p;
        double r189857 = pow(r189855, r189856);
        double r189858 = r189850 - r189855;
        double r189859 = c_n;
        double r189860 = pow(r189858, r189859);
        double r189861 = r189857 * r189860;
        double r189862 = t;
        double r189863 = -r189862;
        double r189864 = exp(r189863);
        double r189865 = r189850 + r189864;
        double r189866 = r189850 / r189865;
        double r189867 = pow(r189866, r189856);
        double r189868 = r189850 - r189866;
        double r189869 = pow(r189868, r189859);
        double r189870 = r189867 * r189869;
        double r189871 = r189861 / r189870;
        return r189871;
}