test.simba

program new;
{$i srl.simba}

function strofchar(s: char; num: Integer): String;
var
  i: Integer;
begin
  Result := '';
  for i := 1 to num do
    Result := Result + s;
end;

function MaxA(a: TIntegerArray): Integer;
var
  L, i: Integer;
begin
  L := High(a);
  if (L < 0) then
  begin
    Result := -1;
    exit;
  end;

  Result := a[0];

  for i := 1 to L do
    if (a[i] > Result) then
      Result := a[i];
end;

procedure plotValues();
var
  x, err, i, ma: Integer;
  vals: TIntegerArray;
  time: Integer;

begin
  setLength(vals, 101);

  time := getSystemTime;
  for i := 1 to 100000 do
  begin

    x := iGaussRange(-50, 50);
    if (inRange(x, -50, 50)) then
      inc(vals[x+50]) else
      inc(err);
  end;
  time := getSystemTime-time;

  ma := MaxA(vals);
  if (ma = 0) then
  begin
    WriteLn('max is 0..');
    Exit;
  end;

  for i := 0 to 100 do
    WriteLn(Format('%4.0d:%s [%d]', [i-50, strofchar('=', round(vals[i]*60.0/ma)), vals[i]]));
  WriteLn(Format('Time taken: %d ms, %.4f/cycle. Errors: %d', 
    [time, time/100000.0, err]));
end;

begin
  setupGauss();

  WriteLn('Gauss Box-Muller Polar:');
  SRL_GaussSettings.genStd := @GaussBMPolar;
  plotValues();

  WriteLn('Gauss n2_gauss');
  SRL_GaussSettings.genStd := @n2_gauss;
  plotValues();


end.