This is a fairly widely-used 2D code for Voronoi diagrams and Delauney triangulations, written in C by Steve Fortune of Bell Laboratories. It is based on Fortune's sweepline algorithm for Voronoi diagrams, and is likely to be the right code to try first.