Add
  calc_elastic_constants(pot, at)
