code: def extremum_icohpvalue(self, summed_spin_channels=True, spin=Spin.up):
          
          if not self._are_coops:
              extremum = sys.float_info.max
          else:
              extremum = -sys.float_info.max
  
          if not self._is_spin_polarized:
              if spin == Spin.down:
                  warnings.warn("This spin channel does not exist. I am switching to Spin.up")
              spin = Spin.up
  
          for value in self._icohplist.values():
              if not value.is_spin_polarized or not summed_spin_channels:
                  if not self._are_coops:
                      if value.icohpvalue(spin) < extremum:
                          extremum = value.icohpvalue(spin)
                          # print(extremum)
                  else:
                      if value.icohpvalue(spin) > extremum:
                          extremum = value.icohpvalue(spin)
                          # print(extremum)
              else:
                  if not self._are_coops:
                      if value.summed_icohp < extremum:
                          extremum = value.summed_icohp
                          # print(extremum)
                  else:
                      if value.summed_icohp > extremum:
                          extremum = value.summed_icohp
                          # print(extremum)
          return extremum
docs: get ICOHP/ICOOP of strongest bond
          Args:
              summed_spin_channels: Boolean to indicate whether the ICOHPs/ICOOPs of both spin channels should be summed
  
              spin: if summed_spin_channels is equal to False, this spin indicates which spin channel should be returned
          Returns:
              lowest ICOHP/largest ICOOP value (i.e. ICOHP/ICOOP value of strongest bond)