The previous version used the same b_tot^2 for all elements, effectively
splitting the sum into the product of two sums. It was not noticed
before because the unit test was for a single element material, Si,
where the sum only has one value.