code: def find_pore_hulls(self, pores=None):
          r
          if pores is None:
              pores = self.pores('delaunay')
          temp = []
          tvals = self['throat.interconnect'].astype(int)
          am = self.create_adjacency_matrix(weights=tvals, fmt='lil',
                                            drop_zeros=True)
          for p in pores:
              Ps = am.rows[p]
              temp.append(Ps)
          return sp.array(temp, dtype=object)
docs: r"""
          Finds the indices of the Voronoi nodes that define the convex hull
          around the given Delaunay nodes.
  
          Parameters
          ----------
          pores : array_like
              The pores whose convex hull are sought.  The given pores should be
              from the 'delaunay' network.  If no pores are given, then the hull
              is found for all 'delaunay' pores.
  
          Notes
          -----
          This metod is not fully optimized as it scans through each pore in a
          for-loop, so could be slow for large networks.