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 r246852 = 1.0;
        double r246853 = s;
        double r246854 = -r246853;
        double r246855 = exp(r246854);
        double r246856 = r246852 + r246855;
        double r246857 = r246852 / r246856;
        double r246858 = c_p;
        double r246859 = pow(r246857, r246858);
        double r246860 = r246852 - r246857;
        double r246861 = c_n;
        double r246862 = pow(r246860, r246861);
        double r246863 = r246859 * r246862;
        double r246864 = t;
        double r246865 = -r246864;
        double r246866 = exp(r246865);
        double r246867 = r246852 + r246866;
        double r246868 = r246852 / r246867;
        double r246869 = pow(r246868, r246858);
        double r246870 = r246852 - r246868;
        double r246871 = pow(r246870, r246861);
        double r246872 = r246869 * r246871;
        double r246873 = r246863 / r246872;
        return r246873;
}