Gravitational wave radiation is only unambiguously defined at future null infinity -- the "location" where global properties of spacetimes can be measured, and also the correct idealization of astrophysical observers. Reaching future null infinity is thus crucial for extracting correct waveforms from numerical relativity simulations of compact binaries, and evolution on hyperboloidal slices -- spacelike, smooth and extending to null infinity -- is a very promising avenue to achieve this. Among current efforts to the hyperboloidal method, I will focus on free evolution of the conformally compactified BSSN / Z4 equations. In this talk I will cover hyperboloidal initial data including black holes and neutron stars in spherical symmetry and ongoing efforts towards simulations in full 3D in this approach.