← Back to branch summary

~ubuntu-branches/ubuntu/gutsy/avr-libc/gutsy

« back to all changes in this revision

Viewing changes to libc/stdlib/rand.c

  • Committer: Bazaar Package Importer
  • Author(s): Hakan Ardo
  • Date: 2005-03-19 11:16:14 UTC
  • mfrom: (1.1.1 upstream) (2.1.2 hoary)
  • Revision ID: james.westby@ubuntu.com-20050319111614-4g01s2ftv5x5nxf3
Tags: 1:1.2.3-3
* Added build depends on netpbm
* Added build depends on tetex-extra

Show diffs side-by-side

added added

removed removed

Lines of Context:
libc/stdlib/rand.c
10
10
 * 2. Redistributions in binary form must reproduce the above copyright
11
11
 *    notice, this list of conditions and the following disclaimer in the
12
12
 *    documentation and/or other materials provided with the distribution.
13
 
 * 3. All advertising materials mentioning features or use of this software
14
 
 *    must display the following acknowledgement:
15
 
 *      This product includes software developed by the University of
16
 
 *      California, Berkeley and its contributors.
17
 
 * 4. Neither the name of the University nor the names of its contributors
 
13
 * 3. Neither the name of the University nor the names of its contributors
18
14
 *    may be used to endorse or promote products derived from this software
19
15
 *    without specific prior written permission.
20
16
 *
31
27
 * SUCH DAMAGE.
32
28
 *
33
29
 * Posix rand_r function added May 1999 by Wes Peters <wes@softweyr.com>.
 
30
 *
 
31
 * $Id: rand.c,v 1.2 2002/10/16 15:52:25 joerg_wunsch Exp $
34
32
 */
35
33
 
36
 
#if defined(LIBC_SCCS) && !defined(lint)
 
34
/*
 
35
 * From:
37
36
static char sccsid[] = "@(#)rand.c      8.1 (Berkeley) 6/14/93";
38
 
#endif /* LIBC_SCCS and not lint */
 
37
*/
39
38
 
40
 
#if 0
41
 
#include <sys/types.h>
42
 
#endif
43
39
#include <stdlib.h>
44
40
 
45
 
#ifdef TEST
46
 
#include <stdio.h>
47
 
#endif /* TEST */
48
 
 
49
41
static int
50
42
do_rand(unsigned long *ctx)
51
43
{
52
 
        return ((*ctx = *ctx * 1103515245 + 12345) % ((unsigned long int)RAND_MAX + 1));
 
44
#ifdef  USE_WEAK_SEEDING
 
45
        /*
 
46
         * Historic implementation compatibility.
 
47
         * The random sequences do not vary much with the seed,
 
48
         * even with overflowing.
 
49
         */
 
50
        return ((*ctx = *ctx * 1103515245L + 12345L) %
 
51
                ((unsigned long)RAND_MAX + 1));
 
52
#else   /* !USE_WEAK_SEEDING */
 
53
        /*
 
54
         * Compute x = (7^5 * x) mod (2^31 - 1)
 
55
         * wihout overflowing 31 bits:
 
56
         *      (2^31 - 1) = 127773 * (7^5) + 2836
 
57
         * From "Random number generators: good ones are hard to find",
 
58
         * Park and Miller, Communications of the ACM, vol. 31, no. 10,
 
59
         * October 1988, p. 1195.
 
60
         */
 
61
        long hi, lo, x;
 
62
 
 
63
        hi = *ctx / 127773L;
 
64
        lo = *ctx % 127773L;
 
65
        x = 16807L * lo - 2836L * hi;
 
66
        if (x <= 0)
 
67
                x += 0x7fffffffL;
 
68
        return ((*ctx = x) % ((unsigned long)RAND_MAX + 1));
 
69
#endif  /* !USE_WEAK_SEEDING */
53
70
}
54
71
 
55
72
 
56
73
int
57
 
rand_r(unsigned int *ctx)
 
74
rand_r(unsigned long *ctx)
58
75
{
59
 
        unsigned long int val = (unsigned long int) *ctx;
60
 
        *ctx = do_rand(&val);
61
 
        return (int) *ctx;
 
76
        return do_rand(ctx);
62
77
}
63
78
 
64
79
 
65
 
static unsigned long int next = 1;
 
80
static unsigned long next = 1;
66
81
 
67
82
int
68
83
rand(void)
76
91
        next = seed;
77
92
}
78
93
 
79
 
#ifdef TEST
80
 
 
81
 
main()
82
 
{
83
 
    int i;
84
 
    unsigned myseed;
85
 
 
86
 
    printf("seeding rand with 0x19610910: \n");
87
 
    srand(0x19610910);
88
 
 
89
 
    printf("generating three pseudo-random numbers:\n");
90
 
    for (i = 0; i < 3; i++)
91
 
    {
92
 
        printf("next random number = %d\n", rand());
93
 
    }
94
 
 
95
 
    printf("generating the same sequence with rand_r:\n");
96
 
    myseed = 0x19610910;
97
 
    for (i = 0; i < 3; i++)
98
 
    {
99
 
        printf("next random number = %d\n", rand_r(&myseed));
100
 
    }
101
 
 
102
 
    return 0;
103
 
}
104
 
 
105
 
#endif /* TEST */
106