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 r139841 = 1.0;
        double r139842 = s;
        double r139843 = -r139842;
        double r139844 = exp(r139843);
        double r139845 = r139841 + r139844;
        double r139846 = r139841 / r139845;
        double r139847 = c_p;
        double r139848 = pow(r139846, r139847);
        double r139849 = r139841 - r139846;
        double r139850 = c_n;
        double r139851 = pow(r139849, r139850);
        double r139852 = r139848 * r139851;
        double r139853 = t;
        double r139854 = -r139853;
        double r139855 = exp(r139854);
        double r139856 = r139841 + r139855;
        double r139857 = r139841 / r139856;
        double r139858 = pow(r139857, r139847);
        double r139859 = r139841 - r139857;
        double r139860 = pow(r139859, r139850);
        double r139861 = r139858 * r139860;
        double r139862 = r139852 / r139861;
        return r139862;
}