Given a configuration of points, a procedure for constructing the corresponding Voronoi diagram is given. The procedure is exact for molecules in the bulk. Polyhedra of surface molecules can be either eliminated or included using a periodic boundary condition. The construction is of interest in astronomy, biology, chemistry, materials science, as well as in physics (with points representing atoms, molecules, ions, etc.). The present method is more efficient than other procedures described in the literature.