The phenomenon of quasi-periodicity in deterministic dynamical systems describes stationary solutions, which neither exhibit a finite period length nor are chaotic. Recently, an increasing demand for robust numerical methods is driven by applied dynamics and industrial applications. In this context, direct time integration proves to be impractical due to extensive integration intervals. Therefore, in a first step, this contribution aims on giving an application oriented survey of the basic theory as well as alternative concepts. In the following, the focus is set on the direct computation of invariant manifolds (surfaces) on which quasi-periodic solutions evolve. This approach offers a unique framework from which classical methods (e.g., the multi-harmonic-balance) can be systematically deduced and mutual similarities between different methods may be revealed. This contribution starts with a brief summary of related mathematical basics, which is followed by an overview of available methods. Subsequently, the computation of invariant manifolds by means of solving a partial differential equation is emphasized. These PDEs may be formulated using different parametrization strategies. Here, the concept of hyper-time parametrization is particularly interesting, since it is a promising starting point for the development of numerical schemes with general applicability in engineering problems. In order to solve the underlying PDE, various methods may be used. The implementation of a Fourier-Galerkin method as well as a finite difference method is presented and compared on the basis of computational results of the van-der-Pol equation (with and without forcing). Moreover, it is demonstrated that both methods apply to periodic as well as quasi-periodic solutions alike. In order to exemplify the practical use, these methods are applied to a generic rotordynamic model problem.