add convergence test of calculate_s2m()

add a variable geometry_residence to test the residence of each round, this variable is to test the convergence