arc4random.c revision 1.9 1 1.9 perry /* $NetBSD: arc4random.c,v 1.9 2005/12/24 21:11:16 perry Exp $ */
2 1.1 itojun /* $OpenBSD: arc4random.c,v 1.6 2001/06/05 05:05:38 pvalchev Exp $ */
3 1.1 itojun
4 1.1 itojun /*
5 1.1 itojun * Arc4 random number generator for OpenBSD.
6 1.1 itojun * Copyright 1996 David Mazieres <dm (at) lcs.mit.edu>.
7 1.1 itojun *
8 1.1 itojun * Modification and redistribution in source and binary forms is
9 1.1 itojun * permitted provided that due credit is given to the author and the
10 1.1 itojun * OpenBSD project by leaving this copyright notice intact.
11 1.1 itojun */
12 1.1 itojun
13 1.1 itojun /*
14 1.1 itojun * This code is derived from section 17.1 of Applied Cryptography,
15 1.1 itojun * second edition, which describes a stream cipher allegedly
16 1.1 itojun * compatible with RSA Labs "RC4" cipher (the actual description of
17 1.1 itojun * which is a trade secret). The same algorithm is used as a stream
18 1.1 itojun * cipher called "arcfour" in Tatu Ylonen's ssh package.
19 1.1 itojun *
20 1.1 itojun * Here the stream cipher has been modified always to include the time
21 1.1 itojun * when initializing the state. That makes it impossible to
22 1.1 itojun * regenerate the same random sequence twice, so this can't be used
23 1.1 itojun * for encryption, but will generate good random numbers.
24 1.1 itojun *
25 1.1 itojun * RC4 is a registered trademark of RSA Laboratories.
26 1.1 itojun */
27 1.1 itojun
28 1.8 lukem #include <sys/cdefs.h>
29 1.8 lukem #if defined(LIBC_SCCS) && !defined(lint)
30 1.9 perry __RCSID("$NetBSD: arc4random.c,v 1.9 2005/12/24 21:11:16 perry Exp $");
31 1.8 lukem #endif /* LIBC_SCCS and not lint */
32 1.8 lukem
33 1.7 kleink #include "namespace.h"
34 1.1 itojun #include <fcntl.h>
35 1.1 itojun #include <stdlib.h>
36 1.1 itojun #include <unistd.h>
37 1.1 itojun #include <sys/types.h>
38 1.1 itojun #include <sys/param.h>
39 1.1 itojun #include <sys/time.h>
40 1.1 itojun #include <sys/sysctl.h>
41 1.1 itojun
42 1.7 kleink #ifdef __weak_alias
43 1.7 kleink __weak_alias(arc4random,_arc4random)
44 1.7 kleink #endif
45 1.7 kleink
46 1.1 itojun struct arc4_stream {
47 1.1 itojun u_int8_t i;
48 1.1 itojun u_int8_t j;
49 1.1 itojun u_int8_t s[256];
50 1.1 itojun };
51 1.1 itojun
52 1.6 kleink static int rs_initialized;
53 1.1 itojun static struct arc4_stream rs;
54 1.1 itojun
55 1.1 itojun static inline void arc4_init(struct arc4_stream *);
56 1.1 itojun static inline void arc4_addrandom(struct arc4_stream *, u_char *, int);
57 1.1 itojun static void arc4_stir(struct arc4_stream *);
58 1.1 itojun static inline u_int8_t arc4_getbyte(struct arc4_stream *);
59 1.1 itojun static inline u_int32_t arc4_getword(struct arc4_stream *);
60 1.1 itojun
61 1.1 itojun static inline void
62 1.1 itojun arc4_init(as)
63 1.1 itojun struct arc4_stream *as;
64 1.1 itojun {
65 1.1 itojun int n;
66 1.1 itojun
67 1.1 itojun for (n = 0; n < 256; n++)
68 1.1 itojun as->s[n] = n;
69 1.1 itojun as->i = 0;
70 1.1 itojun as->j = 0;
71 1.1 itojun }
72 1.1 itojun
73 1.1 itojun static inline void
74 1.1 itojun arc4_addrandom(as, dat, datlen)
75 1.1 itojun struct arc4_stream *as;
76 1.1 itojun u_char *dat;
77 1.1 itojun int datlen;
78 1.1 itojun {
79 1.1 itojun int n;
80 1.1 itojun u_int8_t si;
81 1.1 itojun
82 1.1 itojun as->i--;
83 1.1 itojun for (n = 0; n < 256; n++) {
84 1.1 itojun as->i = (as->i + 1);
85 1.1 itojun si = as->s[as->i];
86 1.1 itojun as->j = (as->j + si + dat[n % datlen]);
87 1.1 itojun as->s[as->i] = as->s[as->j];
88 1.1 itojun as->s[as->j] = si;
89 1.1 itojun }
90 1.1 itojun as->j = as->i;
91 1.1 itojun }
92 1.1 itojun
93 1.1 itojun static void
94 1.1 itojun arc4_stir(as)
95 1.1 itojun struct arc4_stream *as;
96 1.1 itojun {
97 1.1 itojun int fd;
98 1.1 itojun struct {
99 1.1 itojun struct timeval tv;
100 1.1 itojun u_int rnd[(128 - sizeof(struct timeval)) / sizeof(u_int)];
101 1.1 itojun } rdat;
102 1.3 itojun int n;
103 1.1 itojun
104 1.1 itojun gettimeofday(&rdat.tv, NULL);
105 1.2 itojun fd = open("/dev/urandom", O_RDONLY);
106 1.1 itojun if (fd != -1) {
107 1.1 itojun read(fd, rdat.rnd, sizeof(rdat.rnd));
108 1.1 itojun close(fd);
109 1.1 itojun }
110 1.4 itojun #ifdef KERN_URND
111 1.1 itojun else {
112 1.5 thorpej int mib[2];
113 1.5 thorpej u_int i;
114 1.1 itojun size_t len;
115 1.1 itojun
116 1.1 itojun /* Device could not be opened, we might be chrooted, take
117 1.1 itojun * randomness from sysctl. */
118 1.1 itojun
119 1.1 itojun mib[0] = CTL_KERN;
120 1.4 itojun mib[1] = KERN_URND;
121 1.1 itojun
122 1.5 thorpej for (i = 0; i < sizeof(rdat.rnd) / sizeof(u_int); i++) {
123 1.1 itojun len = sizeof(u_int);
124 1.1 itojun if (sysctl(mib, 2, &rdat.rnd[i], &len, NULL, 0) == -1)
125 1.1 itojun break;
126 1.1 itojun }
127 1.1 itojun }
128 1.1 itojun #endif
129 1.1 itojun /* fd < 0 or failed sysctl ? Ah, what the heck. We'll just take
130 1.1 itojun * whatever was on the stack... */
131 1.1 itojun
132 1.1 itojun arc4_addrandom(as, (void *) &rdat, sizeof(rdat));
133 1.3 itojun
134 1.3 itojun /*
135 1.3 itojun * Throw away the first N words of output, as suggested in the
136 1.3 itojun * paper "Weaknesses in the Key Scheduling Algorithm of RC4"
137 1.3 itojun * by Fluher, Mantin, and Shamir. (N = 256 in our case.)
138 1.3 itojun */
139 1.3 itojun for (n = 0; n < 256 * 4; n++)
140 1.3 itojun arc4_getbyte(as);
141 1.1 itojun }
142 1.1 itojun
143 1.1 itojun static inline u_int8_t
144 1.1 itojun arc4_getbyte(as)
145 1.1 itojun struct arc4_stream *as;
146 1.1 itojun {
147 1.1 itojun u_int8_t si, sj;
148 1.1 itojun
149 1.1 itojun as->i = (as->i + 1);
150 1.1 itojun si = as->s[as->i];
151 1.1 itojun as->j = (as->j + si);
152 1.1 itojun sj = as->s[as->j];
153 1.1 itojun as->s[as->i] = sj;
154 1.1 itojun as->s[as->j] = si;
155 1.1 itojun return (as->s[(si + sj) & 0xff]);
156 1.1 itojun }
157 1.1 itojun
158 1.1 itojun static inline u_int32_t
159 1.1 itojun arc4_getword(as)
160 1.1 itojun struct arc4_stream *as;
161 1.1 itojun {
162 1.1 itojun u_int32_t val;
163 1.1 itojun val = arc4_getbyte(as) << 24;
164 1.1 itojun val |= arc4_getbyte(as) << 16;
165 1.1 itojun val |= arc4_getbyte(as) << 8;
166 1.1 itojun val |= arc4_getbyte(as);
167 1.1 itojun return val;
168 1.1 itojun }
169 1.1 itojun
170 1.1 itojun void
171 1.1 itojun arc4random_stir()
172 1.1 itojun {
173 1.1 itojun if (!rs_initialized) {
174 1.1 itojun arc4_init(&rs);
175 1.1 itojun rs_initialized = 1;
176 1.1 itojun }
177 1.1 itojun arc4_stir(&rs);
178 1.1 itojun }
179 1.1 itojun
180 1.1 itojun void
181 1.1 itojun arc4random_addrandom(dat, datlen)
182 1.1 itojun u_char *dat;
183 1.1 itojun int datlen;
184 1.1 itojun {
185 1.1 itojun if (!rs_initialized)
186 1.1 itojun arc4random_stir();
187 1.1 itojun arc4_addrandom(&rs, dat, datlen);
188 1.1 itojun }
189 1.1 itojun
190 1.1 itojun u_int32_t
191 1.1 itojun arc4random()
192 1.1 itojun {
193 1.1 itojun if (!rs_initialized)
194 1.1 itojun arc4random_stir();
195 1.1 itojun return arc4_getword(&rs);
196 1.1 itojun }
197 1.1 itojun
198 1.1 itojun #if 0
199 1.1 itojun /*-------- Test code for i386 --------*/
200 1.1 itojun #include <stdio.h>
201 1.1 itojun #include <machine/pctr.h>
202 1.1 itojun int
203 1.1 itojun main(int argc, char **argv)
204 1.1 itojun {
205 1.1 itojun const int iter = 1000000;
206 1.1 itojun int i;
207 1.1 itojun pctrval v;
208 1.1 itojun
209 1.1 itojun v = rdtsc();
210 1.1 itojun for (i = 0; i < iter; i++)
211 1.1 itojun arc4random();
212 1.1 itojun v = rdtsc() - v;
213 1.1 itojun v /= iter;
214 1.1 itojun
215 1.1 itojun printf("%qd cycles\n", v);
216 1.1 itojun }
217 1.1 itojun #endif
218