From cd3f264a275fc72e8f7758573336ef384d4bfd89 Mon Sep 17 00:00:00 2001
From: Vissarion Fysikopoulos <vissarion.fysikopoulos@oracle.com>
Date: Thu, 31 Aug 2017 11:02:55 +0300
Subject: [PATCH] [formulas] [distance] Activate differential quantities
 computation in distance point segment

---
 include/boost/geometry/formulas/distance_point_segment.hpp | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/include/boost/geometry/formulas/distance_point_segment.hpp b/include/boost/geometry/formulas/distance_point_segment.hpp
index a2f32cd2a..afa4eebd7 100644
--- a/include/boost/geometry/formulas/distance_point_segment.hpp
+++ b/include/boost/geometry/formulas/distance_point_segment.hpp
@@ -25,8 +25,6 @@
       https://arxiv.org/abs/1102.1215
 */
 
-#define BOOST_GEOMETRY_DISTANCE_POINT_SEGMENT_DEBUG
-
 namespace boost { namespace geometry { namespace formula
 {
 
@@ -42,7 +40,7 @@ class distance_point_segment{
 
 public:
 
-    typedef Inverse<CT, true, false, false, false, false> inverse_distance_type;
+    typedef Inverse<CT, true, false, false, true, true> inverse_distance_type;
     typedef Inverse<CT, false, true, false, false, false> inverse_azimuth_type;
     typedef Direct<CT, true, false, false, false> direct_distance_type;
 
@@ -266,6 +264,9 @@ public:
             std::cout << "delta_g4=" << delta_g4  << std::endl;
             std::cout << "g4=" << g4 * math::r2d<CT>() << std::endl;
             std::cout << "der=" << der  << std::endl;
+            std::cout << "M43=" << M43 << std::endl;
+            std::cout << "spherical limit=" << cos(s14/earth_radius) << std::endl;
+            std::cout << "m34=" << m34 << std::endl;
             std::cout << "new_s14=" << s14 << std::endl;
             std::cout << std::setprecision(16) << "dist     =" << res34.distance << std::endl;
             std::cout << "---------end of step " << counter << std::endl<< std::endl;