The two polymorphs of TiO2, rutile and anatase, have been investigated at the ab initio level using different Hamiltonians with all-electron Gaussian and projector augmented plane wave basis sets. Their equilibrium lattice parameters, relative stabilities, binding energies, and band structures have been evaluated. The calculations have been performed at the Hartree-Fock, density functional theory (DFT), and hybrid (B3LYP and PBE0) levels. As regards DFT, the local density and generalized gradient (PBE) approximations have been used. Our results show an excellent agreement with the experimental band structures and binding energies for the B3LYP and PBE0 functionals, while the best structural descriptions are obtained at the PBE0 level. However, no matter which Hamiltonian and method are used, anatase is found more stable than rutile, in contrast with recent experimental reports, although the relative stabilities of the two phases are very close to each other. Nevertheless, based on the overall results, the hybrid PBE0 functional appears as a good compromise to obtain an accurate description of both structural and electronic properties of solids.