These codes recover some of the results in paper PHYSICAL REVIEW X 11, 011030 (2021) ================= Please run the Live Script: note_main_factorizable_ensemble.mlx with Matlab version later than R2016b Using parallel is preferred, otherwise, replace all 'parfor' by 'for'. Or read the note without Matlab: note_main_factorizable_ensemble.pdf