home *** CD-ROM | disk | FTP | other *** search
/ GEMini Atari / GEMini_Atari_CD-ROM_Walnut_Creek_December_1993.iso / files / gnu / g__inc / xhypgeom.h < prev    next >
Encoding:
C/C++ Source or Header  |  1993-07-23  |  2.0 KB  |  74 lines
  1. // This may look like C code, but it is really -*- C++ -*-
  2. /* 
  3. Copyright (C) 1988 Free Software Foundation
  4.     written by Dirk Grunwald (grunwald@cs.uiuc.edu)
  5.  
  6. This file is part of the GNU C++ Library.  This library is free
  7. software; you can redistribute it and/or modify it under the terms of
  8. the GNU Library General Public License as published by the Free
  9. Software Foundation; either version 2 of the License, or (at your
  10. option) any later version.  This library is distributed in the hope
  11. that it will be useful, but WITHOUT ANY WARRANTY; without even the
  12. implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  13. PURPOSE.  See the GNU Library General Public License for more details.
  14. You should have received a copy of the GNU Library General Public
  15. License along with this library; if not, write to the Free Software
  16. Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  17. */
  18. #ifndef _HyperGeometric_h
  19. #ifdef __GNUG__
  20. #pragma once
  21. #pragma interface
  22. #endif
  23. #define _HyperGeometric_h 
  24.  
  25. #include <xrandom.h>
  26.  
  27. class HyperGeometric: public Random {
  28. protected:
  29.     double pMean;
  30.     double pVariance;
  31.     double pP;
  32.     void setState();
  33.  
  34. public:
  35.     HyperGeometric(double mean, double variance, RNG *gen);
  36.  
  37.     double mean();
  38.     double mean(double x);
  39.     double variance();
  40.     double variance(double x);
  41.  
  42.     virtual double operator()();
  43. };
  44.  
  45. #if defined(__OPTIMIZE__) || defined(USE_LIBGXX_INLINES)
  46.  
  47. inline void HyperGeometric::setState() {
  48.   double z = pVariance / (pMean * pMean);
  49.   pP = 0.5 * (1.0 - sqrt((z - 1.0) / ( z + 1.0 )));
  50. }
  51.  
  52. inline HyperGeometric::HyperGeometric(double mean, double variance, RNG *gen)
  53. : (gen) {
  54.   pMean = mean; pVariance = variance;
  55.   setState();
  56. }
  57.  
  58. inline double HyperGeometric::mean() { return pMean; };
  59.  
  60. inline double HyperGeometric::mean(double x) {
  61.   double t = pMean; pMean = x;
  62.   setState(); return t;
  63. }
  64.  
  65. inline double HyperGeometric::variance() { return pVariance; }
  66.  
  67. inline double HyperGeometric::variance(double x) {
  68.   double t = pVariance; pVariance = x;
  69.   setState(); return t;
  70. }
  71.  
  72. #endif
  73. #endif
  74.