Return to multprec_random_numbers.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Numbers

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Numbers/multprec_random_numbers.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
        !             2: with Standard_Random_Numbers;
        !             3:
        !             4: package body Multprec_Random_Numbers is
        !             5:
        !             6:   function Random ( size : natural ) return Natural_Number is
        !             7:
        !             8:     res : Natural_Number;
        !             9:     arn : Array_of_Naturals(0..size);
        !            10:
        !            11:   begin
        !            12:     for i in arn'range loop
        !            13:       arn(i) := Standard_Random_Numbers.Random(0,Base-1);
        !            14:     end loop;
        !            15:     res := Create(arn);
        !            16:     return res;
        !            17:   end Random;
        !            18:
        !            19:   function Random ( size : natural ) return Integer_Number is
        !            20:
        !            21:     res : Integer_Number;
        !            22:     sig : double_float := Standard_Random_Numbers.Random;
        !            23:     rnd : Natural_Number := Random(size);
        !            24:
        !            25:   begin
        !            26:     res := Create(rnd);
        !            27:     if sig < 0.0
        !            28:      then Min(res);
        !            29:     end if;
        !            30:     Clear(rnd);
        !            31:     return res;
        !            32:   end Random;
        !            33:
        !            34:   function Random ( size : natural ) return Floating_Number is
        !            35:
        !            36:     res : Floating_Number;
        !            37:     naturfrac : Natural_Number := Random(size);
        !            38:     frac,expo : Integer_Number;
        !            39:
        !            40:   begin
        !            41:     frac := Create(naturfrac);
        !            42:     Clear(naturfrac);
        !            43:     expo := Create(-Decimal_Places(frac));
        !            44:     res := Create(frac,expo);
        !            45:     Clear(frac); Clear(expo);
        !            46:     return res;
        !            47:   end Random;
        !            48:
        !            49:   function Random ( size : natural ) return Complex_Number is
        !            50:
        !            51:     res : Complex_Number;
        !            52:     re : Floating_Number := Random(size);
        !            53:     im : Floating_Number := Random(size);
        !            54:
        !            55:   begin
        !            56:     res := Create(re,im);
        !            57:     Clear(re); Clear(im);
        !            58:     return res;
        !            59:   end Random;
        !            60:
        !            61: end Multprec_Random_Numbers;