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 r21030478 = 1.0;
        double r21030479 = s;
        double r21030480 = -r21030479;
        double r21030481 = exp(r21030480);
        double r21030482 = r21030478 + r21030481;
        double r21030483 = r21030478 / r21030482;
        double r21030484 = c_p;
        double r21030485 = pow(r21030483, r21030484);
        double r21030486 = r21030478 - r21030483;
        double r21030487 = c_n;
        double r21030488 = pow(r21030486, r21030487);
        double r21030489 = r21030485 * r21030488;
        double r21030490 = t;
        double r21030491 = -r21030490;
        double r21030492 = exp(r21030491);
        double r21030493 = r21030478 + r21030492;
        double r21030494 = r21030478 / r21030493;
        double r21030495 = pow(r21030494, r21030484);
        double r21030496 = r21030478 - r21030494;
        double r21030497 = pow(r21030496, r21030487);
        double r21030498 = r21030495 * r21030497;
        double r21030499 = r21030489 / r21030498;
        return r21030499;
}