Matrix-free GPU implementation of a preconditioned conjugate gradient solver for anisotropic elliptic PDEs