GTRRateMatrix::GTRRateMatrix(StateMap const & staMap, vector_t const & rateParameters, vector_t const & equiFreqs) : BaseTRRateMatrix(rateParameters, equiFreqs, "GTR"), staMap_(staMap) { assert( (staMap.stateCount() * ( staMap.stateCount() - 1) ) / 2 == rateParameterCount_); assert(staMap.stateCount() == equiFreqCount_); }
GTRRateMatrix::GTRRateMatrix(StateMap const & staMap) : BaseTRRateMatrix( (staMap.stateCount() * ( staMap.stateCount() - 1) ) / 2, staMap.stateCount(), "GTR"), staMap_(staMap) {}