fix P matrix calculation for an array of t for GTR

jangevaare · Jul 27, 2017 · da38079 · da38079
1 parent 914acd5
commit da38079
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/substitution_models/gtr.jl b/src/substitution_models/gtr.jl
@@ -83,7 +83,7 @@ function P(gtr::GTR, t::Array{Float64})
     error("t must be positive")
   end
   eig_vals, eig_vecs = eig(Q(mod))
-  return [expm(eig_vecs * (diagm(eig_vals)*i) * eig_vecs') for i in t]
+  return [eig_vecs * expm(diagm(eig_vals)*i) * eig_vecs' for i in t]
 end