code: def energy_density(self, strain, convert_GPa_to_eV=True):
          
          return sum([c.energy_density(strain, convert_GPa_to_eV)
                      for c in self])
docs: Calculates the elastic energy density due to a strain