Skip to content

Commit 9241a98

Browse files
committed
Revert "Merge branch 'main' into main"
This reverts commit 6dbf95a, reversing changes made to 46f86a9.
1 parent 6dbf95a commit 9241a98

19 files changed

+380
-1856
lines changed

alignment/seqregions.cpp

+3-4
Original file line numberDiff line numberDiff line change
@@ -415,11 +415,10 @@ void cmaple::merge_identicalRACGT_TwoLowers(
415415
}
416416

417417
if (seq1_region.type == TYPE_R) {
418-
auto prev_log_lh = log_lh;
419418
log_lh += (total_blength_1 + total_blength_2) *
420419
(cumulative_rate[end_pos + 1] - cumulative_rate[pos]);
421420
} else {
422-
log_lh += model->getDiagonalMutationMatrixEntry(seq1_region.type, end_pos) *
421+
log_lh += model->diagonal_mut_mat[seq1_region.type] *
423422
(total_blength_1 + total_blength_2);
424423
}
425424
}
@@ -465,11 +464,11 @@ void cmaple::calSiteLhs_identicalRACGT(std::vector<RealNumType>&site_lh_contribu
465464
* (cumulative_rate[i + 1] - cumulative_rate[i]);
466465
}
467466
} else {
468-
log_lh += model->getDiagonalMutationMatrixEntry(seq1_region.type, pos) * total_blength;
467+
log_lh += model->diagonal_mut_mat[seq1_region.type] * total_blength;
469468

470469
// compute site lh contributions
471470
site_lh_contributions[static_cast<std::vector<RealNumType>::size_type>(pos)] +=
472-
model->getDiagonalMutationMatrixEntry(seq1_region.type, pos) * total_blength;
471+
model->diagonal_mut_mat[seq1_region.type] * total_blength;
473472
}
474473
}
475474

0 commit comments

Comments
 (0)