“…It is well-known that Stokes equations describe low Reynolds number flow motion and play a fundamental role in the numerical modeling of incompressible viscous flows. Recently, there has been an increasing interest in solving Stokes problems by various meshfree (or meshless) methods [ 1 , 2 , 3 , 4 , 5 ] to alleviate mesh-related dilemmas, including some collocation meshless methods, such as virtual interpolation point method [ 6 ], generalized finite difference method [ 7 ], divergence-free kernel approximation method [ 8 ], as well as some Galerkin meshless methods, such as the moving least square reproducing kernel method [ 9 ], weighted extended B-spline method [ 10 ], Galerkin boundary node method [ 11 ], and the divergence-free meshless local Petrov–Galerkin method [ 12 ].…”