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 r2122799 = 1.0;
        double r2122800 = s;
        double r2122801 = -r2122800;
        double r2122802 = exp(r2122801);
        double r2122803 = r2122799 + r2122802;
        double r2122804 = r2122799 / r2122803;
        double r2122805 = c_p;
        double r2122806 = pow(r2122804, r2122805);
        double r2122807 = r2122799 - r2122804;
        double r2122808 = c_n;
        double r2122809 = pow(r2122807, r2122808);
        double r2122810 = r2122806 * r2122809;
        double r2122811 = t;
        double r2122812 = -r2122811;
        double r2122813 = exp(r2122812);
        double r2122814 = r2122799 + r2122813;
        double r2122815 = r2122799 / r2122814;
        double r2122816 = pow(r2122815, r2122805);
        double r2122817 = r2122799 - r2122815;
        double r2122818 = pow(r2122817, r2122808);
        double r2122819 = r2122816 * r2122818;
        double r2122820 = r2122810 / r2122819;
        return r2122820;
}