Crystal plasticity is governed by the motion of lattice dislocations. Although continuum theories of static dislocation assemblies date back to the 1950s, the line-like character of these defects posed serious problems for the development of a continuum theory of plasticity which is based on the averaged dynamics of dislocation systems. Only recently the geometrical problem of performing meaningful averages over systems of moving, oriented lines has been solved. Such averaging leads to the definition of a dislocation density tensor of second order along with its evolution equation. This tensor can be envisaged as the analogue of the classical dislocation density tensor in an extended space which includes the line orientation as an independent variable. In the current work, we discuss the numerical implementation of a continuum theory of dislocation evolution that is based on this dislocation density measure and apply this to some simple benchmark problems as well as to plane-strain micro-bending.