Should we consider changing the degree of freedom numbering such that
degrees of freedom at a node are located next to each other in the
matrix/vector? This will be needed for parallel assembly. It may also
speed up the assembly of sparse matrices by allowing blocks to be assembled.
