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 r131849 = 1.0;
        double r131850 = s;
        double r131851 = -r131850;
        double r131852 = exp(r131851);
        double r131853 = r131849 + r131852;
        double r131854 = r131849 / r131853;
        double r131855 = c_p;
        double r131856 = pow(r131854, r131855);
        double r131857 = r131849 - r131854;
        double r131858 = c_n;
        double r131859 = pow(r131857, r131858);
        double r131860 = r131856 * r131859;
        double r131861 = t;
        double r131862 = -r131861;
        double r131863 = exp(r131862);
        double r131864 = r131849 + r131863;
        double r131865 = r131849 / r131864;
        double r131866 = pow(r131865, r131855);
        double r131867 = r131849 - r131865;
        double r131868 = pow(r131867, r131858);
        double r131869 = r131866 * r131868;
        double r131870 = r131860 / r131869;
        return r131870;
}