Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r7526214 = 1.0;
        double r7526215 = s;
        double r7526216 = -r7526215;
        double r7526217 = exp(r7526216);
        double r7526218 = r7526214 + r7526217;
        double r7526219 = r7526214 / r7526218;
        double r7526220 = c_p;
        double r7526221 = pow(r7526219, r7526220);
        double r7526222 = r7526214 - r7526219;
        double r7526223 = c_n;
        double r7526224 = pow(r7526222, r7526223);
        double r7526225 = r7526221 * r7526224;
        double r7526226 = t;
        double r7526227 = -r7526226;
        double r7526228 = exp(r7526227);
        double r7526229 = r7526214 + r7526228;
        double r7526230 = r7526214 / r7526229;
        double r7526231 = pow(r7526230, r7526220);
        double r7526232 = r7526214 - r7526230;
        double r7526233 = pow(r7526232, r7526223);
        double r7526234 = r7526231 * r7526233;
        double r7526235 = r7526225 / r7526234;
        return r7526235;
}