  In this paper we study the statistical properties of a reversible cellular
automaton in two out-of-equilibrium settings. In the first part we consider two
instances of the initial value problem, corresponding to the inhomogeneous
quench and the local quench. Our main result is an exact matrix product
expression of the time evolution of the probability distribution, which we use
to determine the time evolution of the density profiles analytically. In the
second part we study the model on a finite lattice coupled with stochastic
boundaries. Once again we derive an exact matrix product expression of the
stationary distribution, as well as the particle current and density profiles
in the stationary state. The exact expressions reveal the existence of
different phases with either ballistic or diffusive transport depending on the
boundary parameters.
