123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // This file is manually converted from PROJ4
- // Copyright (c) 2008-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // This file was modified by Oracle on 2017, 2018.
- // Modifications copyright (c) 2017-2018, Oracle and/or its affiliates.
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- // This file is converted from PROJ4, http://trac.osgeo.org/proj
- // PROJ4 is originally written by Gerald Evenden (then of the USGS)
- // PROJ4 is maintained by Frank Warmerdam
- // PROJ4 is converted to Geometry Library by Barend Gehrels (Geodan, Amsterdam)
- // Original copyright notice:
- // Permission is hereby granted, free of charge, to any person obtaining a
- // copy of this software and associated documentation files (the "Software"),
- // to deal in the Software without restriction, including without limitation
- // the rights to use, copy, modify, merge, publish, distribute, sublicense,
- // and/or sell copies of the Software, and to permit persons to whom the
- // Software is furnished to do so, subject to the following conditions:
- // The above copyright notice and this permission notice shall be included
- // in all copies or substantial portions of the Software.
- // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
- // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
- // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
- // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
- // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
- // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
- // DEALINGS IN THE SOFTWARE.
- #ifndef BOOST_GEOMETRY_PROJECTIONS_PHI2_HPP
- #define BOOST_GEOMETRY_PROJECTIONS_PHI2_HPP
- #include <boost/geometry/srs/projections/exception.hpp>
- #include <boost/geometry/srs/projections/impl/pj_strerrno.hpp>
- #include <boost/geometry/util/math.hpp>
- namespace boost { namespace geometry { namespace projections {
- namespace detail {
- template <typename T>
- inline T pj_phi2(T const& ts, T const& e)
- {
- static const T TOL = 1.0e-10;
- static const int N_ITER = 15;
- T eccnth, Phi, con, dphi;
- int i;
- eccnth = .5 * e;
- Phi = geometry::math::half_pi<T>() - 2. * atan (ts);
- i = N_ITER;
- do {
- con = e * sin (Phi);
- dphi = geometry::math::half_pi<T>() - 2. * atan (ts * math::pow((T(1) - con) /
- (T(1) + con), eccnth)) - Phi;
- Phi += dphi;
- } while ( geometry::math::abs(dphi) > TOL && --i);
- if (i <= 0)
- BOOST_THROW_EXCEPTION( projection_exception(error_non_con_inv_phi2) );
- return Phi;
- }
- } // namespace detail
- }}} // namespace boost::geometry::projections
- #endif
|