Numerical simulations of two-phase flows driven by viscosity (e.g. for bubble motions in glass melting process) rely on the ability to efficiently compute the solutions to discretized Stokes equations. When using boundary element methods to track fluid interfaces, one usually faces the problem of solving linear systems with a dense matrix with a size proportional to the system number of degrees...
