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 r156459 = 1.0;
        double r156460 = s;
        double r156461 = -r156460;
        double r156462 = exp(r156461);
        double r156463 = r156459 + r156462;
        double r156464 = r156459 / r156463;
        double r156465 = c_p;
        double r156466 = pow(r156464, r156465);
        double r156467 = r156459 - r156464;
        double r156468 = c_n;
        double r156469 = pow(r156467, r156468);
        double r156470 = r156466 * r156469;
        double r156471 = t;
        double r156472 = -r156471;
        double r156473 = exp(r156472);
        double r156474 = r156459 + r156473;
        double r156475 = r156459 / r156474;
        double r156476 = pow(r156475, r156465);
        double r156477 = r156459 - r156475;
        double r156478 = pow(r156477, r156468);
        double r156479 = r156476 * r156478;
        double r156480 = r156470 / r156479;
        return r156480;
}