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 r225177 = 1.0;
        double r225178 = s;
        double r225179 = -r225178;
        double r225180 = exp(r225179);
        double r225181 = r225177 + r225180;
        double r225182 = r225177 / r225181;
        double r225183 = c_p;
        double r225184 = pow(r225182, r225183);
        double r225185 = r225177 - r225182;
        double r225186 = c_n;
        double r225187 = pow(r225185, r225186);
        double r225188 = r225184 * r225187;
        double r225189 = t;
        double r225190 = -r225189;
        double r225191 = exp(r225190);
        double r225192 = r225177 + r225191;
        double r225193 = r225177 / r225192;
        double r225194 = pow(r225193, r225183);
        double r225195 = r225177 - r225193;
        double r225196 = pow(r225195, r225186);
        double r225197 = r225194 * r225196;
        double r225198 = r225188 / r225197;
        return r225198;
}