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 r7612632 = 1.0;
        double r7612633 = s;
        double r7612634 = -r7612633;
        double r7612635 = exp(r7612634);
        double r7612636 = r7612632 + r7612635;
        double r7612637 = r7612632 / r7612636;
        double r7612638 = c_p;
        double r7612639 = pow(r7612637, r7612638);
        double r7612640 = r7612632 - r7612637;
        double r7612641 = c_n;
        double r7612642 = pow(r7612640, r7612641);
        double r7612643 = r7612639 * r7612642;
        double r7612644 = t;
        double r7612645 = -r7612644;
        double r7612646 = exp(r7612645);
        double r7612647 = r7612632 + r7612646;
        double r7612648 = r7612632 / r7612647;
        double r7612649 = pow(r7612648, r7612638);
        double r7612650 = r7612632 - r7612648;
        double r7612651 = pow(r7612650, r7612641);
        double r7612652 = r7612649 * r7612651;
        double r7612653 = r7612643 / r7612652;
        return r7612653;
}