VectorFunctionSynchronizer<P_V,P_M,Q_V,Q_M>::VectorFunctionSynchronizer(
const BaseVectorFunction<P_V,P_M,Q_V,Q_M>& inputFunction,
const P_V& auxPVec,
const Q_V& auxQVec)
:
m_env (inputFunction.domainSet().env()),
m_vectorFunction(inputFunction),
m_auxPVec (auxPVec),
m_auxQVec (auxQVec)
{
}
StatisticalForwardProblem<P_V,P_M,Q_V,Q_M>::StatisticalForwardProblem(
/*! The prefix */ const char* prefix,
/*! Options (if no input file) */ const SfpOptionsValues* alternativeOptionsValues, // dakota
/*! The input RV */ const BaseVectorRV <P_V,P_M>& paramRv,
/*! The QoI function */ const BaseVectorFunction<P_V,P_M,Q_V,Q_M>& qoiFunction,
/*! The QoI RV */ GenericVectorRV <Q_V,Q_M>& qoiRv)
:
m_env (paramRv.env()),
m_paramRv (paramRv),
m_qoiFunction (qoiFunction),
m_qoiRv (qoiRv),
m_paramChain (NULL),
m_qoiChain (NULL),
m_mcSeqGenerator (NULL),
m_solutionRealizer (NULL),
#ifdef UQ_ALSO_COMPUTE_MDFS_WITHOUT_KDE
m_subMdfGrids (NULL),
m_subMdfValues (NULL),
#endif
#ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
m_subSolutionMdf (NULL),
m_subCdfGrids (NULL),
m_subCdfValues (NULL),
m_subSolutionCdf (NULL),
m_unifiedCdfGrids (NULL),
m_unifiedCdfValues (NULL),
m_unifiedSolutionCdf (NULL),
#endif
m_solutionPdf (NULL),
m_optionsObj (alternativeOptionsValues),
m_userDidNotProvideOptions(false)
{
if (m_env.subDisplayFile()) {
*m_env.subDisplayFile() << "Entering StatisticalForwardProblem<P_V,P_M,Q_V,Q_M>::constructor()"
<< ": prefix = " << prefix
<< ", alternativeOptionsValues = " << alternativeOptionsValues
<< ", m_env.optionsInputFileName() = " << m_env.optionsInputFileName()
<< std::endl;
}
// If NULL, we create one
if (m_optionsObj == NULL) {
SfpOptionsValues * tempOptions = new SfpOptionsValues(&m_env, prefix);
// We did this dance because scanOptionsValues is not a const method, but
// m_optionsObj is a pointer to const
m_optionsObj = tempOptions;
// We do this so we don't delete the user's object in the dtor
m_userDidNotProvideOptions = true;
}
if (m_optionsObj->m_help != "") {
if (m_env.subDisplayFile()) {
*m_env.subDisplayFile() << (*m_optionsObj) << std::endl;
}
}
queso_require_equal_to_msg(paramRv.imageSet().vectorSpace().dimLocal(), qoiFunction.domainSet().vectorSpace().dimLocal(), "'paramRv' and 'qoiFunction' are related to vector spaces of different dimensions");
queso_require_equal_to_msg(qoiFunction.imageSet().vectorSpace().dimLocal(), qoiRv.imageSet().vectorSpace().dimLocal(), "'qoiFunction' and 'qoiRv' are related to vector spaces of different dimensions");
if (m_env.subDisplayFile()) {
*m_env.subDisplayFile() << "Leaving StatisticalForwardProblem<P_V,P_M,Q_V,Q_M>::constructor()"
<< ": prefix = " << m_optionsObj->m_prefix
<< std::endl;
}
}
