本文整理汇总了C++中Chain::compute_partialJointPosteriorDensity_overSubSample方法的典型用法代码示例。如果您正苦于以下问题：C++ Chain::compute_partialJointPosteriorDensity_overSubSample方法的具体用法？C++ Chain::compute_partialJointPosteriorDensity_overSubSample怎么用？C++ Chain::compute_partialJointPosteriorDensity_overSubSample使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Chain的用法示例。
在下文中一共展示了Chain::compute_partialJointPosteriorDensity_overSubSample方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: computeJointDensity_MPI_overSubSample
// 'demographicPara' is a 1x6 matrix.
double MaxPosterior::computeJointDensity_MPI_overSubSample(Eigen::MatrixXd demographicPara, IM im, popTree* poptree, Chain coldCh, unsigned int nProcs, unsigned int crr_procID)
{
// REMOVE
//if(crr_procID == 0)
// std::cout << "\n In MaxPosterior::computeJointDensity_MPI_overSubSample\n";
poptree->replacePara(demographicPara);
coldCh.compute_partialJointPosteriorDensity_overSubSample(poptree, im, crr_procID, nProcs);
std::vector<long double> expectationOfEachCoalProb = coldCh.get_expectationOfCoalProb();
unsigned int nloci = coldCh.GetNumLoci();
double posterior = 1;
if(nloci != expectationOfEachCoalProb.size())
{
std::cout << "** Error In MaxPosterior::computeJointDensity_MPI_overSubSample() **\n";
std::cout << "nloci = " << nloci
<< " expectationOfEachCoalProb.size() = " << expectationOfEachCoalProb.size() <<"\n\n";
}
for(unsigned int lc=0; lc< nloci; lc++)
{
double each_local = expectationOfEachCoalProb.at(lc);
double each_global =0;
MPI::COMM_WORLD.Barrier();
MPI::COMM_WORLD.Allreduce(&each_local, &each_global, 1, MPI_DOUBLE, MPI_SUM);
MPI::COMM_WORLD.Barrier();
if(crr_procID ==0)
{
posterior *=each_global;
}
MPI::COMM_WORLD.Barrier();
}
if(crr_procID ==0)
{
Eigen::Vector3d paraMax = im.get_paraMax();
double priorPopTree = poptree->computeJointPrior(paraMax);
posterior *=priorPopTree;
}
if(crr_procID ==0)
{
totalComputingTime_eigen += coldCh.get_eachComputingTime_eigen();
totalComputingTime_condiProb += coldCh.get_eachComputingTime_condiProb();
totalNum_eigenFunctionCalls++;
totalNum_condiProbFunctionCalls += coldCh.get_countCondiProbFunctionCalls();
// std::cout << "totalComputingTime_eigen = " << totalComputingTime_eigen.count()/1000 << "milliseconds\n";
}
return posterior;
}