log_double_t correction(const data_partition& P,const vector<int>& nodes) { if (P.variable_alignment()) { // get the lengths of then internal node int length = P.seqlength(nodes[0]); return pow(P.sequence_length_pr(length), 2); } else return 1; }