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 r181125 = 1.0;
        double r181126 = s;
        double r181127 = -r181126;
        double r181128 = exp(r181127);
        double r181129 = r181125 + r181128;
        double r181130 = r181125 / r181129;
        double r181131 = c_p;
        double r181132 = pow(r181130, r181131);
        double r181133 = r181125 - r181130;
        double r181134 = c_n;
        double r181135 = pow(r181133, r181134);
        double r181136 = r181132 * r181135;
        double r181137 = t;
        double r181138 = -r181137;
        double r181139 = exp(r181138);
        double r181140 = r181125 + r181139;
        double r181141 = r181125 / r181140;
        double r181142 = pow(r181141, r181131);
        double r181143 = r181125 - r181141;
        double r181144 = pow(r181143, r181134);
        double r181145 = r181142 * r181144;
        double r181146 = r181136 / r181145;
        return r181146;
}