docs: def find_atomics(formula: Formula) -> Set[PLAtomic]:
      """Finds all the atomic formulas"""
      f = formula
      res = set()
      if isinstance(formula, PLFormula):
          res = formula.find_atomics()
      # elif isinstance(f, PLNot):
      #     res = res.union(find_atomics(f.f))
      # elif isinstance(f, PLBinaryOperator):
      #     for subf in f.formulas:
      #         res = res.union(find_atomics(subf))
      else:
          res.add(f)
      return res