Implementation of a multigrid solver on GPU for Stokes equations with strongly variable viscosity based on Matlab and CUDA