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 r13316146 = 1.0;
        double r13316147 = s;
        double r13316148 = -r13316147;
        double r13316149 = exp(r13316148);
        double r13316150 = r13316146 + r13316149;
        double r13316151 = r13316146 / r13316150;
        double r13316152 = c_p;
        double r13316153 = pow(r13316151, r13316152);
        double r13316154 = r13316146 - r13316151;
        double r13316155 = c_n;
        double r13316156 = pow(r13316154, r13316155);
        double r13316157 = r13316153 * r13316156;
        double r13316158 = t;
        double r13316159 = -r13316158;
        double r13316160 = exp(r13316159);
        double r13316161 = r13316146 + r13316160;
        double r13316162 = r13316146 / r13316161;
        double r13316163 = pow(r13316162, r13316152);
        double r13316164 = r13316146 - r13316162;
        double r13316165 = pow(r13316164, r13316155);
        double r13316166 = r13316163 * r13316165;
        double r13316167 = r13316157 / r13316166;
        return r13316167;
}