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 r3009947 = 1.0;
        double r3009948 = s;
        double r3009949 = -r3009948;
        double r3009950 = exp(r3009949);
        double r3009951 = r3009947 + r3009950;
        double r3009952 = r3009947 / r3009951;
        double r3009953 = c_p;
        double r3009954 = pow(r3009952, r3009953);
        double r3009955 = r3009947 - r3009952;
        double r3009956 = c_n;
        double r3009957 = pow(r3009955, r3009956);
        double r3009958 = r3009954 * r3009957;
        double r3009959 = t;
        double r3009960 = -r3009959;
        double r3009961 = exp(r3009960);
        double r3009962 = r3009947 + r3009961;
        double r3009963 = r3009947 / r3009962;
        double r3009964 = pow(r3009963, r3009953);
        double r3009965 = r3009947 - r3009963;
        double r3009966 = pow(r3009965, r3009956);
        double r3009967 = r3009964 * r3009966;
        double r3009968 = r3009958 / r3009967;
        return r3009968;
}