Although they do not lend themselves to analytical resolution, three-body atomic or molecular systems are still simple enough to allow for very precise theoretical predictions of their energy levels, which makes them attractive candidates for fundamental tests and determination of fundamental physical constants. Focusing on the hydrogen molecular ions (H + 2 , HD + , D + 2 ), we outline the methods which have been used to improve the theoretical accuracy by several orders of magnitude over the last two decades. The three-body Schrödinger equation can be solved with extreme precision by variational methods with trial functions involving exponentials of interparticle distances. Quantum electrodynamics (QED) corrections are evaluated in the framework of nonrelativistic QED (NRQED). The current status of theory and possibilities of further improvement are briefly sketched.