Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r82129216 = 1.0;
        double r82129217 = s;
        double r82129218 = -r82129217;
        double r82129219 = exp(r82129218);
        double r82129220 = r82129216 + r82129219;
        double r82129221 = r82129216 / r82129220;
        double r82129222 = c_p;
        double r82129223 = pow(r82129221, r82129222);
        double r82129224 = r82129216 - r82129221;
        double r82129225 = c_n;
        double r82129226 = pow(r82129224, r82129225);
        double r82129227 = r82129223 * r82129226;
        double r82129228 = t;
        double r82129229 = -r82129228;
        double r82129230 = exp(r82129229);
        double r82129231 = r82129216 + r82129230;
        double r82129232 = r82129216 / r82129231;
        double r82129233 = pow(r82129232, r82129222);
        double r82129234 = r82129216 - r82129232;
        double r82129235 = pow(r82129234, r82129225);
        double r82129236 = r82129233 * r82129235;
        double r82129237 = r82129227 / r82129236;
        return r82129237;
}