Note: We simply switched the linear solver for solving the small size linear systems in topology steps from sparse to dense,
which makes some of the results slightly different due to the numerical differences.
However, we achieve 3x speedup with 2% shorter seams in average.