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 r146174 = 1.0;
        double r146175 = s;
        double r146176 = -r146175;
        double r146177 = exp(r146176);
        double r146178 = r146174 + r146177;
        double r146179 = r146174 / r146178;
        double r146180 = c_p;
        double r146181 = pow(r146179, r146180);
        double r146182 = r146174 - r146179;
        double r146183 = c_n;
        double r146184 = pow(r146182, r146183);
        double r146185 = r146181 * r146184;
        double r146186 = t;
        double r146187 = -r146186;
        double r146188 = exp(r146187);
        double r146189 = r146174 + r146188;
        double r146190 = r146174 / r146189;
        double r146191 = pow(r146190, r146180);
        double r146192 = r146174 - r146190;
        double r146193 = pow(r146192, r146183);
        double r146194 = r146191 * r146193;
        double r146195 = r146185 / r146194;
        return r146195;
}