hxrand.cpp
上传用户:dangjiwu
上传日期:2013-07-19
资源大小:42019k
文件大小:7k
源码类别:

Symbian

开发平台:

Visual C++

  1. /* ***** BEGIN LICENSE BLOCK *****
  2.  * Source last modified: $Id: hxrand.cpp,v 1.2.36.3 2004/07/09 01:48:15 hubbe Exp $
  3.  * 
  4.  * Portions Copyright (c) 1995-2004 RealNetworks, Inc. All Rights Reserved.
  5.  * 
  6.  * The contents of this file, and the files included with this file,
  7.  * are subject to the current version of the RealNetworks Public
  8.  * Source License (the "RPSL") available at
  9.  * http://www.helixcommunity.org/content/rpsl unless you have licensed
  10.  * the file under the current version of the RealNetworks Community
  11.  * Source License (the "RCSL") available at
  12.  * http://www.helixcommunity.org/content/rcsl, in which case the RCSL
  13.  * will apply. You may also obtain the license terms directly from
  14.  * RealNetworks.  You may not use this file except in compliance with
  15.  * the RPSL or, if you have a valid RCSL with RealNetworks applicable
  16.  * to this file, the RCSL.  Please see the applicable RPSL or RCSL for
  17.  * the rights, obligations and limitations governing use of the
  18.  * contents of the file.
  19.  * 
  20.  * Alternatively, the contents of this file may be used under the
  21.  * terms of the GNU General Public License Version 2 or later (the
  22.  * "GPL") in which case the provisions of the GPL are applicable
  23.  * instead of those above. If you wish to allow use of your version of
  24.  * this file only under the terms of the GPL, and not to allow others
  25.  * to use your version of this file under the terms of either the RPSL
  26.  * or RCSL, indicate your decision by deleting the provisions above
  27.  * and replace them with the notice and other provisions required by
  28.  * the GPL. If you do not delete the provisions above, a recipient may
  29.  * use your version of this file under the terms of any one of the
  30.  * RPSL, the RCSL or the GPL.
  31.  * 
  32.  * This file is part of the Helix DNA Technology. RealNetworks is the
  33.  * developer and/or licensor of the Original Code and owns the
  34.  * copyrights in the portions it created.
  35.  * 
  36.  * This file, and the files included with this file, is distributed
  37.  * and made available on an 'AS IS' basis, WITHOUT WARRANTY OF ANY
  38.  * KIND, EITHER EXPRESS OR IMPLIED, AND REALNETWORKS HEREBY DISCLAIMS
  39.  * ALL SUCH WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES
  40.  * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, QUIET
  41.  * ENJOYMENT OR NON-INFRINGEMENT.
  42.  * 
  43.  * Technology Compatibility Kit Test Suite(s) Location:
  44.  *    http://www.helixcommunity.org/content/tck
  45.  * 
  46.  * Contributor(s):
  47.  * 
  48.  * ***** END LICENSE BLOCK ***** */
  49. /*
  50.  * 
  51.  * (c) Copyright 1989 OPEN SOFTWARE FOUNDATION, INC.
  52.  * (c) Copyright 1989 HEWLETT-PACKARD COMPANY
  53.  * (c) Copyright 1989 DIGITAL EQUIPMENT CORPORATION
  54.  * To anyone who acknowledges that this file is provided "AS IS"
  55.  * without any express or implied warranty:
  56.  *                 permission to use, copy, modify, and distribute this
  57.  * file for any purpose is hereby granted without fee, provided that
  58.  * the above copyright notices and this notice appears in all source
  59.  * code copies, and that none of the names of Open Software
  60.  * Foundation, Inc., Hewlett-Packard Company, or Digital Equipment
  61.  * Corporation be used in advertising or publicity pertaining to
  62.  * distribution of the software without specific, written prior
  63.  * permission.  Neither Open Software Foundation, Inc., Hewlett-
  64.  * Packard Company, nor Digital Equipment Corporation makes any
  65.  * representations about the suitability of this software for any
  66.  * purpose.
  67.  * 
  68.  */
  69. /****************************************************************************
  70. **
  71. **    U U I D   T R U E   R A N D O M   N U M B E R   G E N E R A T O R
  72. **
  73. *****************************************************************************
  74. **
  75. ** This random number generator (RNG) was found in the ALGORITHMS Notesfile.
  76. **
  77. ** (Note 16.7, July 7, 1989 by Robert (RDVAX::)Gries, Cambridge Research Lab,
  78. **  Computational Quality Group)
  79. **
  80. ** It is really a "Multiple Prime Random Number Generator" (MPRNG) and is
  81. ** completely discussed in reference #1 (see below).
  82. **
  83. **   References:
  84. **   1) "The Multiple Prime Random Number Generator" by Alexander Hass
  85. **      pp. 368 to 381 in ACM Transactions on Mathematical Software,
  86. **      December, 1987
  87. **   2) "The Art of Computer Programming: Seminumerical Algorithms
  88. **      (vol 2)" by Donald E. Knuth, pp. 39 to 113.
  89. **
  90. ** A summary of the notesfile entry follows:
  91. **
  92. ** Gries discusses the two RNG's available for ULTRIX-C.  The default RNG
  93. ** uses a Linear Congruential Method (very popular) and the second RNG uses
  94. ** a technique known as a linear feedback shift register.
  95. **
  96. ** The first (default) RNG suffers from bit-cycles (patterns/repetition),
  97. ** ie. it's "not that random."
  98. **
  99. ** While the second RNG passes all the emperical tests, there are "states"
  100. ** that become "stable", albeit contrived.
  101. **
  102. ** Gries then presents the MPRNG and says that it passes all emperical
  103. ** tests listed in reference #2.  In addition, the number of calls to the
  104. ** MPRNG before a sequence of bit position repeats appears to have a normal
  105. ** distribution.
  106. **
  107. ** Note (mbs): I have coded the Gries's MPRNG with the same constants that
  108. ** he used in his paper.  I have no way of knowing whether they are "ideal"
  109. ** for the range of numbers we are dealing with.
  110. **
  111. ****************************************************************************/
  112. #include "hxtypes.h"
  113. #include "hxrand.h"
  114. #include "hxtick.h"
  115. #ifdef _WINDOWS
  116. #include <windows.h>
  117. #endif
  118. #include "hxheap.h"
  119. #ifdef _DEBUG
  120. #undef HX_THIS_FILE
  121. static const char HX_THIS_FILE[] = __FILE__;
  122. #endif
  123. /*
  124.  * optimal/recommended starting values according to the reference
  125.  */
  126. #define RAND_M_INIT     971;      
  127. #define RAND_IA_INIT     11113;    
  128. #define RAND_IB_INIT     104322;   
  129. #define RAND_IRAND_INIT     4181;
  130. CMultiplePrimeRandom::CMultiplePrimeRandom()
  131. {
  132.     m_rand_m     = RAND_M_INIT;      
  133.     m_rand_ia     = RAND_IA_INIT;    
  134.     m_rand_ib     = RAND_IB_INIT;   
  135.     m_rand_irand    = RAND_IRAND_INIT;
  136.     SetSeed(HX_GET_TICKCOUNT());
  137. }
  138. CMultiplePrimeRandom::CMultiplePrimeRandom(ULONG32 seed)
  139. {
  140.     m_rand_m     = RAND_M_INIT;      
  141.     m_rand_ia     = RAND_IA_INIT;    
  142.     m_rand_ib     = RAND_IB_INIT;   
  143.     m_rand_irand    = RAND_IRAND_INIT;
  144.     SetSeed(seed);
  145. }
  146. void CMultiplePrimeRandom::SetSeed(ULONG32 seed)
  147. {
  148.     // seed the generator
  149.     m_rand_irand += seed; 
  150. }
  151. ULONG32 CMultiplePrimeRandom::GetRandomNumber(void)
  152. {
  153.     m_rand_m += 7;
  154.     m_rand_ia += 1907;
  155.     m_rand_ib += 73939;
  156.     if (m_rand_m >= 9973) m_rand_m -= 9871;
  157.     if (m_rand_ia >= 99991) m_rand_ia -= 89989;
  158.     if (m_rand_ib >= 224729) m_rand_ib -= 96233;
  159.     m_rand_irand = (m_rand_irand * m_rand_m) + m_rand_ia + m_rand_ib;
  160.     return m_rand_irand;
  161. }