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 r150506 = 1.0;
        double r150507 = s;
        double r150508 = -r150507;
        double r150509 = exp(r150508);
        double r150510 = r150506 + r150509;
        double r150511 = r150506 / r150510;
        double r150512 = c_p;
        double r150513 = pow(r150511, r150512);
        double r150514 = r150506 - r150511;
        double r150515 = c_n;
        double r150516 = pow(r150514, r150515);
        double r150517 = r150513 * r150516;
        double r150518 = t;
        double r150519 = -r150518;
        double r150520 = exp(r150519);
        double r150521 = r150506 + r150520;
        double r150522 = r150506 / r150521;
        double r150523 = pow(r150522, r150512);
        double r150524 = r150506 - r150522;
        double r150525 = pow(r150524, r150515);
        double r150526 = r150523 * r150525;
        double r150527 = r150517 / r150526;
        return r150527;
}