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 r19608204 = 1.0;
        double r19608205 = s;
        double r19608206 = -r19608205;
        double r19608207 = exp(r19608206);
        double r19608208 = r19608204 + r19608207;
        double r19608209 = r19608204 / r19608208;
        double r19608210 = c_p;
        double r19608211 = pow(r19608209, r19608210);
        double r19608212 = r19608204 - r19608209;
        double r19608213 = c_n;
        double r19608214 = pow(r19608212, r19608213);
        double r19608215 = r19608211 * r19608214;
        double r19608216 = t;
        double r19608217 = -r19608216;
        double r19608218 = exp(r19608217);
        double r19608219 = r19608204 + r19608218;
        double r19608220 = r19608204 / r19608219;
        double r19608221 = pow(r19608220, r19608210);
        double r19608222 = r19608204 - r19608220;
        double r19608223 = pow(r19608222, r19608213);
        double r19608224 = r19608221 * r19608223;
        double r19608225 = r19608215 / r19608224;
        return r19608225;
}