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 r138939 = 1.0;
        double r138940 = s;
        double r138941 = -r138940;
        double r138942 = exp(r138941);
        double r138943 = r138939 + r138942;
        double r138944 = r138939 / r138943;
        double r138945 = c_p;
        double r138946 = pow(r138944, r138945);
        double r138947 = r138939 - r138944;
        double r138948 = c_n;
        double r138949 = pow(r138947, r138948);
        double r138950 = r138946 * r138949;
        double r138951 = t;
        double r138952 = -r138951;
        double r138953 = exp(r138952);
        double r138954 = r138939 + r138953;
        double r138955 = r138939 / r138954;
        double r138956 = pow(r138955, r138945);
        double r138957 = r138939 - r138955;
        double r138958 = pow(r138957, r138948);
        double r138959 = r138956 * r138958;
        double r138960 = r138950 / r138959;
        return r138960;
}