code: def find_throat_facets(self, throats=None):
          r
          if throats is None:
              throats = self.throats('delaunay')
          temp = []
          tvals = self['throat.interconnect'].astype(int)
          am = self.create_adjacency_matrix(weights=tvals, fmt='lil',
                                            drop_zeros=True)
          for t in throats:
              P12 = self['throat.conns'][t]
              Ps = list(set(am.rows[P12][0]).intersection(am.rows[P12][1]))
              temp.append(Ps)
          return sp.array(temp, dtype=object)
docs: r"""
          Finds the indicies of the Voronoi nodes that define the facet or
          ridge between the Delaunay nodes connected by the given throat.
  
          Parameters
          ----------
          throats : array_like
              The throats whose facets are sought.  The given throats should be
              from the 'delaunay' network. If no throats are specified, all
              'delaunay' throats are assumed.
  
          Notes
          -----
          The method is not well optimized as it scans through each given throat
          inside a for-loop, so it could be slow for large networks.