DistributionInverter.h

Go to the documentation of this file.

//            QMcBeaver
//
//         Constructed by 
//
//     Michael Todd Feldmann 
//              and 
//   David Randall "Chip" Kent IV
//
// Copyright 2000.  All rights reserved.
//
// drkent@users.sourceforge.net mtfeldmann@users.sourceforge.net

#ifndef DISTRIBUTION_INVERTER_H
#define DISTRIBUTION_INVERTER_H

#define TINY (1e-12)

#include <iostream>
#include <string>
#include <math.h>

#include "Array1D.h"
#include "Random.h"
#include "LinearSpline.h"

using namespace std;

class DistributionInverter
{
public:  
  DistributionInverter(){;} 

  void operator=( const DistributionInverter & rhs );

  void initialize(Array1D<double> &x_input, 
                  Array1D<double> &y_input);

  double random();
  
private:
  //F_inverse=inverse(F); give me a y and I'll give you an x
  LinearSpline F_inverse; 

  void make_F_and_F_inverse(Array1D<double> &x_input, 
                            Array1D<double> &y_input);
};


#endif

---