From e4b057a444c814eab35128e12508d30c2c8e9e02 Mon Sep 17 00:00:00 2001
From: Cora Allen-Coleman <6508099+coraallencoleman@users.noreply.github.com>
Date: Tue, 28 Apr 2020 15:21:25 -0500
Subject: [PATCH] array bounds debug

---
 src/phyLiNCoptimization.jl | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/phyLiNCoptimization.jl b/src/phyLiNCoptimization.jl
index 2c752a925..9976c5fbe 100644
--- a/src/phyLiNCoptimization.jl
+++ b/src/phyLiNCoptimization.jl
@@ -712,6 +712,9 @@ function addhybridedgeLiNC!(obj::SSM, currLik::Float64, maxhybrid::Int,
     # fixroot=true: to restore edge2 if need be, with deletehybridedge!
     isnothing(result) && return nothing
     newhybridnode, newhybridedge = result
+    if maximum([e.number for e in obj.net.edge]) > size(obj.logtrans)[3]
+        @debug "The new hybrid edge number is greater than the allowed maxedge in logtrans"
+    end
     # unzip only at new node and its child edge
     unzipat_canonical!(newhybridnode, getChildEdge(newhybridnode))
     updateSSM!(obj, true; constraints=constraints)