I think your implementation of logsumexp is incorrect (see my comment on Github https://github.com/TuringLang/Turing.jl/issues/1374#issuecomment-673375285).
BTW I’m a bit confused where to reply, just posting here as well in case you missed my reply on Github.