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 r168096 = 1.0;
        double r168097 = s;
        double r168098 = -r168097;
        double r168099 = exp(r168098);
        double r168100 = r168096 + r168099;
        double r168101 = r168096 / r168100;
        double r168102 = c_p;
        double r168103 = pow(r168101, r168102);
        double r168104 = r168096 - r168101;
        double r168105 = c_n;
        double r168106 = pow(r168104, r168105);
        double r168107 = r168103 * r168106;
        double r168108 = t;
        double r168109 = -r168108;
        double r168110 = exp(r168109);
        double r168111 = r168096 + r168110;
        double r168112 = r168096 / r168111;
        double r168113 = pow(r168112, r168102);
        double r168114 = r168096 - r168112;
        double r168115 = pow(r168114, r168105);
        double r168116 = r168113 * r168115;
        double r168117 = r168107 / r168116;
        return r168117;
}