Having calculated the total energy, both the density-kernel
and the expansion coefficients for the localised orbitals
are varied. Because of the non-orthogonality
of the support functions, it is necessary to take note of the tensor
properties of the gradients [163], as noted in section 4.6.