AlphaPhi solver works, and gets the RIGHT ANSWER... which is weird, because it looks like the likelihood is a sum of nonconvex things.
But it works anyways!
No clue why, no clue if this generalizes to real data.