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 r155880 = 1.0;
        double r155881 = s;
        double r155882 = -r155881;
        double r155883 = exp(r155882);
        double r155884 = r155880 + r155883;
        double r155885 = r155880 / r155884;
        double r155886 = c_p;
        double r155887 = pow(r155885, r155886);
        double r155888 = r155880 - r155885;
        double r155889 = c_n;
        double r155890 = pow(r155888, r155889);
        double r155891 = r155887 * r155890;
        double r155892 = t;
        double r155893 = -r155892;
        double r155894 = exp(r155893);
        double r155895 = r155880 + r155894;
        double r155896 = r155880 / r155895;
        double r155897 = pow(r155896, r155886);
        double r155898 = r155880 - r155896;
        double r155899 = pow(r155898, r155889);
        double r155900 = r155897 * r155899;
        double r155901 = r155891 / r155900;
        return r155901;
}