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 r191106 = 1.0;
        double r191107 = s;
        double r191108 = -r191107;
        double r191109 = exp(r191108);
        double r191110 = r191106 + r191109;
        double r191111 = r191106 / r191110;
        double r191112 = c_p;
        double r191113 = pow(r191111, r191112);
        double r191114 = r191106 - r191111;
        double r191115 = c_n;
        double r191116 = pow(r191114, r191115);
        double r191117 = r191113 * r191116;
        double r191118 = t;
        double r191119 = -r191118;
        double r191120 = exp(r191119);
        double r191121 = r191106 + r191120;
        double r191122 = r191106 / r191121;
        double r191123 = pow(r191122, r191112);
        double r191124 = r191106 - r191122;
        double r191125 = pow(r191124, r191115);
        double r191126 = r191123 * r191125;
        double r191127 = r191117 / r191126;
        return r191127;
}