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 r137770 = 1.0;
        double r137771 = s;
        double r137772 = -r137771;
        double r137773 = exp(r137772);
        double r137774 = r137770 + r137773;
        double r137775 = r137770 / r137774;
        double r137776 = c_p;
        double r137777 = pow(r137775, r137776);
        double r137778 = r137770 - r137775;
        double r137779 = c_n;
        double r137780 = pow(r137778, r137779);
        double r137781 = r137777 * r137780;
        double r137782 = t;
        double r137783 = -r137782;
        double r137784 = exp(r137783);
        double r137785 = r137770 + r137784;
        double r137786 = r137770 / r137785;
        double r137787 = pow(r137786, r137776);
        double r137788 = r137770 - r137786;
        double r137789 = pow(r137788, r137779);
        double r137790 = r137787 * r137789;
        double r137791 = r137781 / r137790;
        return r137791;
}