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 r173191 = 1.0;
        double r173192 = s;
        double r173193 = -r173192;
        double r173194 = exp(r173193);
        double r173195 = r173191 + r173194;
        double r173196 = r173191 / r173195;
        double r173197 = c_p;
        double r173198 = pow(r173196, r173197);
        double r173199 = r173191 - r173196;
        double r173200 = c_n;
        double r173201 = pow(r173199, r173200);
        double r173202 = r173198 * r173201;
        double r173203 = t;
        double r173204 = -r173203;
        double r173205 = exp(r173204);
        double r173206 = r173191 + r173205;
        double r173207 = r173191 / r173206;
        double r173208 = pow(r173207, r173197);
        double r173209 = r173191 - r173207;
        double r173210 = pow(r173209, r173200);
        double r173211 = r173208 * r173210;
        double r173212 = r173202 / r173211;
        return r173212;
}