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 r170869 = 1.0;
        double r170870 = s;
        double r170871 = -r170870;
        double r170872 = exp(r170871);
        double r170873 = r170869 + r170872;
        double r170874 = r170869 / r170873;
        double r170875 = c_p;
        double r170876 = pow(r170874, r170875);
        double r170877 = r170869 - r170874;
        double r170878 = c_n;
        double r170879 = pow(r170877, r170878);
        double r170880 = r170876 * r170879;
        double r170881 = t;
        double r170882 = -r170881;
        double r170883 = exp(r170882);
        double r170884 = r170869 + r170883;
        double r170885 = r170869 / r170884;
        double r170886 = pow(r170885, r170875);
        double r170887 = r170869 - r170885;
        double r170888 = pow(r170887, r170878);
        double r170889 = r170886 * r170888;
        double r170890 = r170880 / r170889;
        return r170890;
}