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 r169847 = 1.0;
        double r169848 = s;
        double r169849 = -r169848;
        double r169850 = exp(r169849);
        double r169851 = r169847 + r169850;
        double r169852 = r169847 / r169851;
        double r169853 = c_p;
        double r169854 = pow(r169852, r169853);
        double r169855 = r169847 - r169852;
        double r169856 = c_n;
        double r169857 = pow(r169855, r169856);
        double r169858 = r169854 * r169857;
        double r169859 = t;
        double r169860 = -r169859;
        double r169861 = exp(r169860);
        double r169862 = r169847 + r169861;
        double r169863 = r169847 / r169862;
        double r169864 = pow(r169863, r169853);
        double r169865 = r169847 - r169863;
        double r169866 = pow(r169865, r169856);
        double r169867 = r169864 * r169866;
        double r169868 = r169858 / r169867;
        return r169868;
}