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 r124431 = 1.0;
        double r124432 = s;
        double r124433 = -r124432;
        double r124434 = exp(r124433);
        double r124435 = r124431 + r124434;
        double r124436 = r124431 / r124435;
        double r124437 = c_p;
        double r124438 = pow(r124436, r124437);
        double r124439 = r124431 - r124436;
        double r124440 = c_n;
        double r124441 = pow(r124439, r124440);
        double r124442 = r124438 * r124441;
        double r124443 = t;
        double r124444 = -r124443;
        double r124445 = exp(r124444);
        double r124446 = r124431 + r124445;
        double r124447 = r124431 / r124446;
        double r124448 = pow(r124447, r124437);
        double r124449 = r124431 - r124447;
        double r124450 = pow(r124449, r124440);
        double r124451 = r124448 * r124450;
        double r124452 = r124442 / r124451;
        return r124452;
}