Minor changes in Kwak's optimization
Main simulation 3