home *** CD-ROM | disk | FTP | other *** search
/ GEMini Atari / GEMini_Atari_CD-ROM_Walnut_Creek_December_1993.iso / files / gnu / g__lib / random.cc < prev    next >
C/C++ Source or Header  |  1993-07-23  |  4KB  |  173 lines
  1. /* 
  2. Copyright (C) 1988 Free Software Foundation
  3.     written by Dirk Grunwald (grunwald@cs.uiuc.edu)
  4.  
  5. This file is part of GNU CC.
  6.  
  7. GNU CC is distributed in the hope that it will be useful,
  8. but WITHOUT ANY WARRANTY.  No author or distributor
  9. accepts responsibility to anyone for the consequences of using it
  10. or for whether it serves any particular purpose or works at all,
  11. unless he says so in writing.  Refer to the GNU CC General Public
  12. License for full details.
  13.  
  14. Everyone is granted permission to copy, modify and redistribute
  15. GNU CC, but only under the conditions described in the
  16. GNU CC General Public License.   A copy of this license is
  17. supposed to have been given to you along with GNU CC so you
  18. can know your rights and responsibilities.  It should be in a
  19. file named COPYING.  Among other things, the copyright notice
  20. and this notice must be preserved on all copies.  
  21. */
  22. #include <builtin.h>
  23. #include <Random.h>
  24.  
  25. double Random::operator()()
  26. {
  27.     (*lib_error_handler)("Random", "Unimplemented");
  28. }
  29.  
  30. #include <Binomial.h>
  31.  
  32. double Binomial::operator()()
  33. {
  34.     int s = 0;
  35.     for (int i = 0; i < pN; i++) {
  36.     if (pGenerator -> asDouble() < pU) {
  37.         s++;
  38.     }
  39.     }
  40.     return(double(s));
  41. }
  42.  
  43. #include <Geometric.h>
  44.  
  45. double Geometric::operator()()
  46. {
  47.     int samples;
  48.     for (samples = 1; pGenerator -> asDouble() < pMean; samples++);
  49.     return((double) samples);
  50. }
  51.  
  52. #include <HyperGeometric.h>
  53.  
  54. double HyperGeometric::operator()()
  55. {
  56.     double d = (pGenerator -> asDouble() > pP) ? (1.0 - pP) :  (pP);
  57.     return(-pMean * log(pGenerator -> asDouble()) / (2.0 * d) );
  58. }
  59.  
  60. #include <NegativeExpntl.h>
  61.  
  62. double NegativeExpntl::operator()()
  63. {
  64.     return(-pMean * log(pGenerator -> asDouble()));
  65. }
  66.  
  67. #include <Normal.h>
  68.  
  69. //
  70. //    See Simulation, Modelling & Analysis by Law & Kelton, pp259
  71. //
  72. //    This is the ``polar'' method.
  73. // 
  74.  
  75. double Normal::operator()()
  76. {
  77.     
  78.     if (haveCachedNormal == 1) {
  79.     haveCachedNormal = 0;
  80.     return(cachedNormal * pStdDev + pMean );
  81.     } else {
  82.     
  83.     for(;;) {
  84.         double u1 = pGenerator -> asDouble();
  85.         double u2 = pGenerator -> asDouble();
  86.         double v1 = 2 * u1 - 1;
  87.         double v2 = 2 * u2 - 1;
  88.         double w = (v1 * v1) + (v2 * v2);
  89.         
  90. //
  91. //    We actually generate two IID normal distribution variables.
  92. //    We cache the one & return the other.
  93. // 
  94.         if (w <= 1) {
  95.         double y = sqrt( (-2 * log(w)) / w);
  96.         double x1 = v1 * y;
  97.         double x2 = v2 * y;
  98.         
  99.         haveCachedNormal = 1;
  100.         cachedNormal = x2;
  101.         return(x1 * pStdDev + pMean);
  102.         }
  103.     }
  104.     }
  105. }
  106.  
  107. #include <Poisson.h>
  108.  
  109. double Poisson::operator()()
  110. {
  111.     double bound = exp(-1.0 * pMean);
  112.     int count = 0;
  113.  
  114.     for (double product = 1.0;
  115.      product >= bound;
  116.      product *= pGenerator -> asDouble()) {
  117.          count++;
  118.      }
  119.     return(count - 1);
  120. }
  121.  
  122. #include <DiscreteUniform.h>
  123.  
  124. double DiscreteUniform::operator()()
  125. {
  126.     long tmp = long(floor(delta * pGenerator -> asDouble()));
  127.     return( double(pLow + tmp) );
  128. }
  129.  
  130. #include <Uniform.h>
  131.  
  132. double Uniform::operator()()
  133. {
  134.     return( pLow + delta * pGenerator -> asDouble() );
  135. }
  136.  
  137. #include <Weibull.h>
  138.  
  139. //
  140. //    See Simulation, Modelling & Analysis by Law & Kelton, pp259
  141. //
  142. //    This is the ``polar'' method.
  143. // 
  144.  
  145. double Weibull::operator()()
  146. {
  147.     return( pow(pBeta * ( - log(1 - pGenerator -> asDouble()) ), pInvAlpha) );
  148. }
  149.  
  150. #include <Erlang.h>
  151.  
  152. double Erlang::operator()()
  153. {
  154.     double prod = 1.0;
  155.  
  156.     for (int i = 0; i < k; i++) {
  157.     prod *= pGenerator -> asDouble();
  158.     }
  159.     return(-log(prod)/a);
  160. }
  161.  
  162. #include "LogNormal.h"
  163.  
  164. //
  165. //    See Simulation, Modelling & Analysis by Law & Kelton, pp260
  166. //
  167. // 
  168.  
  169. double LogNormal::operator()()
  170. {
  171.     return( pow(M_E, this->Normal::operator()() ) );
  172. }
  173.