Add euclidean distance computation to point

diff --git a/incl/bcar.hh b/incl/bcar.hh
index 3cef171b475cb1bd381c0aee62445b2f3044879e..a84053ae6320722d4a89936c5a1397f7d94afca7 100644 (file)
--- a/incl/bcar.hh
+++ b/incl/bcar.hh
@@ -65,6 +65,9 @@ public:
        */
        void rotate(Point const& c, double const angl);
 
+       /*! Return Euclidean distance to `p`. */
+       double edist(Point const& p) const;
+
        friend std::ostream& operator<<(std::ostream& out, Point const& p);
 };
 

diff --git a/src/bcar.cc b/src/bcar.cc
index 7c8dd93f6fab1da00e9f590890a911ef1c14757e..54d30b9ebd40e1912d2ac49e3db1ec885deee151 100644 (file)
--- a/src/bcar.cc
+++ b/src/bcar.cc
@@ -108,6 +108,12 @@ Point::rotate(Point const& c, double const angl)
        this->y(ny + c.y());
 }
 
+double
+Point::edist(Point const& p) const
+{
+       return sqrt(pow(p.x() - this->x_, 2.0) + pow(p.y() - this->y_, 2.0));
+}
+
 std::ostream&
 operator<<(std::ostream& out, Point const& p)
 {