adagsl_335d13f0/examples/randist/src/randist.adb

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
with Ada.Text_Io; use Ada.Text_IO ;
with Ada.Long_Float_Text_IO ; use Ada.Long_Float_Text_IO ;
with Ada.Command_Line; use Ada.Command_Line;

with Interfaces.C; use Interfaces.C ;
with Interfaces.C.Strings; use Interfaces.C.Strings ;
with gsl ;
with gsl.rng ;
with gsl.randist ;

procedure randist is

   rng : access gsl.rng.gsl_rng ; 
   sigma : double := 3.0 ;

   rn : double ;
   outfile : File_Type ;

   procedure Create(fn : String ) is
   begin
      Create(outfile, Out_File, fn );
      Set_Output(outfile);
   end Create ;
   procedure Close is
   begin
      Close(outfile);
      Set_Output(Standard_Output);
   end Close;

   procedure TestGaussian is
   begin 
      Put(Interfaces.C.Strings.Value( gsl.rng.name(rng) ));
      New_Line;
      Create (Value(gsl.rng.name(rng)) & ".txt");
      for i in 1..1024
      loop
         rn := gsl.randist.gaussian(rng,sigma);
         Put(Long_Float(rn));
         New_Line;
      end loop ;
      Close ;
   end TestGaussian ;

   procedure TestGaussianPdf is
      x : double := -6.0 ;
      xdelta : double := 2.0 * abs(x) / 1024.0 ;
      pd : double ;
   begin 
      Put(Interfaces.C.Strings.Value( gsl.rng.name(rng) ));
      New_Line;
      Create (Value(gsl.rng.name(rng)) & ".pdf.txt");
      for i in 1..1024
      loop
         pd := gsl.randist.gaussian_pdf(x,sigma);
         Put(Long_Float(x)); Put(";");
         Put(Long_Float(pd));
         New_Line;
         x := x + xdelta ;
      end loop ;
      Close ;
   end TestGaussianPdf ;

begin
   Put_Line("Default Random Number Generator ");
   Put_Line("Setup from environment");
   rng := gsl.rng.alloc(gsl.rng.env_setup);
   if Argument_Count > 0
   then
      sigma := double'Value(Argument(1));
   end if;
   TestGaussian;
   TestGaussianPdf;
   gsl.rng.free(rng);
end randist ;