docs: def contains(ell, p, shell_only=False):
          """
          Check to see whether point is inside
          conic.
  
          :param exact: Only solutions exactly on conic
            are considered (default: False).
          """
          v = augment(p)
          _ = ell.solve(v)
          return N.allclose(_,0) if shell_only else _ <= 0