  An error estimate for a canonical discretization of the harmonic map heat
flow into spheres is derived. The numerical scheme uses standard finite
elements with a nodal treatment of linearized unit-length constraints. The
analysis is based on elementary approximation results and only uses the
discrete weak formulation.
