A class for handling Inverse Gamma joint PDFs. More...

#include <InverseGammaJointPdf.h>

Inheritance diagram for QUESO::InverseGammaJointPdf< V, M >:

Collaboration diagram for QUESO::InverseGammaJointPdf< V, M >:

Public Member Functions
Constructor/Destructor methods
	InverseGammaJointPdf (const char *prefix, const VectorSet< V, M > &domainSet, const V &alpha, const V &beta)
	Constructor. More...

	~InverseGammaJointPdf ()
	Destructor. More...

Math methods
double	actualValue (const V &domainVector, const V domainDirection, V gradVector, M hessianMatrix, V hessianEffect) const
	Actual value of the Gamma PDF. More...

double	lnValue (const V &domainVector, const V domainDirection, V gradVector, M hessianMatrix, V hessianEffect) const
	TODO: Logarithm of the value of the Gamma PDF. More...

double	computeLogOfNormalizationFactor (unsigned int numSamples, bool updateFactorInternally) const
	Computes the logarithm of the normalization factor. More...

Public Member Functions inherited from QUESO::BaseJointPdf< V, M >
	BaseJointPdf (const char *prefix, const VectorSet< V, M > &domainSet)
	Default constructor. More...

virtual	~BaseJointPdf ()
	Destructor. More...

virtual void	setNormalizationStyle (unsigned int value) const
	Sets a value to be used in the normalization style (stored in the protected attribute m_normalizationStyle.) More...

void	setLogOfNormalizationFactor (double value) const
	Sets a logarithmic value to be used in the normalization factor (stored in the protected attribute m_normalizationStyle.) More...

Public Member Functions inherited from QUESO::BaseScalarFunction< V, M >
	BaseScalarFunction (const char *prefix, const VectorSet< V, M > &domainSet)
	Default constructor. More...

virtual	~BaseScalarFunction ()
	Destructor. More...

const VectorSet< V, M > &	domainSet () const
	Access to the protected attribute `m_domainSet:` domain set of the scalar function. More...

Protected Attributes
V	m_alpha

V	m_beta

Protected Attributes inherited from QUESO::BaseJointPdf< V, M >
unsigned int	m_normalizationStyle

double	m_logOfNormalizationFactor

Protected Attributes inherited from QUESO::BaseScalarFunction< V, M >
const BaseEnvironment &	m_env

std::string	m_prefix

const VectorSet< V, M > &	m_domainSet
	Domain set of the scalar function. More...

Additional Inherited Members
Protected Member Functions inherited from QUESO::BaseJointPdf< V, M >
double	commonComputeLogOfNormalizationFactor (unsigned int numSamples, bool updateFactorInternally) const
	Common method (to the derived classes) to compute the logarithm of the normalization factor. More...

Detailed Description

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

A class for handling Inverse Gamma joint PDFs.

This class allows the mathematical definition of an Inverse Gamma Joint PDF.

Definition at line 49 of file InverseGammaJointPdf.h.

Constructor & Destructor Documentation

template<class V , class M >

QUESO::InverseGammaJointPdf< V, M >::InverseGammaJointPdf	(	const char *	prefix,
		const VectorSet< V, M > &	domainSet,
		const V &	alpha,
		const V &	beta
	)

Constructor.

Constructs a new object of the class, given a prefix, the domain set, and the parameters alpha and beta of the Inverse Gamma PDF.

Definition at line 33 of file InverseGammaJointPdf.C.

References QUESO::BaseEnvironment::displayVerbosity(), QUESO::BaseScalarFunction< V, M >::m_env, QUESO::BaseScalarFunction< V, M >::m_prefix, and QUESO::BaseEnvironment::subDisplayFile().

   :
   BaseJointPdf<V,M>(((std::string)(prefix)+"uni").c_str(),domainSet),
   m_alpha(alpha),
   m_beta (beta)
 {
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
     *m_env.subDisplayFile() << "Entering InverseGammaJointPdf<V,M>::constructor()"
                             << ": prefix = " << m_prefix
                             << std::endl;
   }
 
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
     *m_env.subDisplayFile() << "Leaving InverseGammaJointPdf<V,M>::constructor()"
                             << ": prefix = " << m_prefix
                             << std::endl;
   }
 }

template<class V , class M >

QUESO::InverseGammaJointPdf< V, M >::~InverseGammaJointPdf ( )

Destructor.

Definition at line 57 of file InverseGammaJointPdf.C.

58 {

59 }

Member Function Documentation

template<class V , class M >

double QUESO::InverseGammaJointPdf< V, M >::actualValue	(	const V &	domainVector,
		const V *	domainDirection,
		V *	gradVector,
		M *	hessianMatrix,
		V *	hessianEffect
	)		const

virtual

Actual value of the Gamma PDF.

This routine calls method lnValue() and returns the exponent of the returning value of such method.

Implements QUESO::BaseJointPdf< V, M >.

Definition at line 63 of file InverseGammaJointPdf.C.

References UQ_FATAL_TEST_MACRO.

 {
   UQ_FATAL_TEST_MACRO(domainVector.sizeLocal() != this->m_domainSet.vectorSpace().dimLocal(),
                       m_env.worldRank(),
                       "InverseGammaJointPdf<V,M>::actualValue()",
                       "invalid input");
 
   UQ_FATAL_TEST_MACRO((domainDirection || gradVector || hessianMatrix || hessianEffect),
                       m_env.worldRank(),
                       "InverseGammaJointPdf<V,M>::actualValue()",
                       "incomplete code for gradVector, hessianMatrix and hessianEffect calculations");
 
   // No need to multiply by exp(m_logOfNormalizationFactor) because 'lnValue()' is called [PDF-07]
   return exp(this->lnValue(domainVector,domainDirection,gradVector,hessianMatrix,hessianEffect));
 }

template<class V , class M >

double QUESO::InverseGammaJointPdf< V, M >::computeLogOfNormalizationFactor	(	unsigned int	numSamples,
		bool	updateFactorInternally
	)		const

virtual

Computes the logarithm of the normalization factor.

This routine calls BaseJointPdf::commonComputeLogOfNormalizationFactor().

Implements QUESO::BaseJointPdf< V, M >.

Definition at line 113 of file InverseGammaJointPdf.C.

References QUESO::BaseJointPdf< V, M >::commonComputeLogOfNormalizationFactor().

 {
   double value = 0.;
 
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Entering InverseGammaJointPdf<V,M>::computeLogOfNormalizationFactor()"
                             << std::endl;
   }
   value = BaseJointPdf<V,M>::commonComputeLogOfNormalizationFactor(numSamples, updateFactorInternally);
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Leaving InverseGammaJointPdf<V,M>::computeLogOfNormalizationFactor()"
                             << ", m_logOfNormalizationFactor = " << m_logOfNormalizationFactor
                             << std::endl;
   }
 
   return value;
 }

template<class V , class M >

double QUESO::InverseGammaJointPdf< V, M >::lnValue	(	const V &	domainVector,
		const V *	domainDirection,
		V *	gradVector,
		M *	hessianMatrix,
		V *	hessianEffect
	)		const

virtual

TODO: Logarithm of the value of the Gamma PDF.

Todo:: : implement me!

Implements QUESO::BaseJointPdf< V, M >.

Definition at line 86 of file InverseGammaJointPdf.C.

References UQ_FATAL_TEST_MACRO.

 {
   UQ_FATAL_TEST_MACRO((domainDirection || gradVector || hessianMatrix || hessianEffect),
                       m_env.worldRank(),
                       "InverseGammaJointPdf<V,M>::lnValue()",
                       "incomplete code for gradVector, hessianMatrix and hessianEffect calculations");
 
   double result = 0.;
   for (unsigned int i = 0; i < domainVector.sizeLocal(); ++i) {
     result -= (m_alpha[i]+1.)*log(domainVector[i]);
     result -= m_beta[i]/domainVector[i];
     if (m_normalizationStyle == 0) {
       // Code needs to be done yet
     }
   }
   result += m_logOfNormalizationFactor; // [PDF-07]
 
   return result;
 }