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 r149055 = 1.0;
        double r149056 = s;
        double r149057 = -r149056;
        double r149058 = exp(r149057);
        double r149059 = r149055 + r149058;
        double r149060 = r149055 / r149059;
        double r149061 = c_p;
        double r149062 = pow(r149060, r149061);
        double r149063 = r149055 - r149060;
        double r149064 = c_n;
        double r149065 = pow(r149063, r149064);
        double r149066 = r149062 * r149065;
        double r149067 = t;
        double r149068 = -r149067;
        double r149069 = exp(r149068);
        double r149070 = r149055 + r149069;
        double r149071 = r149055 / r149070;
        double r149072 = pow(r149071, r149061);
        double r149073 = r149055 - r149071;
        double r149074 = pow(r149073, r149064);
        double r149075 = r149072 * r149074;
        double r149076 = r149066 / r149075;
        return r149076;
}