Before revision:  The anti-symmetric structure of the weak equation system is preserved by the matrix-vector formulation and it is utilised to prove optimal-order and uniform-in-time error estimates.
After revision:  The anti-symmetric structure of the weak equation system is preserved by the spatial discretisation. A new stability proof, based on this structure, combined with consistency bounds proves optimal-order and uniform-in-time error estimates.