docs: def normalize_to(self, comp, factor=1):
          """
          Normalizes the reaction to one of the compositions.
          By default, normalizes such that the composition given has a
          coefficient of 1. Another factor can be specified.
  
          Args:
              comp (Composition): Composition to normalize to
              factor (float): Factor to normalize to. Defaults to 1.
          """
          scale_factor = abs(1 / self._coeffs[self._all_comp.index(comp)]
                             * factor)
          self._coeffs = [c * scale_factor for c in self._coeffs]