#include <boost/geometry/formulas/spherical.hpp>
#include <boost/geometry/formulas/flattening.hpp>
-#include <boost/geometry/core/srs.hpp>
+
#include <boost/mpl/assert.hpp>
#include <boost/math/special_functions/hypot.hpp>
+ C31 * (sin2_sig3 - sin2_sig1)
+ C32 * (sin4_sig3 - sin4_sig1));
- int sign = c1;
- if (bet3 < c0)
- {
- sign = cminus1;
- }
-
+ CT const sign = bet3 >= c0
+ ? c1
+ : cminus1;
+
CT const dlon_max = omg13 - sign * f * sin_alp0 * I3;
return dlon_max;