| 63 | } |
| 64 | |
| 65 | inline double plnorm_0(double x, int lower_tail, int log_p){ |
| 66 | #ifdef IEEE_754 |
| 67 | if (ISNAN(x)) |
| 68 | return x + 1.0; |
| 69 | #endif |
| 70 | |
| 71 | if (x > 0) |
| 72 | return Rcpp::stats::pnorm_0(::log(x), lower_tail, log_p); |
| 73 | return R_DT_0; |
| 74 | } |
| 75 | |
| 76 | inline double plnorm_1(double x, double meanlog, int lower_tail, int log_p) { |
| 77 | #ifdef IEEE_754 |