home *** CD-ROM | disk | FTP | other *** search
- // This may look like C code, but it is really -*- C++ -*-
- /*
- Copyright (C) 1988 Free Software Foundation
- written by Dirk Grunwald (grunwald@cs.uiuc.edu)
-
- This file is part of GNU CC.
-
- GNU CC is distributed in the hope that it will be useful,
- but WITHOUT ANY WARRANTY. No author or distributor
- accepts responsibility to anyone for the consequences of using it
- or for whether it serves any particular purpose or works at all,
- unless he says so in writing. Refer to the GNU CC General Public
- License for full details.
-
- Everyone is granted permission to copy, modify and redistribute
- GNU CC, but only under the conditions described in the
- GNU CC General Public License. A copy of this license is
- supposed to have been given to you along with GNU CC so you
- can know your rights and responsibilities. It should be in a
- file named COPYING. Among other things, the copyright notice
- and this notice must be preserved on all copies.
- */
- #ifndef _Erlang_h
- #ifdef __GNUG__
- #pragma once
- #pragma interface
- #endif
- #define _Erlang_h 1
-
- #include <_Random.h>
-
- class Erlang: public Random {
- protected:
- double pMean;
- double pVariance;
- int k;
- double a;
- void setState();
- public:
- Erlang(double mean, double variance, RNG *gen);
-
- double mean();
- double mean(double x);
- double variance();
- double variance(double x);
-
- virtual double operator()();
-
- };
-
- #if defined(__OPTIMIZE__) || defined(USE_LIBGXX_INLINES)
-
- inline void Erlang::setState() {
- k = int( (pMean * pMean ) / pVariance + 0.5 );
- k = (k > 0) ? k : 1;
- a = k / pMean;
- }
-
- inline Erlang::Erlang(double mean, double variance, RNG *gen) : (gen)
- {
- pMean = mean; pVariance = variance;
- setState();
- }
-
- inline double Erlang::mean() { return pMean; }
- inline double Erlang::mean(double x) {
- double tmp = pMean; pMean = x; setState(); return tmp;
- };
-
- inline double Erlang::variance() { return pVariance; }
- inline double Erlang::variance(double x) {
- double tmp = pVariance; pVariance = x; setState(); return tmp;
- }
-
-
- #endif
- #endif
-