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 r130195 = 1.0;
        double r130196 = s;
        double r130197 = -r130196;
        double r130198 = exp(r130197);
        double r130199 = r130195 + r130198;
        double r130200 = r130195 / r130199;
        double r130201 = c_p;
        double r130202 = pow(r130200, r130201);
        double r130203 = r130195 - r130200;
        double r130204 = c_n;
        double r130205 = pow(r130203, r130204);
        double r130206 = r130202 * r130205;
        double r130207 = t;
        double r130208 = -r130207;
        double r130209 = exp(r130208);
        double r130210 = r130195 + r130209;
        double r130211 = r130195 / r130210;
        double r130212 = pow(r130211, r130201);
        double r130213 = r130195 - r130211;
        double r130214 = pow(r130213, r130204);
        double r130215 = r130212 * r130214;
        double r130216 = r130206 / r130215;
        return r130216;
}