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 r248919 = 1.0;
        double r248920 = s;
        double r248921 = -r248920;
        double r248922 = exp(r248921);
        double r248923 = r248919 + r248922;
        double r248924 = r248919 / r248923;
        double r248925 = c_p;
        double r248926 = pow(r248924, r248925);
        double r248927 = r248919 - r248924;
        double r248928 = c_n;
        double r248929 = pow(r248927, r248928);
        double r248930 = r248926 * r248929;
        double r248931 = t;
        double r248932 = -r248931;
        double r248933 = exp(r248932);
        double r248934 = r248919 + r248933;
        double r248935 = r248919 / r248934;
        double r248936 = pow(r248935, r248925);
        double r248937 = r248919 - r248935;
        double r248938 = pow(r248937, r248928);
        double r248939 = r248936 * r248938;
        double r248940 = r248930 / r248939;
        return r248940;
}