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 r155021 = 1.0;
        double r155022 = s;
        double r155023 = -r155022;
        double r155024 = exp(r155023);
        double r155025 = r155021 + r155024;
        double r155026 = r155021 / r155025;
        double r155027 = c_p;
        double r155028 = pow(r155026, r155027);
        double r155029 = r155021 - r155026;
        double r155030 = c_n;
        double r155031 = pow(r155029, r155030);
        double r155032 = r155028 * r155031;
        double r155033 = t;
        double r155034 = -r155033;
        double r155035 = exp(r155034);
        double r155036 = r155021 + r155035;
        double r155037 = r155021 / r155036;
        double r155038 = pow(r155037, r155027);
        double r155039 = r155021 - r155037;
        double r155040 = pow(r155039, r155030);
        double r155041 = r155038 * r155040;
        double r155042 = r155032 / r155041;
        return r155042;
}