MurrellGroup / MolecularEvolution.jl

A Julia framework for developing phylogenetic models
MIT License
11 stars 4 forks source link

felsenstein_down! on polytomies #29

Open nossleinad opened 3 months ago

nossleinad commented 3 months ago

felsenstein_down! uses a nested loop over children and siblings for combine! ops. This is fine for binary trees but very slow (quadratic time complexity, with respect to children) for polytomies. Since a polytomy is equivalent to a binary tree with some added internal nodes with 0 branchlength, we can improve here.