F-step is numerically hard to solve. I set BarHomogeneous = 1 to get a suboptimal solution
There are still some numerical problem to solve, but the alternation starts to work now. The numerical problems include
1. Ill scaling of the solution. When the initial seed is bad, the angular momentum is huge, while all other entries in the solution are about 1
2. Using multiple threads will generate different result in GUROBI