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 r164173 = 1.0;
        double r164174 = s;
        double r164175 = -r164174;
        double r164176 = exp(r164175);
        double r164177 = r164173 + r164176;
        double r164178 = r164173 / r164177;
        double r164179 = c_p;
        double r164180 = pow(r164178, r164179);
        double r164181 = r164173 - r164178;
        double r164182 = c_n;
        double r164183 = pow(r164181, r164182);
        double r164184 = r164180 * r164183;
        double r164185 = t;
        double r164186 = -r164185;
        double r164187 = exp(r164186);
        double r164188 = r164173 + r164187;
        double r164189 = r164173 / r164188;
        double r164190 = pow(r164189, r164179);
        double r164191 = r164173 - r164189;
        double r164192 = pow(r164191, r164182);
        double r164193 = r164190 * r164192;
        double r164194 = r164184 / r164193;
        return r164194;
}