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 r175558 = 1.0;
        double r175559 = s;
        double r175560 = -r175559;
        double r175561 = exp(r175560);
        double r175562 = r175558 + r175561;
        double r175563 = r175558 / r175562;
        double r175564 = c_p;
        double r175565 = pow(r175563, r175564);
        double r175566 = r175558 - r175563;
        double r175567 = c_n;
        double r175568 = pow(r175566, r175567);
        double r175569 = r175565 * r175568;
        double r175570 = t;
        double r175571 = -r175570;
        double r175572 = exp(r175571);
        double r175573 = r175558 + r175572;
        double r175574 = r175558 / r175573;
        double r175575 = pow(r175574, r175564);
        double r175576 = r175558 - r175574;
        double r175577 = pow(r175576, r175567);
        double r175578 = r175575 * r175577;
        double r175579 = r175569 / r175578;
        return r175579;
}