This is the code accompanying the paper A concise representation for adaptive hexagonal meshes, submitted for publication (2024).
Adaptive hexagonal meshes were introduced in the paper Hexagonal LOD for interactive terrain rendering by Sußner et al. (2005). See also Hexagonal Delaunay triangulation by Sußner and Greiner (2009) and Hexagon-based all-quadrilateral mesh generation with guaranteed angle bounds by Liang and Zhang (2011).
Running make outputs a sample adaptive hexagonal mesh in eps and csv. The eps file is converted to pdf using pstopdf. Choose the level of refinement by setting R, as in make R=6. Compute a sequence of increasingly refined meshes with sh GRSI.sh.
The file in.csv is a sample adaptive hexagonal mesh generated by random refinement. It can be loaded with make I=in.csv.
The code should work in both Python 2 and Python 3.