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 r133956 = 1.0;
        double r133957 = s;
        double r133958 = -r133957;
        double r133959 = exp(r133958);
        double r133960 = r133956 + r133959;
        double r133961 = r133956 / r133960;
        double r133962 = c_p;
        double r133963 = pow(r133961, r133962);
        double r133964 = r133956 - r133961;
        double r133965 = c_n;
        double r133966 = pow(r133964, r133965);
        double r133967 = r133963 * r133966;
        double r133968 = t;
        double r133969 = -r133968;
        double r133970 = exp(r133969);
        double r133971 = r133956 + r133970;
        double r133972 = r133956 / r133971;
        double r133973 = pow(r133972, r133962);
        double r133974 = r133956 - r133972;
        double r133975 = pow(r133974, r133965);
        double r133976 = r133973 * r133975;
        double r133977 = r133967 / r133976;
        return r133977;
}