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 r99108 = 1.0;
        double r99109 = s;
        double r99110 = -r99109;
        double r99111 = exp(r99110);
        double r99112 = r99108 + r99111;
        double r99113 = r99108 / r99112;
        double r99114 = c_p;
        double r99115 = pow(r99113, r99114);
        double r99116 = r99108 - r99113;
        double r99117 = c_n;
        double r99118 = pow(r99116, r99117);
        double r99119 = r99115 * r99118;
        double r99120 = t;
        double r99121 = -r99120;
        double r99122 = exp(r99121);
        double r99123 = r99108 + r99122;
        double r99124 = r99108 / r99123;
        double r99125 = pow(r99124, r99114);
        double r99126 = r99108 - r99124;
        double r99127 = pow(r99126, r99117);
        double r99128 = r99125 * r99127;
        double r99129 = r99119 / r99128;
        return r99129;
}