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 r3024171 = 1.0;
        double r3024172 = s;
        double r3024173 = -r3024172;
        double r3024174 = exp(r3024173);
        double r3024175 = r3024171 + r3024174;
        double r3024176 = r3024171 / r3024175;
        double r3024177 = c_p;
        double r3024178 = pow(r3024176, r3024177);
        double r3024179 = r3024171 - r3024176;
        double r3024180 = c_n;
        double r3024181 = pow(r3024179, r3024180);
        double r3024182 = r3024178 * r3024181;
        double r3024183 = t;
        double r3024184 = -r3024183;
        double r3024185 = exp(r3024184);
        double r3024186 = r3024171 + r3024185;
        double r3024187 = r3024171 / r3024186;
        double r3024188 = pow(r3024187, r3024177);
        double r3024189 = r3024171 - r3024187;
        double r3024190 = pow(r3024189, r3024180);
        double r3024191 = r3024188 * r3024190;
        double r3024192 = r3024182 / r3024191;
        return r3024192;
}