nprodHHMM.h

#ifndef nprodHHMM_H
#define nprodHHMM_H

#include <vector>
#include <cstdint>
#include <memory>
#include <Eigen/Dense>
#include <map>
#include "baseHHMM.hpp"
#include "prodHHMM.hpp"
#include "myalgorithm.hpp"

using namespace std;
using namespace Eigen;

namespace hhmm{

  template <typename T>
  using up = unique_ptr<T>;

  class TestHHMM;

  class nprodHHMM final:public baseHHMM{
    friend TestHHMM;
  protected:
    MatrixXld transMat;
  public:
    vector<up<baseHHMM>> children;
    map<uint64_t,uint32_t> convert;

    nprodHHMM(vector<uint32_t> const&,uint32_t);
    nprodHHMM(vector<uint32_t> const&,uint32_t,uint32_t,baseHHMM*);
    ~nprodHHMM() noexcept {};
    MatrixXld& trans();
    long double trans(up<baseHHMM> const&);
    long double trans(baseHHMM*);
    long double trans(const up<baseHHMM>&,const up<baseHHMM>&);
    long double trans(baseHHMM*,const up<baseHHMM>&);
    long double trans(const up<baseHHMM>&,baseHHMM*);
    void cpyTransMat(MatrixXld const&);
    long double& setTrans(baseHHMM*,const up<baseHHMM>&);
    void clearParam();
    void initParam(vector<long double> const&);
    void check();
    void log(uint32_t,uint32_t,string const&);
  };
}

#endif