#include <GcmTotalInfo.h>

Collaboration diagram for QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >:

Public Member Functions
	GcmTotalInfo (const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &s)

	GcmTotalInfo (const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &s, const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &e)

	~GcmTotalInfo ()

unsigned int	initializeTotalDim (const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &s)

unsigned int	initializeTotalDim (const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &s, const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &e)

Public Attributes
const BaseEnvironment &	m_env

unsigned int	m_numConstituents

std::vector< const VectorSet < P_V, P_M > * >	m_constitutiveDomains

std::vector< const BaseVectorRV< P_V, P_M > * >	m_constitutivePriorRvs

unsigned int	m_totalDim

VectorSpace< P_V, P_M >	m_totalSpace

double	m_totalDomainVolume

ConcatenationSubset< P_V, P_M >	m_totalDomain

ConcatenatedVectorRV< P_V, P_M >	m_totalPriorRv

GenericVectorRV< P_V, P_M >	m_totalPostRv

P_V	m_like_previousTotal

P_V	m_totalPostMean

P_V	m_totalPostMedian

P_V	m_totalPostMode

double	m_totalPostMaxLnValue

P_V	m_totalMLE

double	m_totalLikeMaxLnValue

VectorSet< P_V, P_M > *	m_solutionDomain

BaseJointPdf< P_V, P_M > *	m_solutionPdf

BaseVectorRealizer< P_V, P_M > *	m_solutionRealizer

MetropolisHastingsSG< P_V, P_M > *	m_mhSeqGenerator

MLSampling< P_V, P_M > *	m_mlSampler

BaseVectorSequence< P_V, P_M > *	m_chain

Private Member Functions
void	commonConstructor ()

Detailed Description

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
class QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >

Definition at line 44 of file GcmTotalInfo.h.

Constructor & Destructor Documentation

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmTotalInfo ( const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > & s )

Definition at line 32 of file GcmTotalInfo.C.

References QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::commonConstructor(), QUESO::BaseEnvironment::displayVerbosity(), QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_env, and QUESO::BaseEnvironment::subDisplayFile().

   :
   m_env                 (s.m_env),
   m_numConstituents     (4),
   m_constitutiveDomains (m_numConstituents,(const VectorSet   <P_V,P_M>*) NULL),
   m_constitutivePriorRvs(m_numConstituents,(const BaseVectorRV<P_V,P_M>*) NULL),
   m_totalDim            (initializeTotalDim(s)),
   m_totalSpace          (m_env,"total_",m_totalDim,NULL),
   m_totalDomainVolume   (1.),
   m_totalDomain         ("total_",m_totalSpace,m_totalDomainVolume,m_constitutiveDomains),
   m_totalPriorRv        ("total_prior_",m_constitutivePriorRvs,m_totalDomain),
   m_totalPostRv         ("total_post_",m_totalSpace),
   m_like_previousTotal  (m_totalSpace.zeroVector()),
   m_totalPostMean       (m_totalSpace.zeroVector()),
   m_totalPostMedian     (m_totalSpace.zeroVector()),
   m_totalPostMode       (m_totalSpace.zeroVector()),
   m_totalPostMaxLnValue (-INFINITY),
   m_totalMLE            (m_totalSpace.zeroVector()),
   m_totalLikeMaxLnValue (-INFINITY),
   m_solutionDomain      (NULL),
   m_solutionPdf         (NULL),
   m_solutionRealizer    (NULL),
   m_mhSeqGenerator      (NULL),
   m_mlSampler           (NULL),
   m_chain               (NULL)
 {
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Entering GcmTotalInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor(1)"
                             << std::endl;
   }
 
   commonConstructor();
 
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Leaving GcmTotalInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor(1)"
                             << std::endl;
   }
 }

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmTotalInfo	(	const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &	s,
		const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &	e
	)

Definition at line 73 of file GcmTotalInfo.C.

References QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::commonConstructor(), QUESO::BaseEnvironment::displayVerbosity(), QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_env, and QUESO::BaseEnvironment::subDisplayFile().

   :
   m_env                 (s.m_env),
   m_numConstituents     (8),
   m_constitutiveDomains (m_numConstituents,(const VectorSet<P_V,P_M>*) NULL),
   m_constitutivePriorRvs(m_numConstituents,(const BaseVectorRV<P_V,P_M>*) NULL),
   m_totalDim            (initializeTotalDim(s,e)),
   m_totalSpace          (m_env, "total_", m_totalDim, NULL),
   m_totalDomainVolume   (1.),
   m_totalDomain         ("total_",m_totalSpace,m_totalDomainVolume,m_constitutiveDomains),
   m_totalPriorRv        ("total_prior_",m_constitutivePriorRvs,m_totalDomain),
   m_totalPostRv         ("total_post_",m_totalSpace),
   m_like_previousTotal  (m_totalSpace.zeroVector()),
   m_totalPostMean       (m_totalSpace.zeroVector()),
   m_totalPostMedian     (m_totalSpace.zeroVector()),
   m_totalPostMode       (m_totalSpace.zeroVector()),
   m_totalPostMaxLnValue (-INFINITY),
   m_totalMLE            (m_totalSpace.zeroVector()),
   m_totalLikeMaxLnValue (-INFINITY),
   m_solutionDomain      (NULL),
   m_solutionPdf         (NULL),
   m_solutionRealizer    (NULL),
   m_mhSeqGenerator      (NULL),
   m_mlSampler           (NULL),
   m_chain               (NULL)
 {
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Entering GcmTotalInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor(2)"
                             << std::endl;
   }
 
   commonConstructor();
 
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "Leaving GcmTotalInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor(2)"
                             << std::endl;
   }
 }

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::~GcmTotalInfo ( )

Definition at line 115 of file GcmTotalInfo.C.

 {
   if (m_chain) {
     m_chain->clear();
     delete m_chain;
   }
   if (m_mlSampler       ) delete m_mlSampler;
   if (m_mhSeqGenerator  ) delete m_mhSeqGenerator;
   if (m_solutionRealizer) delete m_solutionRealizer;
   if (m_solutionPdf     ) delete m_solutionPdf;
   if (m_solutionDomain  ) delete m_solutionDomain;
 }

Member Function Documentation

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

void QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::commonConstructor ( )

private

Definition at line 175 of file GcmTotalInfo.C.

Referenced by QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmTotalInfo().

 {
   m_totalPriorRv.pdf().setNormalizationStyle(0); // CSRI - 2013-Aug-06 - with Laura
   for (unsigned int i = 0; i < m_numConstituents; ++i) {
     m_totalDomainVolume *= m_constitutiveDomains[i]->volume();
   }
 
   //********************************************************************************
   // Display information
   //********************************************************************************
   if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
     *m_env.subDisplayFile() << "In GcmTotalInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::commonConstructor()"
                             << "\n  m_totalDim          = "  << m_totalDim
                             << "\n  m_numConstituents   = "  << m_numConstituents
                             << "\n  m_totalDomainVolume = "  << m_totalDomainVolume
                             << std::endl;
   }
 
   return;
 }

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

unsigned int QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::initializeTotalDim ( const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > & s )

Definition at line 130 of file GcmTotalInfo.C.

 {
   m_constitutiveDomains[0] = &(s.m_1lambdaEtaDomain);
   m_constitutiveDomains[1] = &(s.m_2lambdaWDomain);
   m_constitutiveDomains[2] = &(s.m_3rhoWDomain);
   m_constitutiveDomains[3] = &(s.m_4lambdaSDomain);
 
   m_constitutivePriorRvs[0] = &(s.m_1lambdaEtaPriorRv);
   m_constitutivePriorRvs[1] = &(s.m_2lambdaWPriorRv);
   m_constitutivePriorRvs[2] = &(s.m_3rhoWPriorRv);
   m_constitutivePriorRvs[3] = &(s.m_4lambdaSPriorRv);
 
   return (s.m_1lambdaEtaDim + s.m_2lambdaWDim + s.m_3rhoWDim + s.m_4lambdaSDim);
 }

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >

unsigned int QUESO::GcmTotalInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::initializeTotalDim	(	const GcmSimulationInfo< S_V, S_M, P_V, P_M, Q_V, Q_M > &	s,
		const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &	e
	)

Definition at line 148 of file GcmTotalInfo.C.

 {
   m_constitutiveDomains[0] = &(s.m_1lambdaEtaDomain);
   m_constitutiveDomains[1] = &(s.m_2lambdaWDomain);
   m_constitutiveDomains[2] = &(s.m_3rhoWDomain);
   m_constitutiveDomains[3] = &(s.m_4lambdaSDomain);
   m_constitutiveDomains[4] = &(e.m_5lambdaYDomain);
   m_constitutiveDomains[5] = &(e.m_6lambdaVDomain);
   m_constitutiveDomains[6] = &(e.m_7rhoVDomain);
   m_constitutiveDomains[7] = &(e.m_8thetaPriorRv.imageSet());
 
   m_constitutivePriorRvs[0] = &(s.m_1lambdaEtaPriorRv);
   m_constitutivePriorRvs[1] = &(s.m_2lambdaWPriorRv);
   m_constitutivePriorRvs[2] = &(s.m_3rhoWPriorRv);
   m_constitutivePriorRvs[3] = &(s.m_4lambdaSPriorRv);
   m_constitutivePriorRvs[4] = &(e.m_5lambdaYPriorRv);
   m_constitutivePriorRvs[5] = &(e.m_6lambdaVPriorRv);
   m_constitutivePriorRvs[6] = &(e.m_7rhoVPriorRv);
   m_constitutivePriorRvs[7] = &(e.m_8thetaPriorRv);
 
   return (s.m_1lambdaEtaDim + s.m_2lambdaWDim + s.m_3rhoWDim + s.m_4lambdaSDim + e.m_5lambdaYDim + e.m_6lambdaVDim + e.m_7rhoVDim + e.m_8thetaDim);
 }