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 r6348607 = 1.0;
        double r6348608 = s;
        double r6348609 = -r6348608;
        double r6348610 = exp(r6348609);
        double r6348611 = r6348607 + r6348610;
        double r6348612 = r6348607 / r6348611;
        double r6348613 = c_p;
        double r6348614 = pow(r6348612, r6348613);
        double r6348615 = r6348607 - r6348612;
        double r6348616 = c_n;
        double r6348617 = pow(r6348615, r6348616);
        double r6348618 = r6348614 * r6348617;
        double r6348619 = t;
        double r6348620 = -r6348619;
        double r6348621 = exp(r6348620);
        double r6348622 = r6348607 + r6348621;
        double r6348623 = r6348607 / r6348622;
        double r6348624 = pow(r6348623, r6348613);
        double r6348625 = r6348607 - r6348623;
        double r6348626 = pow(r6348625, r6348616);
        double r6348627 = r6348624 * r6348626;
        double r6348628 = r6348618 / r6348627;
        return r6348628;
}