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 r153979 = 1.0;
        double r153980 = s;
        double r153981 = -r153980;
        double r153982 = exp(r153981);
        double r153983 = r153979 + r153982;
        double r153984 = r153979 / r153983;
        double r153985 = c_p;
        double r153986 = pow(r153984, r153985);
        double r153987 = r153979 - r153984;
        double r153988 = c_n;
        double r153989 = pow(r153987, r153988);
        double r153990 = r153986 * r153989;
        double r153991 = t;
        double r153992 = -r153991;
        double r153993 = exp(r153992);
        double r153994 = r153979 + r153993;
        double r153995 = r153979 / r153994;
        double r153996 = pow(r153995, r153985);
        double r153997 = r153979 - r153995;
        double r153998 = pow(r153997, r153988);
        double r153999 = r153996 * r153998;
        double r154000 = r153990 / r153999;
        return r154000;
}