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 r178960 = 1.0;
        double r178961 = s;
        double r178962 = -r178961;
        double r178963 = exp(r178962);
        double r178964 = r178960 + r178963;
        double r178965 = r178960 / r178964;
        double r178966 = c_p;
        double r178967 = pow(r178965, r178966);
        double r178968 = r178960 - r178965;
        double r178969 = c_n;
        double r178970 = pow(r178968, r178969);
        double r178971 = r178967 * r178970;
        double r178972 = t;
        double r178973 = -r178972;
        double r178974 = exp(r178973);
        double r178975 = r178960 + r178974;
        double r178976 = r178960 / r178975;
        double r178977 = pow(r178976, r178966);
        double r178978 = r178960 - r178976;
        double r178979 = pow(r178978, r178969);
        double r178980 = r178977 * r178979;
        double r178981 = r178971 / r178980;
        return r178981;
}