As is well known, a working MG has three main components: (i) Smoothing via an iterative method; (ii) Restriction from a fine grid to a coarse grid; (iii) Interpolation from a coarse grid to a fine one. On the coarsest grid, an effective unilevel solver is used for accurate solution; here, we shall use the semi-implicit method represented by Equation (19). Without reducing the importance of the restriction and interpolation operators, the efficiency of every MG method strongly relies on the efficiency of the relaxation method, also known as smoother, used at each level in reducing