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 r166182 = 1.0;
        double r166183 = s;
        double r166184 = -r166183;
        double r166185 = exp(r166184);
        double r166186 = r166182 + r166185;
        double r166187 = r166182 / r166186;
        double r166188 = c_p;
        double r166189 = pow(r166187, r166188);
        double r166190 = r166182 - r166187;
        double r166191 = c_n;
        double r166192 = pow(r166190, r166191);
        double r166193 = r166189 * r166192;
        double r166194 = t;
        double r166195 = -r166194;
        double r166196 = exp(r166195);
        double r166197 = r166182 + r166196;
        double r166198 = r166182 / r166197;
        double r166199 = pow(r166198, r166188);
        double r166200 = r166182 - r166198;
        double r166201 = pow(r166200, r166191);
        double r166202 = r166199 * r166201;
        double r166203 = r166193 / r166202;
        return r166203;
}