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 r173168 = 1.0;
        double r173169 = s;
        double r173170 = -r173169;
        double r173171 = exp(r173170);
        double r173172 = r173168 + r173171;
        double r173173 = r173168 / r173172;
        double r173174 = c_p;
        double r173175 = pow(r173173, r173174);
        double r173176 = r173168 - r173173;
        double r173177 = c_n;
        double r173178 = pow(r173176, r173177);
        double r173179 = r173175 * r173178;
        double r173180 = t;
        double r173181 = -r173180;
        double r173182 = exp(r173181);
        double r173183 = r173168 + r173182;
        double r173184 = r173168 / r173183;
        double r173185 = pow(r173184, r173174);
        double r173186 = r173168 - r173184;
        double r173187 = pow(r173186, r173177);
        double r173188 = r173185 * r173187;
        double r173189 = r173179 / r173188;
        return r173189;
}