md_rand.c 15.9 KB
Newer Older
1
/* crypto/rand/md_rand.c */
2
/* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com)
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 48 49 50 51 52 53 54 55 56 57 58
 * All rights reserved.
 *
 * This package is an SSL implementation written
 * by Eric Young (eay@cryptsoft.com).
 * The implementation was written so as to conform with Netscapes SSL.
 * 
 * This library is free for commercial and non-commercial use as long as
 * the following conditions are aheared to.  The following conditions
 * apply to all code found in this distribution, be it the RC4, RSA,
 * lhash, DES, etc., code; not just the SSL code.  The SSL documentation
 * included with this distribution is covered by the same copyright terms
 * except that the holder is Tim Hudson (tjh@cryptsoft.com).
 * 
 * Copyright remains Eric Young's, and as such any Copyright notices in
 * the code are not to be removed.
 * If this package is used in a product, Eric Young should be given attribution
 * as the author of the parts of the library used.
 * This can be in the form of a textual message at program startup or
 * in documentation (online or textual) provided with the package.
 * 
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 * 1. Redistributions of source code must retain the copyright
 *    notice, this list of conditions and the following disclaimer.
 * 2. Redistributions in binary form must reproduce the above copyright
 *    notice, this list of conditions and the following disclaimer in the
 *    documentation and/or other materials provided with the distribution.
 * 3. All advertising materials mentioning features or use of this software
 *    must display the following acknowledgement:
 *    "This product includes cryptographic software written by
 *     Eric Young (eay@cryptsoft.com)"
 *    The word 'cryptographic' can be left out if the rouines from the library
 *    being used are not cryptographic related :-).
 * 4. If you include any Windows specific code (or a derivative thereof) from 
 *    the apps directory (application code) you must include an acknowledgement:
 *    "This product includes software written by Tim Hudson (tjh@cryptsoft.com)"
 * 
 * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND
 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
 * SUCH DAMAGE.
 * 
 * The licence and distribution terms for any publically available version or
 * derivative of this code cannot be changed.  i.e. this code cannot simply be
 * copied and put under another distribution licence
 * [including the GNU Public Licence.]
 */

59 60
#define ENTROPY_NEEDED 32  /* require 128 bits of randomness */

B
Bodo Möller 已提交
61 62 63 64 65 66 67
#ifndef MD_RAND_DEBUG
# ifndef NDEBUG
#   define NDEBUG
# endif
#endif

#include <assert.h>
68 69
#include <stdio.h>
#include <time.h>
B
Ben Laurie 已提交
70
#include <string.h>
71

72
#include "openssl/e_os.h"
73

74
#include <openssl/crypto.h>
75
#include <openssl/err.h>
76 77

#if !defined(USE_MD5_RAND) && !defined(USE_SHA1_RAND) && !defined(USE_MDC2_RAND) && !defined(USE_MD2_RAND)
U
Ulf Möller 已提交
78
#if !defined(NO_SHA) && !defined(NO_SHA1)
79
#define USE_SHA1_RAND
U
Ulf Möller 已提交
80 81 82
#elif !defined(NO_MD5)
#define USE_MD5_RAND
#elif !defined(NO_MDC2) && !defined(NO_DES)
83 84 85 86
#define USE_MDC2_RAND
#elif !defined(NO_MD2)
#define USE_MD2_RAND
#else
U
Ulf Möller 已提交
87
#error No message digest algorithm available
88 89 90 91 92 93 94 95 96 97
#endif
#endif

/* Changed how the state buffer used.  I now attempt to 'wrap' such
 * that I don't run over the same locations the next time  go through
 * the 1023 bytes - many thanks to
 * Robert J. LeBlanc <rjl@renaissoft.com> for his comments
 */

#if defined(USE_MD5_RAND)
98
#include <openssl/md5.h>
99 100 101 102 103
#define MD_DIGEST_LENGTH	MD5_DIGEST_LENGTH
#define MD_CTX			MD5_CTX
#define MD_Init(a)		MD5_Init(a)
#define MD_Update(a,b,c)	MD5_Update(a,b,c)
#define	MD_Final(a,b)		MD5_Final(a,b)
104
#define	MD(a,b,c)		MD5(a,b,c)
105
#elif defined(USE_SHA1_RAND)
106
#include <openssl/sha.h>
107 108 109 110 111
#define MD_DIGEST_LENGTH	SHA_DIGEST_LENGTH
#define MD_CTX			SHA_CTX
#define MD_Init(a)		SHA1_Init(a)
#define MD_Update(a,b,c)	SHA1_Update(a,b,c)
#define	MD_Final(a,b)		SHA1_Final(a,b)
112
#define	MD(a,b,c)		SHA1(a,b,c)
113
#elif defined(USE_MDC2_RAND)
114
#include <openssl/mdc2.h>
115 116 117 118 119
#define MD_DIGEST_LENGTH	MDC2_DIGEST_LENGTH
#define MD_CTX			MDC2_CTX
#define MD_Init(a)		MDC2_Init(a)
#define MD_Update(a,b,c)	MDC2_Update(a,b,c)
#define	MD_Final(a,b)		MDC2_Final(a,b)
120
#define	MD(a,b,c)		MDC2(a,b,c)
121
#elif defined(USE_MD2_RAND)
122
#include <openssl/md2.h>
123 124 125 126 127
#define MD_DIGEST_LENGTH	MD2_DIGEST_LENGTH
#define MD_CTX			MD2_CTX
#define MD_Init(a)		MD2_Init(a)
#define MD_Update(a,b,c)	MD2_Update(a,b,c)
#define	MD_Final(a,b)		MD2_Final(a,b)
128
#define	MD(a,b,c)		MD2(a,b,c)
129 130
#endif

131
#include <openssl/rand.h>
132

133 134
/* #define NORAND	1 */
/* #define PREDICT	1 */
135 136 137

#define STATE_SIZE	1023
static int state_num=0,state_index=0;
138
static unsigned char state[STATE_SIZE+MD_DIGEST_LENGTH];
139
static unsigned char md[MD_DIGEST_LENGTH];
140
static long md_count[2]={0,0};
141
static int entropy=0;
142

B
Ben Laurie 已提交
143
const char *RAND_version="RAND" OPENSSL_VERSION_PTEXT;
144

145
static void ssleay_rand_cleanup(void);
146
static void ssleay_rand_seed(const void *buf, int num);
B
Bodo Möller 已提交
147
static void ssleay_rand_add(const void *buf, int num, int add_entropy);
148
static int ssleay_rand_bytes(unsigned char *buf, int num);
149

150
RAND_METHOD rand_ssleay_meth={
151 152 153
	ssleay_rand_seed,
	ssleay_rand_bytes,
	ssleay_rand_cleanup,
154
	ssleay_rand_add,
155 156
	}; 

U
Ulf Möller 已提交
157
RAND_METHOD *RAND_SSLeay(void)
158
	{
159
	return(&rand_ssleay_meth);
160 161
	}

U
Ulf Möller 已提交
162
static void ssleay_rand_cleanup(void)
163
	{
164
	memset(state,0,sizeof(state));
165 166 167
	state_num=0;
	state_index=0;
	memset(md,0,MD_DIGEST_LENGTH);
168 169
	md_count[0]=0;
	md_count[1]=0;
170
	entropy=0;
171 172
	}

173
static void ssleay_rand_add(const void *buf, int num, int add)
174
	{
B
Bodo Möller 已提交
175 176 177
	int i,j,k,st_idx;
	long md_c[2];
	unsigned char local_md[MD_DIGEST_LENGTH];
178 179 180 181 182 183
	MD_CTX m;

#ifdef NORAND
	return;
#endif

B
Bodo Möller 已提交
184 185 186 187 188 189 190 191 192 193 194 195 196 197
	/*
	 * (Based on doc/ssleay.txt, section rand.doc:)
	 *
	 * The input is chopped up into units of 16 bytes (or less for
	 * the last block).  Each of these blocks is run through the MD5
	 * message digest as follow:  The data passed to the MD5 digest
	 * is the current 'md', the same number of bytes from the 'state'
	 * (the location determined by in incremented looping index) as
	 * the current 'block', the new key data 'block', and 'count'
	 * (which is incremented after each use).
	 * The result of this is kept in 'md' and also xored into the
	 * 'state' at the same locations that were used as input into the MD5.
	 */

198 199 200
	CRYPTO_w_lock(CRYPTO_LOCK_RAND);
	st_idx=state_index;

B
Bodo Möller 已提交
201 202 203 204 205 206 207 208 209 210 211
	/* use our own copies of the counters so that even
	 * if a concurrent thread seeds with exactly the
	 * same data and uses the same subarray there's _some_
	 * difference */
	md_c[0] = md_count[0];
	md_c[1] = md_count[1];

	memcpy(local_md, md, sizeof md);

	/* state_index <= state_num <= STATE_SIZE */
	state_index += num;
212
	if (state_index >= STATE_SIZE)
213 214 215 216 217 218 219 220 221
		{
		state_index%=STATE_SIZE;
		state_num=STATE_SIZE;
		}
	else if (state_num < STATE_SIZE)	
		{
		if (state_index > state_num)
			state_num=state_index;
		}
B
Bodo Möller 已提交
222 223 224 225 226 227 228 229
	/* state_index <= state_num <= STATE_SIZE */

	/* state[st_idx], ..., state[(st_idx + num - 1) % STATE_SIZE]
	 * are what we will use now, but other threads may use them
	 * as well */

	md_count[1] += (num / MD_DIGEST_LENGTH) + (num % MD_DIGEST_LENGTH > 0);

230 231 232 233 234 235 236 237
	CRYPTO_w_unlock(CRYPTO_LOCK_RAND);

	for (i=0; i<num; i+=MD_DIGEST_LENGTH)
		{
		j=(num-i);
		j=(j > MD_DIGEST_LENGTH)?MD_DIGEST_LENGTH:j;

		MD_Init(&m);
B
Bodo Möller 已提交
238
		MD_Update(&m,local_md,MD_DIGEST_LENGTH);
239 240 241 242 243 244 245 246 247 248
		k=(st_idx+j)-STATE_SIZE;
		if (k > 0)
			{
			MD_Update(&m,&(state[st_idx]),j-k);
			MD_Update(&m,&(state[0]),k);
			}
		else
			MD_Update(&m,&(state[st_idx]),j);
			
		MD_Update(&m,buf,j);
B
Bodo Möller 已提交
249 250 251
		MD_Update(&m,(unsigned char *)&(md_c[0]),sizeof(md_c));
		MD_Final(local_md,&m);
		md_c[1]++;
252

B
Ben Laurie 已提交
253
		buf=(const char *)buf + j;
254 255 256

		for (k=0; k<j; k++)
			{
B
Bodo Möller 已提交
257 258 259 260 261 262 263 264 265
			/* Parallel threads may interfere with this,
			 * but always each byte of the new state is
			 * the XOR of some previous value of its
			 * and local_md (itermediate values may be lost).
			 * Alway using locking could hurt performance more
			 * than necessary given that conflicts occur only
			 * when the total seeding is longer than the random
			 * state. */
			state[st_idx++]^=local_md[k];
266 267 268 269 270
			if (st_idx >= STATE_SIZE)
				st_idx=0;
			}
		}
	memset((char *)&m,0,sizeof(m));
B
Bodo Möller 已提交
271 272 273 274 275 276 277 278 279 280 281 282 283 284 285

	CRYPTO_w_lock(CRYPTO_LOCK_RAND);
	/* Don't just copy back local_md into md -- this could mean that
	 * other thread's seeding remains without effect (except for
	 * the incremented counter).  By XORing it we keep at least as
	 * much entropy as fits into md. */
	for (k = 0; k < sizeof md; k++)
		{
		md[k] ^= local_md[k];
		}
	CRYPTO_w_unlock(CRYPTO_LOCK_RAND);
	
#ifndef THREADS	
	assert(md_c[1] == md_count[1]);
#endif
286 287 288 289 290 291
	entropy += add;
	}

static void ssleay_rand_seed(const void *buf, int num)
	{
	ssleay_rand_add(buf, num, num);
292 293
	}

294
static int ssleay_rand_bytes(unsigned char *buf, int num)
295 296
	{
	int i,j,k,st_num,st_idx;
297
	int ok;
B
Bodo Möller 已提交
298 299
	long md_c[2];
	unsigned char local_md[MD_DIGEST_LENGTH];
300 301 302
	MD_CTX m;
	static int init=1;
	unsigned long l;
A
Andy Polyakov 已提交
303
#ifndef GETPID_IS_MEANINGLESS
304
	pid_t curr_pid = getpid();
B
Bodo Möller 已提交
305
#endif
306 307 308 309 310 311 312 313 314 315
#ifdef DEVRANDOM
	FILE *fh;
#endif

#ifdef PREDICT
	{
	static unsigned char val=0;

	for (i=0; i<num; i++)
		buf[i]=val++;
316
	return(1);
317 318 319
	}
#endif

B
Bodo Möller 已提交
320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335
	/*
	 * (Based on doc/ssleay.txt, section rand.doc:)
	 *
	 * For each group of 8 bytes (or less), we do the following,
	 *
	 * Input into MD5, the top 8 bytes from 'md', the byte that are
	 * to be overwritten by the random bytes and bytes from the
	 * 'state' (incrementing looping index).  From this digest output
	 * (which is kept in 'md'), the top (upto) 8 bytes are
	 * returned to the caller and the bottom (upto) 8 bytes are xored
	 * into the 'state'.
	 * Finally, after we have finished 'num' random bytes for the
	 * caller, 'count' (which is incremented) and the local and globl 'md'
	 * are fed into MD5 and the results are kept in the global 'md'.
	 */

336 337 338 339 340 341 342
	CRYPTO_w_lock(CRYPTO_LOCK_RAND);

	if (init)
		{
		CRYPTO_w_unlock(CRYPTO_LOCK_RAND);
		/* put in some default random data, we need more than
		 * just this */
343
		RAND_add(&m,sizeof(m),0);
A
Andy Polyakov 已提交
344
#ifndef GETPID_IS_MEANINGLESS
345
		l=curr_pid;
346
		RAND_add(&l,sizeof(l),0);
347
		l=getuid();
348
		RAND_add(&l,sizeof(l),0);
349 350
#endif
		l=time(NULL);
351
		RAND_add(&l,sizeof(l),0);
352

353
#ifdef DEVRANDOM
354 355 356 357 358 359 360 361 362
		/* 
		 * Use a random entropy pool device.
		 * Linux 1.3.x and FreeBSD-Current has 
		 * this. Use /dev/urandom if you can
		 * as /dev/random will block if it runs out
		 * of random entries.
		 */
		if ((fh = fopen(DEVRANDOM, "r")) != NULL)
			{
U
Ulf Möller 已提交
363
			unsigned char tmpbuf[ENTROPY_NEEDED];
B
Bodo Möller 已提交
364
			int n;
365

B
Bodo Möller 已提交
366
			n=fread((unsigned char *)tmpbuf,1,ENTROPY_NEEDED,fh);
367
			fclose(fh);
B
Bodo Möller 已提交
368 369
			RAND_add(tmpbuf,sizeof tmpbuf,n);
			memset(tmpbuf,0,n);
370
			}
371
#endif
372 373 374 375 376
#ifdef PURIFY
		memset(state,0,STATE_SIZE);
		memset(md,0,MD_DIGEST_LENGTH);
#endif
		CRYPTO_w_lock(CRYPTO_LOCK_RAND);
B
Bodo Möller 已提交
377
		init=0;
378 379
		}

380 381
	ok = (entropy >= ENTROPY_NEEDED);

382 383
	st_idx=state_index;
	st_num=state_num;
B
Bodo Möller 已提交
384 385 386 387
	md_c[0] = md_count[0];
	md_c[1] = md_count[1];
	memcpy(local_md, md, sizeof md);

388 389
	state_index+=num;
	if (state_index > state_num)
B
Bodo Möller 已提交
390 391 392 393
		state_index %= state_num;

	/* state[st_idx], ..., state[(st_idx + num - 1) % st_num]
	 * are now ours (but other threads may use them too) */
394

B
Bodo Möller 已提交
395
	md_count[0] += 1;
396 397 398 399 400 401 402
	CRYPTO_w_unlock(CRYPTO_LOCK_RAND);

	while (num > 0)
		{
		j=(num >= MD_DIGEST_LENGTH/2)?MD_DIGEST_LENGTH/2:num;
		num-=j;
		MD_Init(&m);
A
Andy Polyakov 已提交
403
#ifndef GETPID_IS_MEANINGLESS
404 405 406 407 408 409
		if (curr_pid) /* just in the first iteration to save time */
			{
			MD_Update(&m,(unsigned char*)&curr_pid,sizeof curr_pid);
			curr_pid = 0;
			}
#endif
B
Bodo Möller 已提交
410 411
		MD_Update(&m,&(local_md[MD_DIGEST_LENGTH/2]),MD_DIGEST_LENGTH/2);
		MD_Update(&m,(unsigned char *)&(md_c[0]),sizeof(md_c));
412 413 414 415 416 417 418 419 420 421 422
#ifndef PURIFY
		MD_Update(&m,buf,j); /* purify complains */
#endif
		k=(st_idx+j)-st_num;
		if (k > 0)
			{
			MD_Update(&m,&(state[st_idx]),j-k);
			MD_Update(&m,&(state[0]),k);
			}
		else
			MD_Update(&m,&(state[st_idx]),j);
B
Bodo Möller 已提交
423
		MD_Final(local_md,&m);
424 425 426

		for (i=0; i<j; i++)
			{
B
Bodo Möller 已提交
427 428
			state[st_idx++]^=local_md[i]; /* may compete with other threads */
			*(buf++)=local_md[i+MD_DIGEST_LENGTH/2];
429 430 431 432 433 434
			if (st_idx >= st_num)
				st_idx=0;
			}
		}

	MD_Init(&m);
B
Bodo Möller 已提交
435 436 437
	MD_Update(&m,(unsigned char *)&(md_c[0]),sizeof(md_c));
	MD_Update(&m,local_md,MD_DIGEST_LENGTH);
	CRYPTO_w_lock(CRYPTO_LOCK_RAND);
438 439
	MD_Update(&m,md,MD_DIGEST_LENGTH);
	MD_Final(md,&m);
B
Bodo Möller 已提交
440 441
	CRYPTO_w_unlock(CRYPTO_LOCK_RAND);

442
	memset(&m,0,sizeof(m));
443 444 445 446 447 448 449
	if (ok)
		return(1);
	else
		{
		RANDerr(RAND_F_SSLEAY_RAND_BYTES,RAND_R_PRNG_NOT_SEEDED);
		return(0);
		}
450 451 452 453
	}

#ifdef WINDOWS
#include <windows.h>
454
#include <openssl/rand.h>
455 456 457 458 459 460 461 462 463 464 465

/*****************************************************************************
 * Initialisation function for the SSL random generator.  Takes the contents
 * of the screen as random seed.
 *
 * Created 960901 by Gertjan van Oosten, gertjan@West.NL, West Consulting B.V.
 *
 * Code adapted from
 * <URL:http://www.microsoft.com/kb/developr/win_dk/q97193.htm>;
 * the original copyright message is:
 *
466 467 468 469 470 471 472
 *   (C) Copyright Microsoft Corp. 1993.  All rights reserved.
 *
 *   You have a royalty-free right to use, modify, reproduce and
 *   distribute the Sample Files (and/or any modified version) in
 *   any way you find useful, provided that you agree that
 *   Microsoft has no warranty obligations or liability for any
 *   Sample Application Files which are modified.
473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530
 */
/*
 * I have modified the loading of bytes via RAND_seed() mechanism since
 * the origional would have been very very CPU intensive since RAND_seed()
 * does an MD5 per 16 bytes of input.  The cost to digest 16 bytes is the same
 * as that to digest 56 bytes.  So under the old system, a screen of
 * 1024*768*256 would have been CPU cost of approximatly 49,000 56 byte MD5
 * digests or digesting 2.7 mbytes.  What I have put in place would
 * be 48 16k MD5 digests, or efectivly 48*16+48 MD5 bytes or 816 kbytes
 * or about 3.5 times as much.
 * - eric 
 */
void RAND_screen(void)
{
  HDC		hScrDC;		/* screen DC */
  HDC		hMemDC;		/* memory DC */
  HBITMAP	hBitmap;	/* handle for our bitmap */
  HBITMAP	hOldBitmap;	/* handle for previous bitmap */
  BITMAP	bm;		/* bitmap properties */
  unsigned int	size;		/* size of bitmap */
  char		*bmbits;	/* contents of bitmap */
  int		w;		/* screen width */
  int		h;		/* screen height */
  int		y;		/* y-coordinate of screen lines to grab */
  int		n = 16;		/* number of screen lines to grab at a time */

  /* Create a screen DC and a memory DC compatible to screen DC */
  hScrDC = CreateDC("DISPLAY", NULL, NULL, NULL);
  hMemDC = CreateCompatibleDC(hScrDC);

  /* Get screen resolution */
  w = GetDeviceCaps(hScrDC, HORZRES);
  h = GetDeviceCaps(hScrDC, VERTRES);

  /* Create a bitmap compatible with the screen DC */
  hBitmap = CreateCompatibleBitmap(hScrDC, w, n);

  /* Select new bitmap into memory DC */
  hOldBitmap = SelectObject(hMemDC, hBitmap);

  /* Get bitmap properties */
  GetObject(hBitmap, sizeof(BITMAP), (LPSTR)&bm);
  size = (unsigned int)bm.bmWidthBytes * bm.bmHeight * bm.bmPlanes;

  bmbits = Malloc(size);
  if (bmbits) {
    /* Now go through the whole screen, repeatedly grabbing n lines */
    for (y = 0; y < h-n; y += n)
    	{
	unsigned char md[MD_DIGEST_LENGTH];

	/* Bitblt screen DC to memory DC */
	BitBlt(hMemDC, 0, 0, w, n, hScrDC, 0, y, SRCCOPY);

	/* Copy bitmap bits from memory DC to bmbits */
	GetBitmapBits(hBitmap, size, bmbits);

	/* Get the MD5 of the bitmap */
531
	MD(bmbits,size,md);
532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548

	/* Seed the random generator with the MD5 digest */
	RAND_seed(md, MD_DIGEST_LENGTH);
	}

    Free(bmbits);
  }

  /* Select old bitmap back into memory DC */
  hBitmap = SelectObject(hMemDC, hOldBitmap);

  /* Clean up */
  DeleteObject(hBitmap);
  DeleteDC(hMemDC);
  DeleteDC(hScrDC);
}
#endif