Vornoi tessellations

Hi everyone,

I’m working on a crystal plasticity simulation for a polycrystal and need to generate a Voronoi tessellation. Does anyone know how I can create Voronoi diagrams using FreeFEM++?

Thanks for any advice or suggestions!