| 85 | } |
| 86 | |
| 87 | inline double qlnorm_0(double p, int lower_tail, int log_p){ |
| 88 | #ifdef IEEE_754 |
| 89 | if (ISNAN(p)) |
| 90 | return p + 1.0; |
| 91 | #endif |
| 92 | R_Q_P01_boundaries(p, 0, ML_POSINF); |
| 93 | |
| 94 | return ::exp(::Rf_qnorm5(p, 0.0, 1.0, lower_tail, log_p)); |
| 95 | } |
| 96 | |
| 97 | inline double qlnorm_1(double p, double meanlog, int lower_tail, int log_p){ |
| 98 | #ifdef IEEE_754 |