Changeset 18912 for issm/trunk-jpl/src/c/classes/kriging/Observation.cpp

Timestamp:

12/03/14 11:37:24 (10 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added covertree nearest neighbor search (to be debugged)

File:

• issm/trunk-jpl/src/c/classes/kriging/Observation.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/kriging/Observation.cpp

@@ -4 +4 @@
 
 #include <stdlib.h>
+#include <cmath>
+#include <utility>
 #include "../classes.h"
 
@@ -23 +25 @@
 }
 /*}}}*/
+Observation::Observation(double x_in, double y_in,double value_in){
+        this->x = x_in;
+        this->y = y_in;
+        this->value = value_in;
+
+        this->xi     = 0;
+        this->yi     = 0;
+        this->index  = 0;
+        this->weight = 0.;
+}
 Observation::~Observation(){/*{{{*/
         return;
@@ -59 +71 @@
 }
 /*}}}*/
+
+/*Covertree*/
+void Observation::print(void) const{/*{{{*/
+
+        _printf_("Observation\n");
+        _printf_("   x     : " << this->x << "\n");
+        _printf_("   y     : " << this->y << "\n");
+        _printf_("   value : " << this->value << "\n");
+}
+/*}}}*/
+double Observation::distance(const Observation& ob) const
+{
+        return std::sqrt( (std::pow( (ob.getXY().first - this->x), 2 ) + std::pow((ob.getXY().second - this->y), 2) ));
+}
+
+const std::pair<double, double>& Observation::getXY() const
+{
+        return std::make_pair(this->x, this->y);
+}
+bool Observation::operator==(const Observation& ob) const
+{
+        return (ob.getXY().first == this->x && ob.getXY().second == this->y && ob.value == this->value);
+}
+
+void Observation::WriteXYObs(const Observation& ob, double* px, double* py, double* pobs){
+    *px   = ob.getXY().first;
+    *py   = ob.getXY().second;
+    *pobs = ob.value;
+}