Merge from ReForm: Optimize MISolver::solve_multiple_rounding()
- Improve the performance for dense meshes of the solve_multiple_rounding() which O(n^2) by using more efficient data structures and processing.
- Clean up the code, improve formatting and const correctness.