AlignedPhyloLikelihoodHmm.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: The Bio++ Development Group
//
// SPDX-License-Identifier: CECILL-2.1
 
#include <Bpp/Numeric/Hmm/FullHmmTransitionMatrix.h>
 
#include "HmmLikelihood_DF.h"
#include "AlignedPhyloLikelihoodHmm.h"
 
using namespace std;
using namespace bpp;
 
/******************************************************************************/
 
AlignedPhyloLikelihoodHmm::AlignedPhyloLikelihoodHmm(
    Context& context,
    std::shared_ptr<PhyloLikelihoodContainer> pC,
    const std::vector<size_t>& nPhylo,
    bool inCollection) :
  AbstractPhyloLikelihood(context),
  AbstractParametrizable(""),
  AbstractPhyloLikelihoodSet(context, pC, {}, inCollection),
  AbstractAlignedPhyloLikelihood(context, 0),
  AbstractAlignedPhyloLikelihoodSet(context, pC, nPhylo, inCollection),
  hma_(),
  htm_(),
  hpep_(),
  hmm_()
{
  hma_ = make_shared<HmmPhyloAlphabet>(*this);
 
  htm_ = make_shared<FullHmmTransitionMatrix>(hma_, "HMM.");
 
  hpep_ = make_shared<HmmPhyloEmissionProbabilities>(hma_);
 
  hmm_ = shared_ptr<HmmLikelihood_DF>(new HmmLikelihood_DF(context, hma_, htm_, hpep_));
 
  addParameters_(htm_->getParameters());
}
 
void AlignedPhyloLikelihoodHmm::setNamespace(const std::string& nameSpace)
{
  AbstractAlignedPhyloLikelihoodSet::setNamespace(nameSpace);
  hmm_->setNamespace(nameSpace);
}
 
void AlignedPhyloLikelihoodHmm::fireParameterChanged(const ParameterList& parameters)
{
  AbstractAlignedPhyloLikelihoodSet::fireParameterChanged(parameters);
  // hpep_->update();
  hmm_->matchParametersValues(parameters);
}