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.