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 r113449 = 1.0;
        double r113450 = s;
        double r113451 = -r113450;
        double r113452 = exp(r113451);
        double r113453 = r113449 + r113452;
        double r113454 = r113449 / r113453;
        double r113455 = c_p;
        double r113456 = pow(r113454, r113455);
        double r113457 = r113449 - r113454;
        double r113458 = c_n;
        double r113459 = pow(r113457, r113458);
        double r113460 = r113456 * r113459;
        double r113461 = t;
        double r113462 = -r113461;
        double r113463 = exp(r113462);
        double r113464 = r113449 + r113463;
        double r113465 = r113449 / r113464;
        double r113466 = pow(r113465, r113455);
        double r113467 = r113449 - r113465;
        double r113468 = pow(r113467, r113458);
        double r113469 = r113466 * r113468;
        double r113470 = r113460 / r113469;
        return r113470;
}