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