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 r181141 = 1.0;
        double r181142 = s;
        double r181143 = -r181142;
        double r181144 = exp(r181143);
        double r181145 = r181141 + r181144;
        double r181146 = r181141 / r181145;
        double r181147 = c_p;
        double r181148 = pow(r181146, r181147);
        double r181149 = r181141 - r181146;
        double r181150 = c_n;
        double r181151 = pow(r181149, r181150);
        double r181152 = r181148 * r181151;
        double r181153 = t;
        double r181154 = -r181153;
        double r181155 = exp(r181154);
        double r181156 = r181141 + r181155;
        double r181157 = r181141 / r181156;
        double r181158 = pow(r181157, r181147);
        double r181159 = r181141 - r181157;
        double r181160 = pow(r181159, r181150);
        double r181161 = r181158 * r181160;
        double r181162 = r181152 / r181161;
        return r181162;
}