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 r144346 = 1.0;
        double r144347 = s;
        double r144348 = -r144347;
        double r144349 = exp(r144348);
        double r144350 = r144346 + r144349;
        double r144351 = r144346 / r144350;
        double r144352 = c_p;
        double r144353 = pow(r144351, r144352);
        double r144354 = r144346 - r144351;
        double r144355 = c_n;
        double r144356 = pow(r144354, r144355);
        double r144357 = r144353 * r144356;
        double r144358 = t;
        double r144359 = -r144358;
        double r144360 = exp(r144359);
        double r144361 = r144346 + r144360;
        double r144362 = r144346 / r144361;
        double r144363 = pow(r144362, r144352);
        double r144364 = r144346 - r144362;
        double r144365 = pow(r144364, r144355);
        double r144366 = r144363 * r144365;
        double r144367 = r144357 / r144366;
        return r144367;
}