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 r164130 = 1.0;
        double r164131 = s;
        double r164132 = -r164131;
        double r164133 = exp(r164132);
        double r164134 = r164130 + r164133;
        double r164135 = r164130 / r164134;
        double r164136 = c_p;
        double r164137 = pow(r164135, r164136);
        double r164138 = r164130 - r164135;
        double r164139 = c_n;
        double r164140 = pow(r164138, r164139);
        double r164141 = r164137 * r164140;
        double r164142 = t;
        double r164143 = -r164142;
        double r164144 = exp(r164143);
        double r164145 = r164130 + r164144;
        double r164146 = r164130 / r164145;
        double r164147 = pow(r164146, r164136);
        double r164148 = r164130 - r164146;
        double r164149 = pow(r164148, r164139);
        double r164150 = r164147 * r164149;
        double r164151 = r164141 / r164150;
        return r164151;
}