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 r45835162 = 1.0;
        double r45835163 = s;
        double r45835164 = -r45835163;
        double r45835165 = exp(r45835164);
        double r45835166 = r45835162 + r45835165;
        double r45835167 = r45835162 / r45835166;
        double r45835168 = c_p;
        double r45835169 = pow(r45835167, r45835168);
        double r45835170 = r45835162 - r45835167;
        double r45835171 = c_n;
        double r45835172 = pow(r45835170, r45835171);
        double r45835173 = r45835169 * r45835172;
        double r45835174 = t;
        double r45835175 = -r45835174;
        double r45835176 = exp(r45835175);
        double r45835177 = r45835162 + r45835176;
        double r45835178 = r45835162 / r45835177;
        double r45835179 = pow(r45835178, r45835168);
        double r45835180 = r45835162 - r45835178;
        double r45835181 = pow(r45835180, r45835171);
        double r45835182 = r45835179 * r45835181;
        double r45835183 = r45835173 / r45835182;
        return r45835183;
}