Question about Voronoi Mesh Implementation in FreeFEM++ (2D)

Dear FreeFEM++ Community,

I hope you are doing well.

I am currently working on a 2D numerical simulation and am interested in using a Voronoi-type mesh for my HHO problem. As far as I understand, FreeFEM++ primarily supports triangular (Delaunay) meshes in 2D, and there is no direct built-in functionality for generating Voronoi meshes.

I would like to ask for your guidance on the following points:

  1. Is there any recommended approach within FreeFEM++ to construct or approximate a Voronoi mesh in 2D and 3D?.

  2. Has anyone implemented the dual mesh (Voronoi diagram) from a Delaunay triangulation inside FreeFEM++?. If so, could you share an example or outline the method?.

  3. Are there any scripts, plugins, or workflows (possibly involving external tools) that can generate Voronoi meshes and import them into FreeFEM++?.

If anyone has sample code, references, or suggestions, I would greatly appreciate your help.

Thank you very much for your time and support.

With best regards,
Monirul