add test to Delaunay

code/Delaunay.cc → code/Delaunay.cpp

@@ -35,7 +35,7 @@ vector<triple> delaunayTriangulation(vector<T>& x, vector<T>& y) {
   }
   return ret;
 }
-int main()
+int test()
 {
   T xs[]={0, 0, 1, 0.9};
   T ys[]={0, 1, 0, 0.9};

code/Delaunay_test.cpp

@@ -0,0 +1,18 @@
+#include <bits/stdc++.h>
+using namespace std;
+
+#include "Delaunay.cpp"
+
+int main() {
+  T xs[]={0, 0, 1, 0.9};
+  T ys[]={0, 1, 0, 0.9};
+  vector<T> x(&xs[0], &xs[4]), y(&ys[0], &ys[4]);
+  vector<triple> tri = delaunayTriangulation(x, y);
+  assert(tri[0].i == 0);
+  assert(tri[0].j == 1);
+  assert(tri[0].k == 3);
+  assert(tri[1].i == 0);
+  assert(tri[1].j == 3);
+  assert(tri[1].k == 2);
+  return 0;
+}