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 r154247 = 1.0;
        double r154248 = s;
        double r154249 = -r154248;
        double r154250 = exp(r154249);
        double r154251 = r154247 + r154250;
        double r154252 = r154247 / r154251;
        double r154253 = c_p;
        double r154254 = pow(r154252, r154253);
        double r154255 = r154247 - r154252;
        double r154256 = c_n;
        double r154257 = pow(r154255, r154256);
        double r154258 = r154254 * r154257;
        double r154259 = t;
        double r154260 = -r154259;
        double r154261 = exp(r154260);
        double r154262 = r154247 + r154261;
        double r154263 = r154247 / r154262;
        double r154264 = pow(r154263, r154253);
        double r154265 = r154247 - r154263;
        double r154266 = pow(r154265, r154256);
        double r154267 = r154264 * r154266;
        double r154268 = r154258 / r154267;
        return r154268;
}