more comments fixed

imops/cpp/interp2d/delaunator/delaunator.cpp

@@ -91,9 +91,9 @@ inline bool counterclockwise(const Point& p0, const Point& p1, const Point& p2)
     Point v0 = Point::vector(p0, p1);
     Point v1 = Point::vector(p0, p2);
     double det = Point::determinant(v0, v1);
-    double dist = v0.magnitude2() + v1.magnitude2();
     if (det == 0)
         return false;
+    double dist = v0.magnitude2() + v1.magnitude2();
     double reldet = std::abs(dist / det);
     if (reldet > 1e14)
         return false;
@@ -127,6 +127,8 @@ inline Point circumcenter(
     //ABELL - This is suspect for div-by-0.
     const double d = dx * ey - dy * ex;
 
+    assert(("Divide by zero exception", d == 0.0));
+
     const double x = ax + (ey * bl - dy * cl) * 0.5 / d;
     const double y = ay + (dx * cl - ex * bl) * 0.5 / d;
 

imops/cpp/main.cpp

@@ -1,3 +1,4 @@
+#include <cassert>
 #include <iostream>
 #include <stdexcept>
 #include "interp2d/interpolator.h"

imops/interp2d.py

@@ -58,6 +58,12 @@ def __init__(
 
         if not isinstance(points, np.ndarray):
             raise ValueError(f'Wrong type of `points` argument, expected np.ndarray. Got {type(points)}')
+
+        if not isinstance(self.values, np.ndarray):
+            raise ValueError(f'Wrong type of `values` argument, expected np.ndarray. Got {type(self.values)}')
+
+        if self.values.ndim > 1:
+            raise ValueError(f'Wrong shape of `values` argument, expected ndim=1. Got shape {self.values.shape}')
 
         super().__init__(points, num_threads, triangles)
         self.kdtree = KDTree(data=points, **kwargs)