Implementation of the Bowyer-Watson algorithm for computing Delaunay trianglulations and Voronoi diagrams