Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
TUnuranMultiContDist.cxx
Go to the documentation of this file.
1// @(#)root/unuran:$Id$
2// Authors: L. Moneta, J. Leydold Wed Feb 28 2007
3
4/**********************************************************************
5 * *
6 * Copyright (c) 2006 LCG ROOT Math Team, CERN/PH-SFT *
7 * *
8 * *
9 **********************************************************************/
10
11// Implementation file for class TUnuranMultiContDist
12
15
16#include "TF1.h"
17#include <cassert>
18
19
21 fPdf(&pdf),
22 fIsLogPdf(isLogPdf),
23 fOwnFunc(false)
24{
25 //Constructor from generic function interfaces
26}
27
28
29TUnuranMultiContDist::TUnuranMultiContDist (TF1 * func, unsigned int dim, bool isLogPdf) :
30 fPdf(0),
31 fIsLogPdf(isLogPdf),
32 fOwnFunc(false)
33{
34 //Constructor from a TF1 objects
35 if (func) {
36 fPdf = new ROOT::Math::WrappedMultiTF1( *func, dim);
37 fOwnFunc = true;
38 }
39}
40
41
42
45 fPdf(0)
46{
47 // Implementation of copy ctor using assignment operator
48 operator=(rhs);
49}
50
52{
53 // Implementation of assignment operator (copy only the function pointer not the function itself)
54 if (this == &rhs) return *this; // time saving self-test
55 fXmin = rhs.fXmin;
56 fXmax = rhs.fXmax;
57 fMode = rhs.fMode;
58 fIsLogPdf = rhs.fIsLogPdf;
59 fOwnFunc = rhs.fOwnFunc;
60 if (!fOwnFunc)
61 fPdf = rhs.fPdf;
62 else {
63 if (fPdf) delete fPdf;
64 fPdf = (rhs.fPdf) ? rhs.fPdf->Clone() : 0;
65 }
66 return *this;
67}
68
70 // destructor implementation
71 if (fOwnFunc && fPdf) delete fPdf;
72}
73
74
75double TUnuranMultiContDist::Pdf ( const double * x) const {
76 // evaluate the distribution
77 assert(fPdf != 0);
78 return (*fPdf)(x);
79}
80
81
82void TUnuranMultiContDist::Gradient( const double * x, double * grad) const {
83 // do numerical derivation and return gradient in vector grad
84 // grad must point to a vector of at least ndim size
85 unsigned int ndim = NDim();
86 for (unsigned int i = 0; i < ndim; ++i)
87 grad[i] = Derivative(x,i);
88
89 return;
90}
91
92double TUnuranMultiContDist::Derivative( const double * x, int coord) const {
93 // do numerical derivation of gradient using 5 point rule
94 // use 5 point rule
95
96 //double eps = 0.001;
97 //const double kC1 = 8*std::numeric_limits<double>::epsilon();
98 assert(fPdf != 0);
99
100 double h = 0.001;
101
102 std::vector<double> xx(NDim() );
103
104 xx[coord] = x[coord]+h; double f1 = (*fPdf)(&xx.front());
105 xx[coord] = x[coord]-h; double f2 = (*fPdf)(&xx.front());
106
107 xx[coord] = x[coord]+h/2; double g1 = (*fPdf)(&xx.front());
108 xx[coord] = x[coord]-h/2; double g2 = (*fPdf)(&xx.front());
109
110 //compute the central differences
111 double h2 = 1/(2.*h);
112 double d0 = f1 - f2;
113 double d2 = 2*(g1 - g2);
114 //double error = kC1*h2*fx; //compute the error
115 double deriv = h2*(4*d2 - d0)/3.;
116 return deriv;
117}
118
119
120
#define h(i)
Definition RSha256.hxx:106
Documentation for the abstract class IBaseFunctionMultiDim.
Definition IFunction.h:62
virtual IBaseFunctionMultiDimTempl< T > * Clone() const =0
Clone a function.
1-Dim function class
Definition TF1.h:213
TUnuranBaseDist, base class for Unuran distribution classees such as TUnuranContDist (for one-dimensi...
TUnuranMultiContDist class describing multi dimensional continuous distributions.
double Pdf(const double *x) const
evaluate the probability density function, used by UnuRan
double Derivative(const double *x, int icoord) const
evaluate the partial derivative for the given coordinate.
std::vector< double > fXmax
void Gradient(const double *x, double *grad) const
evaluate the gradient vector of the Pdf.
std::vector< double > fXmin
std::vector< double > fMode
TUnuranMultiContDist & operator=(const TUnuranMultiContDist &rhs)
Assignment operator.
const ROOT::Math::IMultiGenFunction * fPdf
unsigned int NDim() const
get number of dimension of the distribution
virtual ~TUnuranMultiContDist()
Destructor.
TUnuranMultiContDist(TF1 *func=0, unsigned int dim=0, bool isLogPdf=false)
Constructor from a TF1 object representing the Probability density function.
Double_t x[n]
Definition legend1.C:17
TF1 * f1
Definition legend1.C:11
WrappedMultiTF1Templ< double > WrappedMultiTF1