code: def computeDelaunayTriangulation(points):
      
      siteList = SiteList(points)
      context  = Context()
      context.triangulate = True
      voronoi(siteList,context)
      return context.triangles
docs: Takes a list of point objects (which must have x and y fields).
          Returns a list of 3-tuples: the indices of the points that form a
          Delaunay triangle.