Iterative multiscale gradient computation for heterogeneous subsurface flow