  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.
