Re: Discrete random with given distribution ?

[Jeffrey Carter <jrcarter@acm.org>]

Wes Groleau wrote:
> 
>    Threshold : array (Out_Type) of Float :=
>      (A => 0.50,
>       B => 0.75,   -- previous plus Probability (B)
>       C => 1.00);  -- previous plus Probability (C)
> 
>    F := FR.Random;
>    Loop_With_Jump_Out:
>      for I in Out_Type loop
>         if F > Threshold (I) then
>            return I;
>         end if;
>      end loop Loop_With_Jump_Out;

This can never return C. F <= 1.0, and for many generators, F < 1.0.
Also, if F = 0.9, this returns A.

I suspect you need

F < Threshold (I)

-- 
Jeff Carter
"Oh Lord, bless this thy hand grenade, that with it thou
mayst blow thine enemies to tiny bits, in thy mercy."
Monty Python and the Holy Grail