Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r4230899 = 1.0;
        double r4230900 = s;
        double r4230901 = -r4230900;
        double r4230902 = exp(r4230901);
        double r4230903 = r4230899 + r4230902;
        double r4230904 = r4230899 / r4230903;
        double r4230905 = c_p;
        double r4230906 = pow(r4230904, r4230905);
        double r4230907 = r4230899 - r4230904;
        double r4230908 = c_n;
        double r4230909 = pow(r4230907, r4230908);
        double r4230910 = r4230906 * r4230909;
        double r4230911 = t;
        double r4230912 = -r4230911;
        double r4230913 = exp(r4230912);
        double r4230914 = r4230899 + r4230913;
        double r4230915 = r4230899 / r4230914;
        double r4230916 = pow(r4230915, r4230905);
        double r4230917 = r4230899 - r4230915;
        double r4230918 = pow(r4230917, r4230908);
        double r4230919 = r4230916 * r4230918;
        double r4230920 = r4230910 / r4230919;
        return r4230920;
}