A really fast C++ library for Delaunay triangulation of 2D points.
delaunator-cpp is a C++ port from https://github.com/mapbox/delaunator a JavaScript implementation of very fast 2D Delaunay algorithm.
- Probably the fastest C++ open source 2D Delaunay implementation
- Example showing triangulation of GeoJson points
examples/basic.cpp
#include <delaunator.hpp>
#include <cstdio>
int main() {
/* x0, y0, x1, y1, ... */
std::vector<double> coords = {-1, 1, 1, 1, 1, -1, -1, -1};
//triangulation happens here
delaunator::Delaunator d(coords);
for(std::size_t i = 0; i < d.triangles.size(); i+=3) {
printf(
"Triangle points: [[%f, %f], [%f, %f], [%f, %f]]\n",
d.coords[2 * d.triangles[i]], //tx0
d.coords[2 * d.triangles[i] + 1], //ty0
d.coords[2 * d.triangles[i + 1]], //tx1
d.coords[2 * d.triangles[i + 1] + 1],//ty1
d.coords[2 * d.triangles[i + 2]], //tx2
d.coords[2 * d.triangles[i + 2] + 1] //ty2
);
}
}
Run on (4 X 2300 MHz CPU s)
2018-09-19 09:40:29
------------------------------------------------------------
Benchmark Time CPU Iterations
------------------------------------------------------------
BM_45K_geojson_nodes 23 ms 23 ms 32
BM_2K_uniform 1 ms 1 ms 970
BM_100K_uniform 63 ms 63 ms 10
BM_200K_uniform 142 ms 141 ms 4
BM_500K_uniform 415 ms 412 ms 2
BM_1M_uniform 1016 ms 1010 ms 1
Library is ~10% faster then JS version for 1M uniform points (details)