Scalable robust solvers for unstructured FE geodynamic modeling applications: Solving the Stokes equation for models with large localized viscosity contrasts: A SCALABLE ROBUST SOLVER FOR THE STOKES EQUATION