Tests and minor changes. One change was made to increase clocklikeness when the substitutions rate is constant. the function is working well but is slow.