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 r111864 = 1.0;
        double r111865 = s;
        double r111866 = -r111865;
        double r111867 = exp(r111866);
        double r111868 = r111864 + r111867;
        double r111869 = r111864 / r111868;
        double r111870 = c_p;
        double r111871 = pow(r111869, r111870);
        double r111872 = r111864 - r111869;
        double r111873 = c_n;
        double r111874 = pow(r111872, r111873);
        double r111875 = r111871 * r111874;
        double r111876 = t;
        double r111877 = -r111876;
        double r111878 = exp(r111877);
        double r111879 = r111864 + r111878;
        double r111880 = r111864 / r111879;
        double r111881 = pow(r111880, r111870);
        double r111882 = r111864 - r111880;
        double r111883 = pow(r111882, r111873);
        double r111884 = r111881 * r111883;
        double r111885 = r111875 / r111884;
        return r111885;
}