Petascale elliptic solvers for anisotropic PDEs on GPU clusters