We study the relativistic heat equation in one space dimension. We prove a local regularity result when the initial datum is locally Lipschitz in its support. We propose a numerical scheme that captures the known features of the solutions and allows for analysing further properties of their qualitative behavior.