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 r4552175 = 1.0;
        double r4552176 = s;
        double r4552177 = -r4552176;
        double r4552178 = exp(r4552177);
        double r4552179 = r4552175 + r4552178;
        double r4552180 = r4552175 / r4552179;
        double r4552181 = c_p;
        double r4552182 = pow(r4552180, r4552181);
        double r4552183 = r4552175 - r4552180;
        double r4552184 = c_n;
        double r4552185 = pow(r4552183, r4552184);
        double r4552186 = r4552182 * r4552185;
        double r4552187 = t;
        double r4552188 = -r4552187;
        double r4552189 = exp(r4552188);
        double r4552190 = r4552175 + r4552189;
        double r4552191 = r4552175 / r4552190;
        double r4552192 = pow(r4552191, r4552181);
        double r4552193 = r4552175 - r4552191;
        double r4552194 = pow(r4552193, r4552184);
        double r4552195 = r4552192 * r4552194;
        double r4552196 = r4552186 / r4552195;
        return r4552196;
}