queso-0.57.1
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Pages
BasicPdfsGsl.C
Go to the documentation of this file.
1 //-----------------------------------------------------------------------bl-
2 //--------------------------------------------------------------------------
3 //
4 // QUESO - a library to support the Quantification of Uncertainty
5 // for Estimation, Simulation and Optimization
6 //
7 // Copyright (C) 2008-2017 The PECOS Development Team
8 //
9 // This library is free software; you can redistribute it and/or
10 // modify it under the terms of the Version 2.1 GNU Lesser General
11 // Public License as published by the Free Software Foundation.
12 //
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Lesser General Public License for more details.
17 //
18 // You should have received a copy of the GNU Lesser General Public
19 // License along with this library; if not, write to the Free Software
20 // Foundation, Inc. 51 Franklin Street, Fifth Floor,
21 // Boston, MA 02110-1301 USA
22 //
23 //-----------------------------------------------------------------------el-
24 
25 #include <queso/BasicPdfsGsl.h>
26 #include <queso/Defines.h>
27 #include <gsl/gsl_randist.h>
28 #include <math.h>
29 
30 namespace QUESO {
31 
34  :
35  BasicPdfsBase(worldRank)
36 {
37 }
38 
39 // Destructor ---------------------------------------
41 {
42 }
43 
44 // --------------------------------------------------
45 double
46 BasicPdfsGsl::betaPdfActualValue(double x, double alpha, double beta) const
47 {
48  double result = gsl_ran_beta_pdf(x,alpha,beta);
49  if (isinf(result)) { // CSRI - 2013-aug-06, with Laura
50  std::cerr << "In BasicPdfsGsl::betaPdfActualValue(): hitting inf"
51  << ", x = " << x
52  << ", alpha = " << alpha
53  << ", beta = " << beta
54  << std::endl;
55  //result = gsl_ran_beta_pdf(x,alpha,beta);
56  result = 0.;
57  }
58  return result;
59 }
60 
61 // --------------------------------------------------
62 double
63 BasicPdfsGsl::gammaPdfActualValue(double x, double a, double b) const
64 {
65  return gsl_ran_gamma_pdf(x,a,b);
66 }
67 
68 } // End namespace QUESO
double gammaPdfActualValue(double x, double a, double b) const
TODO: Actual value of the Gamma PDF.
Definition: BasicPdfsGsl.C:63
~BasicPdfsGsl()
Destructor.
Definition: BasicPdfsGsl.C:40
TODO: Base class for basic PDFs (via either GSL or Boost).
Definition: BasicPdfsBase.h:44
double betaPdfActualValue(double x, double alpha, double beta) const
TODO: Actual value of the Beta PDF.
Definition: BasicPdfsGsl.C:46
BasicPdfsGsl()
Default constructor.

Generated on Tue Jun 5 2018 19:48:54 for queso-0.57.1 by  doxygen 1.8.5