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 r6421922 = 1.0;
        double r6421923 = s;
        double r6421924 = -r6421923;
        double r6421925 = exp(r6421924);
        double r6421926 = r6421922 + r6421925;
        double r6421927 = r6421922 / r6421926;
        double r6421928 = c_p;
        double r6421929 = pow(r6421927, r6421928);
        double r6421930 = r6421922 - r6421927;
        double r6421931 = c_n;
        double r6421932 = pow(r6421930, r6421931);
        double r6421933 = r6421929 * r6421932;
        double r6421934 = t;
        double r6421935 = -r6421934;
        double r6421936 = exp(r6421935);
        double r6421937 = r6421922 + r6421936;
        double r6421938 = r6421922 / r6421937;
        double r6421939 = pow(r6421938, r6421928);
        double r6421940 = r6421922 - r6421938;
        double r6421941 = pow(r6421940, r6421931);
        double r6421942 = r6421939 * r6421941;
        double r6421943 = r6421933 / r6421942;
        return r6421943;
}