A lattice Boltzmann formulation for relativistic fluids is presented and numerically validated through quantitative comparison with recent hydrodynamic simulations of relativistic fluids. In order to illustrate its capability to handle complex geometries, the scheme is also applied to the case of a three-dimensional relativistic shock wave, generated by a supernova explosion, impacting on a massive interstellar cloud. This formulation opens up the possibility of exporting the proven advantages of lattice Boltzmann methods, namely, computational efficiency and easy handling of complex geometries, to the context of (mildly) relativistic fluid dynamics at large, from quark-gluon plasmas up to supernovae with relativistic outflows.