| 283 | } |
| 284 | |
| 285 | double randp(double scale, double shape) { |
| 286 | std::exponential_distribution<double> d(shape); |
| 287 | return scale * exp(d(detail::default_generator())); |
| 288 | } |
| 289 | |
| 290 | std::vector<double> randp(size_t n, double scale, double shape) { |
| 291 | std::vector<double> x(n); |