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 r150452 = 1.0;
        double r150453 = s;
        double r150454 = -r150453;
        double r150455 = exp(r150454);
        double r150456 = r150452 + r150455;
        double r150457 = r150452 / r150456;
        double r150458 = c_p;
        double r150459 = pow(r150457, r150458);
        double r150460 = r150452 - r150457;
        double r150461 = c_n;
        double r150462 = pow(r150460, r150461);
        double r150463 = r150459 * r150462;
        double r150464 = t;
        double r150465 = -r150464;
        double r150466 = exp(r150465);
        double r150467 = r150452 + r150466;
        double r150468 = r150452 / r150467;
        double r150469 = pow(r150468, r150458);
        double r150470 = r150452 - r150468;
        double r150471 = pow(r150470, r150461);
        double r150472 = r150469 * r150471;
        double r150473 = r150463 / r150472;
        return r150473;
}