avc.py 文件源码

def polyhedron(self, coords, j, L):
        """find the polyhedron for center molecule"""
        vor = Voronoi(coords)
        #get the vertices
        points = [vor.vertices[x] for x in vor.regions[vor.point_region[j]] if x != -1]
        return points