JCnuc.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: The Bio++ Development Group
//
// SPDX-License-Identifier: CECILL-2.1
 
#ifndef BPP_PHYL_MODEL_NUCLEOTIDE_JCNUC_H
#define BPP_PHYL_MODEL_NUCLEOTIDE_JCNUC_H
 
 
#include "../AbstractSubstitutionModel.h"
#include "NucleotideSubstitutionModel.h"
 
namespace bpp
{
class JCnuc :
  public AbstractReversibleNucleotideSubstitutionModel
{
private:
  mutable double exp_;
  mutable RowMatrix<double> p_;
 
public:
  JCnuc(std::shared_ptr<const NucleicAlphabet> alpha);
 
  virtual ~JCnuc() {}
 
  JCnuc* clone() const override { return new JCnuc(*this); }
 
public:
  double Pij_t    (size_t i, size_t j, double d) const override;
  double dPij_dt  (size_t i, size_t j, double d) const override;
  double d2Pij_dt2(size_t i, size_t j, double d) const override;
  const Matrix<double>& getPij_t    (double d) const override;
  const Matrix<double>& getdPij_dt  (double d) const override;
  const Matrix<double>& getd2Pij_dt2(double d) const override;
 
  std::string getName() const override { return "JC69"; }
 
  void setFreqFromData(const SequenceDataInterface& data, double pseudoCount = 0) override {}
 
protected:
  void updateMatrices_() override;
};
} // end of namespace bpp.
#endif // BPP_PHYL_MODEL_NUCLEOTIDE_JCNUC_H