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 r142980 = 1.0;
        double r142981 = s;
        double r142982 = -r142981;
        double r142983 = exp(r142982);
        double r142984 = r142980 + r142983;
        double r142985 = r142980 / r142984;
        double r142986 = c_p;
        double r142987 = pow(r142985, r142986);
        double r142988 = r142980 - r142985;
        double r142989 = c_n;
        double r142990 = pow(r142988, r142989);
        double r142991 = r142987 * r142990;
        double r142992 = t;
        double r142993 = -r142992;
        double r142994 = exp(r142993);
        double r142995 = r142980 + r142994;
        double r142996 = r142980 / r142995;
        double r142997 = pow(r142996, r142986);
        double r142998 = r142980 - r142996;
        double r142999 = pow(r142998, r142989);
        double r143000 = r142997 * r142999;
        double r143001 = r142991 / r143000;
        return r143001;
}