| 74 | } |
| 75 | |
| 76 | inline double plnorm_1(double x, double meanlog, int lower_tail, int log_p) { |
| 77 | #ifdef IEEE_754 |
| 78 | if (ISNAN(x) || ISNAN(meanlog)) |
| 79 | return x + meanlog + 1.0; |
| 80 | #endif |
| 81 | |
| 82 | if (x > 0) |
| 83 | return Rcpp::stats::pnorm_1(::log(x), meanlog, lower_tail, log_p); |
| 84 | return R_DT_0; |
| 85 | } |
| 86 | |
| 87 | inline double qlnorm_0(double p, int lower_tail, int log_p){ |
| 88 | #ifdef IEEE_754 |