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 r663447 = 1.0;
        double r663448 = s;
        double r663449 = -r663448;
        double r663450 = exp(r663449);
        double r663451 = r663447 + r663450;
        double r663452 = r663447 / r663451;
        double r663453 = c_p;
        double r663454 = pow(r663452, r663453);
        double r663455 = r663447 - r663452;
        double r663456 = c_n;
        double r663457 = pow(r663455, r663456);
        double r663458 = r663454 * r663457;
        double r663459 = t;
        double r663460 = -r663459;
        double r663461 = exp(r663460);
        double r663462 = r663447 + r663461;
        double r663463 = r663447 / r663462;
        double r663464 = pow(r663463, r663453);
        double r663465 = r663447 - r663463;
        double r663466 = pow(r663465, r663456);
        double r663467 = r663464 * r663466;
        double r663468 = r663458 / r663467;
        return r663468;
}