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 r150565 = 1.0;
        double r150566 = s;
        double r150567 = -r150566;
        double r150568 = exp(r150567);
        double r150569 = r150565 + r150568;
        double r150570 = r150565 / r150569;
        double r150571 = c_p;
        double r150572 = pow(r150570, r150571);
        double r150573 = r150565 - r150570;
        double r150574 = c_n;
        double r150575 = pow(r150573, r150574);
        double r150576 = r150572 * r150575;
        double r150577 = t;
        double r150578 = -r150577;
        double r150579 = exp(r150578);
        double r150580 = r150565 + r150579;
        double r150581 = r150565 / r150580;
        double r150582 = pow(r150581, r150571);
        double r150583 = r150565 - r150581;
        double r150584 = pow(r150583, r150574);
        double r150585 = r150582 * r150584;
        double r150586 = r150576 / r150585;
        return r150586;
}