Many different simulation methods for Stokes flow problems involve a common computationally intense task -- the summation of kernel function over $O(N^2)$ pairs points. One popular technique is Kernel Independent Fast Multipole Method (KIFMM), which constructs spatial adaptive octree all points and places small number equivalent multipole local around each box, completes sum with $O(N)$ cost, u...