The problem of viscoelastic relaxation of the Earth solved by a matrix eigenvalue approach based on discretization in grid space