queso-0.51.1
List of all members
QUESO::LogNormalVectorRV< V, M > Class Template Reference

A class representing a LogNormal vector RV. More...

#include <LogNormalVectorRV.h>

Inheritance diagram for QUESO::LogNormalVectorRV< V, M >:
Inheritance graph
[legend]
Collaboration diagram for QUESO::LogNormalVectorRV< V, M >:
Collaboration graph
[legend]

Public Member Functions

Constructor/Destructor methods
 LogNormalVectorRV (const char *prefix, const VectorSet< V, M > &imageSet, const V &lawExpVector, const V &lawVarVector)
 Constructor. More...
 
virtual ~LogNormalVectorRV ()
 Virtual destructor. More...
 
I/O methods
void print (std::ostream &os) const
 TODO: Prints the vector RV. More...
 
- Public Member Functions inherited from QUESO::BaseVectorRV< V, M >
 BaseVectorRV (const char *prefix, const VectorSet< V, M > &imageSet)
 Constructor. More...
 
virtual ~BaseVectorRV ()
 Virtual destructor. More...
 
const BaseEnvironment & env () const
 QUESO environment; access to private attribute m_env. More...
 
const VectorSet< V, M > & imageSet () const
 Image set of the vector RV; access to private attribute m_imageSet. More...
 
const BaseJointPdf< V, M > & pdf () const
 Posterior Density Function of the vector RV; access to private attribute m_pdf. More...
 
const BaseVectorRealizer< V, M > & realizer () const
 Finds a realization (sample) of the PDF of this vector RV; access to private attribute m_realizer. More...
 
const BaseVectorCdf< V, M > & subCdf () const
 Finds the Cumulative Distribution Function of this vector RV, considering only the sub-sequence of data; access to private attribute m_subCdf. More...
 
const BaseVectorCdf< V, M > & unifiedCdf () const
 Finds the Cumulative Distribution Function of this vector RV, considering the unified sequence of data; access to private attribute m_unifiedCdf. More...
 
const BaseVectorMdf< V, M > & mdf () const
 Finds the Mass Density Function of this vector RV; access to private attribute m_mdf. More...
 

Additional Inherited Members

- Protected Attributes inherited from QUESO::BaseVectorRV< V, M >
const BaseEnvironment & m_env
 
std::string m_prefix
 
const VectorSet< V, M > & m_imageSet
 
BaseJointPdf< V, M > * m_pdf
 
BaseVectorRealizer< V, M > * m_realizer
 
const BaseVectorCdf< V, M > * m_subCdf
 
const BaseVectorCdf< V, M > * m_unifiedCdf
 
const BaseVectorMdf< V, M > * m_mdf
 

Detailed Description

template<class V, class M>
class QUESO::LogNormalVectorRV< V, M >

A class representing a LogNormal vector RV.

This class allows the user to compute the value of a LogNormal PDF and to generate realizations (samples) from it.
The probability density function of a log-normal distribution is:

\[ f(x|\mu,\sigma) = \frac{1}{x \sigma \sqrt{2 \pi}}\, e^{-\frac{(\ln x - \mu)^2}{2\sigma^2}}, x>0 \]

where the parameters denoted $ \mu $ and $ \sigma $ are, respectively, the mean and standard deviation of the variable's natural logarithm; and x>0.

Definition at line 56 of file LogNormalVectorRV.h.

Constructor & Destructor Documentation

template<class V , class M >
QUESO::LogNormalVectorRV< V, M >::LogNormalVectorRV ( const char *  prefix,
const VectorSet< V, M > &  imageSet,
const V &  lawExpVector,
const V &  lawVarVector 
)

Constructor.

Construct a LogNormal vector RV with mean lawExpVector and diagonal covariance matrix lawVarVector whose variates live in imageSet.

Definition at line 35 of file LogNormalVectorRV.C.

References QUESO::BaseEnvironment::displayVerbosity(), QUESO::BaseVectorRV< V, M >::imageSet(), QUESO::BaseVectorRV< V, M >::m_env, QUESO::BaseVectorRV< V, M >::m_imageSet, QUESO::BaseVectorRV< V, M >::m_mdf, QUESO::BaseVectorRV< V, M >::m_pdf, QUESO::BaseVectorRV< V, M >::m_prefix, QUESO::BaseVectorRV< V, M >::m_realizer, QUESO::BaseVectorRV< V, M >::m_subCdf, QUESO::BaseVectorRV< V, M >::m_unifiedCdf, QUESO::BoxSubset< V, M >::maxValues(), QUESO::BoxSubset< V, M >::minValues(), QUESO::BaseEnvironment::subDisplayFile(), UQ_FATAL_TEST_MACRO, and QUESO::BaseEnvironment::worldRank().

40  :
41  BaseVectorRV<V,M>(((std::string)(prefix)+"gau").c_str(),imageSet)
42 {
43  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
44  *m_env.subDisplayFile() << "Entering LogNormalVectorRV<V,M>::constructor() [1]"
45  << ": prefix = " << m_prefix
46  << std::endl;
47  }
48 
49 // begin kemelli 2013-April-23 --------------------------
50 // LogNormal dist is defined only in (0,inf)
51 
52  const BoxSubset<V,M>* imageBox = dynamic_cast<const BoxSubset<V,M>* >(&imageSet);
53  double smallerOfMaxValues = imageBox->maxValues().getMinValue();
54  double smallerOfMinValues = imageBox->minValues().getMinValue();
55 
56  if( smallerOfMinValues < 0 )
57  {
58  std::cerr << "In LogNormalVectorRV<V,M>::constructor()\n"
59  << "LogNormal distribution is only defined in (0, infinity).\n"
60  << "The data provided is: \n"
61  << *imageBox
62  << "Sampling will not cover all interval.\n"
63  << std::endl;
64 
65 
66  UQ_FATAL_TEST_MACRO(smallerOfMaxValues < 0,
67  m_env.worldRank(),
68  "LogNormalVectorRealizer<V,M>::constructor()",
69  "invalid input: LogNormal distribution is only defined in (0, infinity), and min(m_maxValues)<0");
70 
71  }
72 // end kemelli 2013-April-23 --------------------------
73 
74  m_pdf = new LogNormalJointPdf<V,M>(m_prefix.c_str(),
75  m_imageSet,
76  lawExpVector,
77  lawVarVector);
78 
79  M lowerCholLawCovMatrix(lawVarVector);
80  int iRC = lowerCholLawCovMatrix.chol();
81  lowerCholLawCovMatrix.zeroUpper(false);
82  if (iRC) {
83  std::cerr << "In LogNormalVectorRV<V,M>::constructor() [1]: chol failed, will use svd\n";
84  if (m_env.subDisplayFile()) {
85  *m_env.subDisplayFile() << "In LogNormalVectorRV<V,M>::constructor() [1]: chol failed; will use svd; lawVarVector contents are\n";
86  *m_env.subDisplayFile() << lawVarVector; // FIX ME: might demand parallelism
87  *m_env.subDisplayFile() << std::endl;
88  }
89  M matLaw(lawVarVector);
90  M matU (lawVarVector);
91  M matVt (m_imageSet.vectorSpace().zeroVector());
92  V vecS (m_imageSet.vectorSpace().zeroVector());
93  iRC = matLaw.svd(matU,vecS,matVt);
94  UQ_FATAL_TEST_MACRO(iRC,
95  m_env.worldRank(),
96  "LogNormalVectorRV<V,M>::constructor() [1]",
97  "Cholesky decomposition of covariance matrix failed.");
98 
99  vecS.cwSqrt();
100  m_realizer = new LogNormalVectorRealizer<V,M>(m_prefix.c_str(),
101  m_imageSet,
102  lawExpVector,
103  matU,
104  vecS, // already square rooted
105  matVt);
106  }
107  else {
108  m_realizer = new LogNormalVectorRealizer<V,M>(m_prefix.c_str(),
109  m_imageSet,
110  lawExpVector,
111  lowerCholLawCovMatrix);
112  }
113 
114  m_subCdf = NULL; // FIX ME: complete code
115  m_unifiedCdf = NULL; // FIX ME: complete code
116  m_mdf = NULL; // FIX ME: complete code
117 
118  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
119  *m_env.subDisplayFile() << "Leaving LogNormalVectorRV<V,M>::constructor() [1]"
120  << ": prefix = " << m_prefix
121  << std::endl;
122  }
123 }
QUESO::BaseVectorRV::m_unifiedCdf
const BaseVectorCdf< V, M > * m_unifiedCdf
Definition: VectorRV.h:118
QUESO::BaseEnvironment::worldRank
int worldRank() const
Returns the process world rank.
Definition: Environment.C:235
QUESO::BaseEnvironment::subDisplayFile
std::ofstream * subDisplayFile() const
Access function for m_subDisplayFile (displays file on stream).
Definition: Environment.C:305
QUESO::BaseVectorRV::m_imageSet
const VectorSet< V, M > & m_imageSet
Definition: VectorRV.h:114
QUESO::BaseVectorRV::m_env
const BaseEnvironment & m_env
Definition: VectorRV.h:112
QUESO::BaseVectorRV::m_subCdf
const BaseVectorCdf< V, M > * m_subCdf
Definition: VectorRV.h:117
QUESO::BaseVectorRV::m_realizer
BaseVectorRealizer< V, M > * m_realizer
Definition: VectorRV.h:116
QUESO::BaseVectorRV::imageSet
const VectorSet< V, M > & imageSet() const
Image set of the vector RV; access to private attribute m_imageSet.
Definition: VectorRV.C:79
QUESO::BaseEnvironment::displayVerbosity
unsigned int displayVerbosity() const
Definition: Environment.C:436
QUESO::BaseVectorRV::m_mdf
const BaseVectorMdf< V, M > * m_mdf
Definition: VectorRV.h:119
UQ_FATAL_TEST_MACRO
#define UQ_FATAL_TEST_MACRO(test, givenRank, where, what)
Definition: Defines.h:223
QUESO::BaseVectorRV::m_prefix
std::string m_prefix
Definition: VectorRV.h:113
QUESO::BaseVectorRV::m_pdf
BaseJointPdf< V, M > * m_pdf
Definition: VectorRV.h:115
template<class V , class M >
QUESO::LogNormalVectorRV< V, M >::~LogNormalVectorRV ( )
virtual

Virtual destructor.

Definition at line 126 of file LogNormalVectorRV.C.

127 {
128  delete m_mdf;
129  delete m_unifiedCdf;
130  delete m_subCdf;
131  delete m_realizer;
132  delete m_pdf;
133 }
QUESO::BaseVectorRV::m_unifiedCdf
const BaseVectorCdf< V, M > * m_unifiedCdf
Definition: VectorRV.h:118
QUESO::BaseVectorRV::m_subCdf
const BaseVectorCdf< V, M > * m_subCdf
Definition: VectorRV.h:117
QUESO::BaseVectorRV::m_realizer
BaseVectorRealizer< V, M > * m_realizer
Definition: VectorRV.h:116
QUESO::BaseVectorRV::m_mdf
const BaseVectorMdf< V, M > * m_mdf
Definition: VectorRV.h:119
QUESO::BaseVectorRV::m_pdf
BaseJointPdf< V, M > * m_pdf
Definition: VectorRV.h:115

Member Function Documentation

template<class V , class M >
void QUESO::LogNormalVectorRV< V, M >::print ( std::ostream &  os) const
virtual

TODO: Prints the vector RV.

Todo:
: implement me!

Implements QUESO::BaseVectorRV< V, M >.

Definition at line 137 of file LogNormalVectorRV.C.

138 {
139  os << "LogNormalVectorRV<V,M>::print() says, 'Please implement me.'" << std::endl;
140  return;
141 }
The documentation for this class was generated from the following files:
Generated on Thu Apr 23 2015 19:26:18 for queso-0.51.1 by  doxygen 1.8.5