An implementation of the Bowyer-Watson for computing Delaunay trianglulations and Voronoi diagrams