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 r160859 = 1.0;
        double r160860 = s;
        double r160861 = -r160860;
        double r160862 = exp(r160861);
        double r160863 = r160859 + r160862;
        double r160864 = r160859 / r160863;
        double r160865 = c_p;
        double r160866 = pow(r160864, r160865);
        double r160867 = r160859 - r160864;
        double r160868 = c_n;
        double r160869 = pow(r160867, r160868);
        double r160870 = r160866 * r160869;
        double r160871 = t;
        double r160872 = -r160871;
        double r160873 = exp(r160872);
        double r160874 = r160859 + r160873;
        double r160875 = r160859 / r160874;
        double r160876 = pow(r160875, r160865);
        double r160877 = r160859 - r160875;
        double r160878 = pow(r160877, r160868);
        double r160879 = r160876 * r160878;
        double r160880 = r160870 / r160879;
        return r160880;
}