blob: 2cf7a7b2d5fe92bef89320a8de33f93f8dc22054 (
plain) (
blame)
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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
# URAND -- Uniform random number generator. Write a list of random numbers
# to the standard output. The number of columns per line of output and the
# number of significant digits in each number are parameterized.
procedure t_urand()
int nlines # number of lines of output
int ncols # number of random numbers per line
int ndigits # number of digits of precision
long seed # seed for the random number generator
real scale_factor # scale output numbers by this factor
int n, i
long seed1
int clgeti()
long clgetl(), clktime()
real clgetr(), urand()
begin
# Get parameters from the CL.
nlines = clgeti ("nlines")
ncols = clgeti ("ncols")
ndigits = clgeti ("ndigits")
seed = clgetl ("seed")
scale_factor = clgetr ("scale_factor")
# Compute the random numbers and print on the standard output as
# a list, "ncols" numbers per output line. The output format
# is dependent on the ndigits of precision, set by the user.
if (IS_INDEFL(seed))
seed1 = seed1 + clktime (long(0))
else
seed1 = seed
for (n=1; n <= nlines; n=n+1) {
do i = 1, ncols {
call printf ("%*.*g ")
call pargi (ndigits + 2) # field width
call pargi (ndigits) # precision
call pargr (urand (seed1) * scale_factor)
}
call printf ("\n")
}
end
|
