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 r151259 = 1.0;
        double r151260 = s;
        double r151261 = -r151260;
        double r151262 = exp(r151261);
        double r151263 = r151259 + r151262;
        double r151264 = r151259 / r151263;
        double r151265 = c_p;
        double r151266 = pow(r151264, r151265);
        double r151267 = r151259 - r151264;
        double r151268 = c_n;
        double r151269 = pow(r151267, r151268);
        double r151270 = r151266 * r151269;
        double r151271 = t;
        double r151272 = -r151271;
        double r151273 = exp(r151272);
        double r151274 = r151259 + r151273;
        double r151275 = r151259 / r151274;
        double r151276 = pow(r151275, r151265);
        double r151277 = r151259 - r151275;
        double r151278 = pow(r151277, r151268);
        double r151279 = r151276 * r151278;
        double r151280 = r151270 / r151279;
        return r151280;
}