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 r197140 = 1.0;
        double r197141 = s;
        double r197142 = -r197141;
        double r197143 = exp(r197142);
        double r197144 = r197140 + r197143;
        double r197145 = r197140 / r197144;
        double r197146 = c_p;
        double r197147 = pow(r197145, r197146);
        double r197148 = r197140 - r197145;
        double r197149 = c_n;
        double r197150 = pow(r197148, r197149);
        double r197151 = r197147 * r197150;
        double r197152 = t;
        double r197153 = -r197152;
        double r197154 = exp(r197153);
        double r197155 = r197140 + r197154;
        double r197156 = r197140 / r197155;
        double r197157 = pow(r197156, r197146);
        double r197158 = r197140 - r197156;
        double r197159 = pow(r197158, r197149);
        double r197160 = r197157 * r197159;
        double r197161 = r197151 / r197160;
        return r197161;
}