boost/geometry
2
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2026-02-01 08:32:10 +00:00
Code Issues Projects Releases Wiki Activity

[algorithms] Add iterations counter to vincenty_inverse for robustness.

Browse Source
This commit is contained in:
Adam Wulkiewicz
2014-12-06 19:34:56 +01:00
parent e3a7bce127
commit ece229983f
1 changed files with 11 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
include/boost/geometry/algorithms/detail/vincenty_inverse.hpp
View File

@@ -25,6 +25,11 @@
#include <boost/geometry/algorithms/detail/flattening.hpp>
#ifndef BOOST_GEOMETRY_DETAIL_VINCENTY_MAX_STEPS
#define BOOST_GEOMETRY_DETAIL_VINCENTY_MAX_STEPS 1000
#endif
namespace boost { namespace geometry { namespace detail
{
@@ -103,6 +108,8 @@ public:
CT previous_lambda;
int counter = 0; // robustness
do
{
previous_lambda = lambda; // (13)
@@ -119,8 +126,11 @@ public:
lambda = L + (c1 - C) * flattening * sin_alpha *
(sigma + C * sin_sigma * ( cos2_sigma_m + C * cos_sigma * (-c1 + c2 * math::sqr(cos2_sigma_m)))); // (11)
++counter; // robustness
} while ( geometry::math::abs(previous_lambda - lambda) > c_e_12
&& geometry::math::abs(lambda) < pi );
&& geometry::math::abs(lambda) < pi
&& counter < BOOST_GEOMETRY_DETAIL_VINCENTY_MAX_STEPS ); // robustness
}
inline CT distance() const