...I have "borrowed" the code Jochen
Voss wrote for gsl to implement the ziggurat method for generating
random normal values. Here's what I'm seeing after making this
change:
loop vectorized
octave 2.59 0.53
gretl 0.87 0.39
ox 0.85 0.56
Wow!
--
“Remember not only to say the right thing in the right place, but far
more difficult still, to leave unsaid the wrong thing at the tempting
moment.” - Benjamin Franklin (1706-1790)
--