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 r78927 = 1.0;
        double r78928 = s;
        double r78929 = -r78928;
        double r78930 = exp(r78929);
        double r78931 = r78927 + r78930;
        double r78932 = r78927 / r78931;
        double r78933 = c_p;
        double r78934 = pow(r78932, r78933);
        double r78935 = r78927 - r78932;
        double r78936 = c_n;
        double r78937 = pow(r78935, r78936);
        double r78938 = r78934 * r78937;
        double r78939 = t;
        double r78940 = -r78939;
        double r78941 = exp(r78940);
        double r78942 = r78927 + r78941;
        double r78943 = r78927 / r78942;
        double r78944 = pow(r78943, r78933);
        double r78945 = r78927 - r78943;
        double r78946 = pow(r78945, r78936);
        double r78947 = r78944 * r78946;
        double r78948 = r78938 / r78947;
        return r78948;
}