I used a Voronoi function instead of "Divide surface" in order to get the hexagonal tiling. 

also instead of curve connect function, Beziers curves were used. it is more complex, but G2 connection makes the difference. curvature continuity more beautiful than Tangent continuity.

let me know if any interest for the algorithm.

data inside

​​​​​​​