code: def parse_dim(features, check=True):
      
      # try:
      dim = features[0].shape[1]
      # except IndexError:
      #     dim = 1
  
      if check and not dim > 0:
          raise IOError('features dimension must be strictly positive')
      if check and not all([d == dim for d in [x.shape[1] for x in features]]):
          raise IOError('all files must have the same feature dimension')
      return dim
docs: Return the features dimension, raise if error
  
      Raise IOError if features have not all the same positive
      dimension.  Return dim (int), the features dimension.