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 r129631 = 1.0;
        double r129632 = s;
        double r129633 = -r129632;
        double r129634 = exp(r129633);
        double r129635 = r129631 + r129634;
        double r129636 = r129631 / r129635;
        double r129637 = c_p;
        double r129638 = pow(r129636, r129637);
        double r129639 = r129631 - r129636;
        double r129640 = c_n;
        double r129641 = pow(r129639, r129640);
        double r129642 = r129638 * r129641;
        double r129643 = t;
        double r129644 = -r129643;
        double r129645 = exp(r129644);
        double r129646 = r129631 + r129645;
        double r129647 = r129631 / r129646;
        double r129648 = pow(r129647, r129637);
        double r129649 = r129631 - r129647;
        double r129650 = pow(r129649, r129640);
        double r129651 = r129648 * r129650;
        double r129652 = r129642 / r129651;
        return r129652;
}