Bond order efficiency fix

Instead of operating on n_atoms size arrays of gradients, forces due to
bond order gradients are now calculated and stored directly when
gradients are evaluated.