code: def is_coord_subset(subset, superset, atol=1e-8):
      
      c1 = np.array(subset)
      c2 = np.array(superset)
      is_close = np.all(np.abs(c1[:, None, :] - c2[None, :, :]) < atol, axis=-1)
      any_close = np.any(is_close, axis=-1)
      return np.all(any_close)
docs: Tests if all coords in subset are contained in superset.
      Doesn't use periodic boundary conditions
  
      Args:
          subset, superset: List of coords
  
      Returns:
          True if all of subset is in superset.