ocsRanNumGen.hh

Go to the documentation of this file.

// $Id: ocsRanNumGen.hh,v 1.2 2008/08/08 02:23:11 zweck Exp $

//###################################################################
//
//           Optical Communication Systems Simulator
//
//       Copyright (2000):
//       Optical Fiber Communications Laboratory (OFCL)
//       Computer Science & Electrical Engineering Department (CSEE)
//       University of Maryland Baltimore County (UMBC)
//
//###################################################################



//## Method for generation of random number for importance 
//#  sampling by Ivan Lima

#ifndef _OCS_RAN_NUM_GEN_HH_
#define  _OCS_RAN_NUM_GEN_HH_

#ifndef _OCS_SYSTEM_HDRS_
#define _OCS_SYSTEM_HDRS_

#include <iostream.h>
#include <string.h>
#include <math.h> 
#include <fstream.h>
#include <complex.h>
#include <stdlib.h>
#include <stdio.h>
#include <strstream>

#include <string>
  using namespace std;


#endif /* _OCS_SYSTEM_HDRS_ */

#include "ocsConst.hh"
#include "ocsReducedModelTools.hh"
#include "ocsReadWriteTools.hh"



class RanNumGen
{

public:

  
  RanNumGen(string InFileName);


  ~RanNumGen();


  void SaveCurrentState(void);

  

  void ResetRanNumGen(void);

  float GetRanNum(void);




  float GetRanCosineBiased(double alpha, double beta, 
                           double *LikelihoodRatio);
  float GetLikelihoodRatioForRanCosineBiased(double alpha,
                               double ValFunction);                        
                           
  float GetRanSignBiased(double alpha, 
                         double *LikelihoodRatio);


  float GetRanCosineBiased(double alpha, double *LikelihoodRatio);


  float GetGaussianDeviate(void);


  void ReadRanSeqParameters(string InFileName);


  void WriteRanSeqParameters(string OutFileName);


private:

  static const long IM1 = 2147483563;
  static const long IM2 = 2147483399;
  static const double AM = 1.0/2147483563;
  static const long IMM1 = IM1 - 1;
  static const long IA1 =  40014;
  static const long IA2 = 40692;
  static const long IQ1 = 53668;
  static const long IQ2 = 52774;
  static const long IR1 = 12211;
  static const long IR2 = 3791;
  static const int NTAB = 32;
  static const long NDIV = (1+IMM1/NTAB);
  //static const double EPS = 1.2e-7;
  //static const double RNMX = 1.0-EPS;

  static const double RNMX = 1.0 - 1.2e-7;

  long idum,  idumStart;
  long idum2, idum2Start;
  long iy, iyStart;
  long iv[NTAB], ivStart[NTAB];


  string ContinuationFileName;


  string RestartFileName;


  int DebugLevel;

};

#endif  /* _OCS_RAN_NUM_GEN_HH_ */