Oracles + python

include/cartesian_geom/point.h

@@ -10,13 +10,15 @@
 #define POINT_H
 
 #include <iostream>
+#include <cmath>
+#include <vector>
 
 template <class K>
 class point
 {
 private:
     unsigned int d;
-    typedef std::vector<typename K::FT> Coeff;
+    typedef typename std::vector<typename K::FT> Coeff;
     Coeff coeffs;
     typedef typename std::vector<typename K::FT>::iterator iter;
 public:
@@ -33,13 +35,30 @@ class point
         d = dim;
         coeffs = Coeff(begin,endit);
     }
+
+    point(const unsigned int dim, FT* data) {
+        d = dim;
+        coeffs.reserve(d);
+        for (uint i=0; i<dim; i++) {
+            coeffs[i] = data[i];
+        }
+    }
+
+    FT* data() {
+        return coeffs.data();
+    }
+
+    Coeff& get_coeffs() {
+        return coeffs;
+    }
 
     int dimension() {
         return d;
     }
 
     void set_dimension(const unsigned int dim) {
         d = dim;
+        coeffs.reserve(d);
     }
 
     void set_coord(const unsigned int i, FT coord) {
@@ -76,6 +95,16 @@ class point
         return temp;
     }
 
+    double squared_distance(point& p) {
+
+        double dist = 0.0;
+        for (auto mit = coeffs.begin(), pit=p.iter_begin(); pit < p.iter_end(); ++pit, ++mit) {
+            double tmp = (*mit)-(*pit);
+            dist += tmp*tmp;
+        }
+        return dist;
+    }
+
     point operator* (const FT& k) {
         point temp(d, iter_begin(), iter_end());
 
@@ -121,7 +150,13 @@ class point
         std::cout<<"\n";
         #endif
     }
-
+
+    void normalize() {
+        FT lsq=std::sqrt(squared_length());
+        for (auto mit=coeffs.begin(); mit!=coeffs.end(); mit++){
+            (*mit) = (*mit)/lsq;
+        }
+    }
 
     iter iter_begin() {
         return coeffs.begin();