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 r147335 = 1.0;
        double r147336 = s;
        double r147337 = -r147336;
        double r147338 = exp(r147337);
        double r147339 = r147335 + r147338;
        double r147340 = r147335 / r147339;
        double r147341 = c_p;
        double r147342 = pow(r147340, r147341);
        double r147343 = r147335 - r147340;
        double r147344 = c_n;
        double r147345 = pow(r147343, r147344);
        double r147346 = r147342 * r147345;
        double r147347 = t;
        double r147348 = -r147347;
        double r147349 = exp(r147348);
        double r147350 = r147335 + r147349;
        double r147351 = r147335 / r147350;
        double r147352 = pow(r147351, r147341);
        double r147353 = r147335 - r147351;
        double r147354 = pow(r147353, r147344);
        double r147355 = r147352 * r147354;
        double r147356 = r147346 / r147355;
        return r147356;
}