3 * Secure random number generator
5 * (c) 1998 Straylight/Edgeware
8 /*----- Licensing notice --------------------------------------------------*
10 * This file is part of Catacomb.
12 * Catacomb is free software; you can redistribute it and/or modify
13 * it under the terms of the GNU Library General Public License as
14 * published by the Free Software Foundation; either version 2 of the
15 * License, or (at your option) any later version.
17 * Catacomb is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 * GNU Library General Public License for more details.
22 * You should have received a copy of the GNU Library General Public
23 * License along with Catacomb; if not, write to the Free
24 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
28 /*----- Header files ------------------------------------------------------*/
36 #include <mLib/bits.h>
48 #include "twofish-counter.h"
51 #define CIPHER_CTX twofish_counterctx
52 #define CIPHER_INIT twofish_counterinit
53 #define CIPHER_ENCRYPT twofish_counterencrypt
54 #define CIPHER_IVSZ TWOFISH_BLKSZ
55 #define CIPHER_KEYSZ TWOFISH_KEYSZ
57 #define HASH_CTX sha256_ctx
58 #define HASH_INIT sha256_init
59 #define HASH sha256_hash
60 #define HASH_DONE sha256_done
61 #define HASH_SZ SHA256_HASHSZ
63 /*----- Static variables --------------------------------------------------*/
65 static const grand_ops gops;
67 typedef struct rand__gctx {
76 { "Catacomb global random byte pool" },
80 /*----- Macros ------------------------------------------------------------*/
82 #define RAND_RESOLVE(r) \
83 do { if ((r) == RAND_GLOBAL) r = &rand_global.p; } while (0)
85 #define GENCHECK(r) do { \
86 unsigned gen = rand_generation(); \
87 if (r->gen != gen) { r->gen = gen; rand_gate(r); } \
90 static int quick(rand_pool *);
91 #define QUICK(r) do { \
93 if ((r)->s && (r)->s->timer) (r)->s->timer(r); \
96 /*----- Main code ---------------------------------------------------------*/
98 /* --- @rand_init@ --- *
100 * Arguments: @rand_pool *r@ = pointer to a randomness pool
104 * Use: Initializes a randomness pool. The pool doesn't start out
105 * very random: that's your job to sort out. A good suggestion
106 * would be to attach an appropriate noise source and call
110 void rand_init(rand_pool *r)
113 memset(r->pool, 0, sizeof(r->pool));
114 memset(r->buf, 0, sizeof(r->buf));
115 r->gen = rand_generation();
118 r->ibits = r->obits = 0;
120 r->s = &noise_source;
125 /* --- @rand_noisesrc@ --- *
127 * Arguments: @rand_pool *r@ = pointer to a randomness pool
128 * @const rand_source *s@ = pointer to source definition
132 * Use: Sets a noise source for a randomness pool. When the pool's
133 * estimate of good random bits falls to zero, the @getnoise@
134 * function is called, passing the pool handle as an argument.
135 * It is expected to increase the number of good bits by at
136 * least one, because it'll be called over and over again until
137 * there are enough bits to satisfy the caller. The @timer@
138 * function is called frequently throughout the generator's
142 void rand_noisesrc(rand_pool *r, const rand_source *s)
148 /* --- @rand_quick@ --- *
150 * Arguments: @rand_pool *r@ = pointer to a randomness pool
152 * Returns: Zero on success; @-1@ on failure.
154 * Use Attempts to use some machine-specific `quick' source of
155 * entropy to top up @r@. This may not do anything at all on
159 CPU_DISPATCH(static, return, int, quick, (rand_pool *r), (r),
160 pick_quick, trivial_quick);
162 static int trivial_quick(rand_pool *r) { return (-1); }
164 #if __GNUC__ && (CPUFAM_X86 || CPUFAM_AMD64)
165 static int rdrand_quick(rand_pool *r)
170 __asm__ ("0: rdrand %0; jc 9f; dec %1; jnz 0b; 9:"
171 : "=r" (rr), "=r" (i) : "1" (i) : "cc");
173 rand_add(r, &rr, sizeof(rr), 8*sizeof(rr));
178 static quick__functype *pick_quick(void)
180 #if __GNUC__ && (CPUFAM_X86 || CPUFAM_AMD64)
181 DISPATCH_PICK_COND(rand_quick, rdrand_quick,
182 cpu_feature_p(CPUFEAT_X86_RDRAND));
184 DISPATCH_PICK_FALLBACK(rand_quick, trivial_quick);
187 int rand_quick(rand_pool *r) { RAND_RESOLVE(r); return (quick(r)); }
189 /* --- @rand_seed@ --- *
191 * Arguments: @rand_pool *r@ = pointer to a randomness pool
192 * @unsigned bits@ = number of bits to ensure
196 * Use: Ensures that there are at least @bits@ good bits of entropy
197 * in the pool. It is recommended that you call this after
198 * initializing a new pool. Requesting @bits > RAND_IBITS@ is
199 * doomed to failure (and is an error).
202 void rand_seed(rand_pool *r, unsigned bits)
206 assert(((void)"bits pointlessly large in rand_seed", bits <= RAND_IBITS));
207 assert(((void)"no noise source in rand_seed", r->s));
209 while (r->ibits < bits)
214 /* --- @rand_key@ --- *
216 * Arguments: @rand_pool *r@ = pointer to a randomness pool
217 * @const void *k@ = pointer to key data
218 * @size_t sz@ = size of key data
222 * Use: Sets the secret key for a randomness pool. The key is used
223 * when mixing in new random bits.
226 void rand_key(rand_pool *r, const void *k, size_t sz)
230 static const char label[] = "Catacomb random pool key";
234 assert(HASH_SZ >= RAND_KEYSZ);
236 HASH(&hc, label, sizeof(label));
237 if (sz) HASH(&hc, k, sz);
239 memcpy(r->k.k, h, RAND_KEYSZ);
242 /* --- @rand_add@ --- *
244 * Arguments: @rand_pool *r@ = pointer to a randomness pool
245 * @const void *p@ = pointer a buffer of data to add
246 * @size_t sz@ = size of the data buffer
247 * @unsigned goodbits@ = number of good bits estimated in buffer
251 * Use: Mixes the data in the buffer with the contents of the
252 * pool. The estimate of the number of good bits is added to
253 * the pool's own count. The mixing operation is not
254 * cryptographically strong. However, data in the input pool
255 * isn't output directly, only through the one-way gating
256 * operation, so that shouldn't matter.
259 void rand_add(rand_pool *r, const void *p, size_t sz, unsigned goodbits)
264 #if RAND_POOLSZ != 128
265 # error Polynomial in rand_add is out of date. Fix it.
270 i = r->i; rot = r->irot;
274 r->pool[i] ^= (ROL8(o, rot) ^
275 r->pool[(i + 1) % RAND_POOLSZ] ^
276 r->pool[(i + 2) % RAND_POOLSZ] ^
277 r->pool[(i + 7) % RAND_POOLSZ]);
279 i++; if (i >= RAND_POOLSZ) i -= RAND_POOLSZ;
285 r->ibits += goodbits;
286 if (r->ibits > RAND_IBITS)
287 r->ibits = RAND_IBITS;
290 /* --- @rand_goodbits@ --- *
292 * Arguments: @rand_pool *r@ = pointer to a randomness pool
294 * Returns: Estimate of the number of good bits remaining in the pool.
297 unsigned rand_goodbits(rand_pool *r)
300 return (r->ibits + r->obits);
303 /* --- @rand_gate@ --- *
305 * Arguments: @rand_pool *r@ = pointer to a randomness pool
309 * Use: Mixes up the entire state of the generator in a nonreversible
313 void rand_gate(rand_pool *r)
315 octet h[HASH_SZ], g[4];
322 /* --- Hash up all the data in the pool --- */
325 STORE32(g, r->gen); HASH(&hc, g, sizeof(g));
326 HASH(&hc, r->pool, RAND_POOLSZ);
327 HASH(&hc, r->buf, RAND_BUFSZ);
331 /* --- Now mangle all of the data based on the hash --- */
333 assert(CIPHER_KEYSZ <= HASH_SZ);
334 CIPHER_INIT(&cc, h, CIPHER_KEYSZ, 0);
335 CIPHER_ENCRYPT(&cc, r->pool, r->pool, RAND_POOLSZ);
336 CIPHER_ENCRYPT(&cc, r->buf, r->buf, RAND_BUFSZ);
339 /* --- Reset the various state variables --- */
342 r->obits += r->ibits;
343 if (r->obits > RAND_OBITS) {
344 r->ibits = r->obits - r->ibits;
345 r->obits = RAND_OBITS;
351 /* --- @rand_stretch@ --- *
353 * Arguments: @rand_pool *r@ = pointer to a randomness pool
357 * Use: Stretches the contents of the output buffer by transforming
358 * it in a nonreversible way. This doesn't add any entropy
359 * worth speaking about, but it works well enough when the
360 * caller doesn't care about that sort of thing.
363 void rand_stretch(rand_pool *r)
365 octet h[HASH_SZ], g[4];
372 /* --- Hash up all the data in the buffer --- */
375 STORE32(g, r->gen); HASH(&hc, g, sizeof(g));
376 HASH(&hc, r->pool, RAND_POOLSZ);
377 HASH(&hc, r->buf, RAND_BUFSZ);
381 /* --- Now mangle the buffer based on the hash --- */
383 assert(CIPHER_KEYSZ <= HASH_SZ);
384 CIPHER_INIT(&cc, h, CIPHER_KEYSZ, 0);
385 CIPHER_ENCRYPT(&cc, r->buf, r->buf, RAND_BUFSZ);
388 /* --- Reset the various state variables --- */
394 /* --- @rand_get@ --- *
396 * Arguments: @rand_pool *r@ = pointer to a randomness pool
397 * @void *p@ = pointer to output buffer
398 * @size_t sz@ = size of output buffer
402 * Use: Gets random data from the pool. The pool's contents can't be
403 * determined from the output of this function; nor can the
404 * output data be determined from a knowledge of the data input
405 * to the pool wihtout also having knowledge of the secret key.
406 * The good bits counter is decremented, although no special
407 * action is taken if it reaches zero.
410 void rand_get(rand_pool *r, void *p, size_t sz)
421 if (r->o + sz <= RAND_BUFSZ) {
422 memcpy(o, r->buf + r->o, sz);
426 size_t chunk = RAND_BUFSZ - r->o;
428 memcpy(o, r->buf + r->o, chunk);
436 if (r->obits > sz * 8)
442 /* --- @rand_getgood@ --- *
444 * Arguments: @rand_pool *r@ = pointer to a randomness pool
445 * @void *p@ = pointer to output buffer
446 * @size_t sz@ = size of output buffer
450 * Use: Gets random data from the pool, ensuring that there are
451 * enough good bits. This interface isn't recommended: it makes
452 * the generator slow, and doesn't provide much more security
453 * than @rand_get@, assuming you've previously done a
457 void rand_getgood(rand_pool *r, void *p, size_t sz)
465 if (!r->s || !r->s->getnoise) {
475 if (chunk * 8 > r->obits) {
476 if (chunk * 8 > r->ibits + r->obits)
477 do r->s->getnoise(r); while (r->ibits + r->obits < 256);
479 if (chunk * 8 > r->obits)
480 chunk = r->obits / 8;
483 if (chunk + r->o > RAND_BUFSZ)
484 chunk = RAND_BUFSZ - r->o;
486 memcpy(o, r->buf + r->o, chunk);
488 r->obits -= chunk * 8;
494 /*----- Generic random number generator interface -------------------------*/
496 static void gdestroy(grand *r)
499 if (g != &rand_global) {
505 static int gmisc(grand *r, unsigned op, ...)
514 switch (va_arg(ap, unsigned)) {
517 case GRAND_SEEDUINT32:
518 case GRAND_SEEDBLOCK:
535 case GRAND_SEEDINT: {
536 unsigned u = va_arg(ap, unsigned);
537 rand_add(&g->p, &u, sizeof(u), sizeof(u));
539 case GRAND_SEEDUINT32: {
540 uint32 i = va_arg(ap, uint32);
541 rand_add(&g->p, &i, sizeof(i), 4);
543 case GRAND_SEEDBLOCK: {
544 const void *p = va_arg(ap, const void *);
545 size_t sz = va_arg(ap, size_t);
546 rand_add(&g->p, p, sz, sz);
548 case GRAND_SEEDRAND: {
549 grand *rr = va_arg(ap, grand *);
551 rr->ops->fill(rr, buf, sizeof(buf));
552 rand_add(&g->p, buf, sizeof(buf), 8);
561 const void *k = va_arg(ap, const void *);
562 size_t sz = va_arg(ap, size_t);
563 rand_key(&g->p, k, sz);
566 rand_noisesrc(&g->p, va_arg(ap, const rand_source *));
569 rand_seed(&g->p, va_arg(ap, unsigned));
575 rc = rand_goodbits(&g->p);
578 const void *p = va_arg(ap, const void *);
579 size_t sz = va_arg(ap, size_t);
580 unsigned goodbits = va_arg(ap, unsigned);
581 rand_add(&g->p, p, sz, goodbits);
592 static octet gbyte(grand *r)
596 rand_getgood(&g->p, &o, 1);
600 static uint32 gword(grand *r)
604 rand_getgood(&g->p, &b, sizeof(b));
608 static void gfill(grand *r, void *p, size_t sz)
611 rand_get(&g->p, p, sz);
614 static const grand_ops gops = {
618 gword, gbyte, gword, grand_defaultrange, gfill
621 /* --- @rand_create@ --- *
625 * Returns: Pointer to a generic generator.
627 * Use: Constructs a generic generator interface over a Catacomb
628 * entropy pool generator.
631 grand *rand_create(void)
633 gctx *g = S_CREATE(gctx);
639 /*----- That's all, folks -------------------------------------------------*/