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 r125362 = 1.0;
        double r125363 = s;
        double r125364 = -r125363;
        double r125365 = exp(r125364);
        double r125366 = r125362 + r125365;
        double r125367 = r125362 / r125366;
        double r125368 = c_p;
        double r125369 = pow(r125367, r125368);
        double r125370 = r125362 - r125367;
        double r125371 = c_n;
        double r125372 = pow(r125370, r125371);
        double r125373 = r125369 * r125372;
        double r125374 = t;
        double r125375 = -r125374;
        double r125376 = exp(r125375);
        double r125377 = r125362 + r125376;
        double r125378 = r125362 / r125377;
        double r125379 = pow(r125378, r125368);
        double r125380 = r125362 - r125378;
        double r125381 = pow(r125380, r125371);
        double r125382 = r125379 * r125381;
        double r125383 = r125373 / r125382;
        return r125383;
}