poisson_rnd() excruciatingly slow

[mddejong]

Is there a reason why poisson_rnd() is so *glacially slow* for high values
of lambda? On a 2.0 GHz Athlon system, I measured poisson_rnd(200000) to
take approximately 18 seconds, making this function next to useless for
practical use. Compare that to the 0.16 seconds for binomial_rnd(400000,
0.5) which effectively computes the same thing, which I still consider to
be a sizable amount, by the way.

If my questions have already been answered on various lists, then I
apologize for asking them, but as the search-function of said lists seems
to be down, I cannot check myself.


Kind regards,
Maarten



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------