RanlibCPP.h

Go to the documentation of this file.

// Contents: Wrapper for ranlib functions
//           
//
// 
//
// Principal Author: Alexei Miasnikov (2002)
//
// Status: in progress
//

//
//
//   $Id: RanlibCPP.h,v 1.5 2007/09/27 16:44:03 amiasnik Exp $
//
//

#ifndef _RANLIB_H_
#define _RANLIB_H_


#include "ranlib.h"
#include "cdflib.h"
#include <time.h>
#include <unistd.h>
#include <stdlib.h>
#include <iostream>


class RandLibURG
{
 public:

  RandLibURG(long i1,long i2) {
    setall(i1,i2);
  }

  RandLibURG() {
    setall(time(NULL),time(NULL));
    long low = irand(0,2147483560);
    long high = irand(0,2147483560);
    setall(low,high);
  }

  void reset(long i1,long i2) {
    setall(i1,i2);
  }

  void reset(){
    setall(time(NULL),time(NULL));
    long low = irand(0,2147483560);
    long high = irand(0,2147483560);
    setall(low,high); 
  }

  void setSeedPID(){
    long pid  = getpid();
    long low  = labs( (pid + (pid % 131) * time(NULL)) % 2147483560); // pid +  ... is to make shure that 
    long high = labs( (pid + (pid % 17)  * time(NULL)) % 2147483560); // that seed will not reduce to 0
    setall(low,high); 
    low = irand(0,2147483560); 
    high = irand(0,2147483560); 
    setall(low,high);  
  } 

  void getseed( long& s1, long& s2 ){
    long iseed1,iseed2;
    getsd(&iseed1,&iseed2);
    s1 = iseed1; s2 = iseed2;
  }

  float rand() {
    return ranf(); 
  }

  
  float rand(float low,float high) {
    return genunf( low, high);
  }

  
  long irand(long low,long high){
   return  ignuin( low, high);
  }
};


class RandLib
{
 public:

  static RandLibURG ur;

  
  static double chiPValue( int idf, double alpha ){
    int which = 2;
    double p = 1.0 - alpha;
    double q = alpha;
    double x;
    double df = idf;
    int status;
    double bound;
   
    cdfchi(&which,&p,&q,&x,&df,&status,&bound);

    return x;
  }

};
#endif