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 r153627 = 1.0;
        double r153628 = s;
        double r153629 = -r153628;
        double r153630 = exp(r153629);
        double r153631 = r153627 + r153630;
        double r153632 = r153627 / r153631;
        double r153633 = c_p;
        double r153634 = pow(r153632, r153633);
        double r153635 = r153627 - r153632;
        double r153636 = c_n;
        double r153637 = pow(r153635, r153636);
        double r153638 = r153634 * r153637;
        double r153639 = t;
        double r153640 = -r153639;
        double r153641 = exp(r153640);
        double r153642 = r153627 + r153641;
        double r153643 = r153627 / r153642;
        double r153644 = pow(r153643, r153633);
        double r153645 = r153627 - r153643;
        double r153646 = pow(r153645, r153636);
        double r153647 = r153644 * r153646;
        double r153648 = r153638 / r153647;
        return r153648;
}