Cannot sample enough valid points. (more)

\[0 \lt c_p \land 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 r6440604 = 1.0;
        double r6440605 = s;
        double r6440606 = -r6440605;
        double r6440607 = exp(r6440606);
        double r6440608 = r6440604 + r6440607;
        double r6440609 = r6440604 / r6440608;
        double r6440610 = c_p;
        double r6440611 = pow(r6440609, r6440610);
        double r6440612 = r6440604 - r6440609;
        double r6440613 = c_n;
        double r6440614 = pow(r6440612, r6440613);
        double r6440615 = r6440611 * r6440614;
        double r6440616 = t;
        double r6440617 = -r6440616;
        double r6440618 = exp(r6440617);
        double r6440619 = r6440604 + r6440618;
        double r6440620 = r6440604 / r6440619;
        double r6440621 = pow(r6440620, r6440610);
        double r6440622 = r6440604 - r6440620;
        double r6440623 = pow(r6440622, r6440613);
        double r6440624 = r6440621 * r6440623;
        double r6440625 = r6440615 / r6440624;
        return r6440625;
}