code: def get_max_bond_lengths(structure, el_radius_updates=None):
      
      #jmc = JMolCoordFinder(el_radius_updates)
      jmnn = JmolNN(el_radius_updates=el_radius_updates)
  
      bonds_lens = {}
      els = sorted(structure.composition.elements, key=lambda x: x.Z)
  
      for i1 in range(len(els)):
          for i2 in range(len(els) - i1):
              bonds_lens[els[i1], els[i1 + i2]] = jmnn.get_max_bond_distance(
                  els[i1].symbol, els[i1 + i2].symbol)
  
      return bonds_lens
docs: Provides max bond length estimates for a structure based on the JMol
      table and algorithms.
      
      Args:
          structure: (structure)
          el_radius_updates: (dict) symbol->float to update atomic radii
      
      Returns: (dict) - (Element1, Element2) -> float. The two elements are 
          ordered by Z.