summaryrefslogtreecommitdiff log msg author committer range
blob: e64a46e4a8a36987d9ba1e71f6ffe15864d752e9 (plain)
 ```1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 ``` ``````// Boost.Geometry (aka GGL, Generic Geometry Library) // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands. // Copyright (c) 2008-2012 Bruno Lalande, Paris, France. // Copyright (c) 2009-2012 Mateusz Loskot, London, UK. // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands. // 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) #ifndef BOOST_GEOMETRY_STRATEGIES_TRANSFORM_INVERSE_TRANSFORMER_HPP #define BOOST_GEOMETRY_STRATEGIES_TRANSFORM_INVERSE_TRANSFORMER_HPP // Remove the ublas checking, otherwise the inverse might fail // (while nothing seems to be wrong) #ifdef BOOST_UBLAS_TYPE_CHECK #undef BOOST_UBLAS_TYPE_CHECK #endif #define BOOST_UBLAS_TYPE_CHECK 0 #include #include #include namespace boost { namespace geometry { namespace strategy { namespace transform { /*! \brief Transformation strategy to do an inverse transformation in a Cartesian coordinate system \ingroup strategies */ template < typename CalculationType, std::size_t Dimension1, std::size_t Dimension2 > class inverse_transformer : public ublas_transformer { public : template inline inverse_transformer(Transformer const& input) { typedef boost::numeric::ublas::matrix matrix_type; // create a working copy of the input matrix_type copy(input.matrix()); // create a permutation matrix for the LU-factorization typedef boost::numeric::ublas::permutation_matrix<> permutation_matrix; permutation_matrix pm(copy.size1()); // perform LU-factorization int res = boost::numeric::ublas::lu_factorize(copy, pm); if( res == 0 ) { // create identity matrix this->m_matrix.assign(boost::numeric::ublas::identity_matrix(copy.size1())); // backsubstitute to get the inverse boost::numeric::ublas::lu_substitute(copy, pm, this->m_matrix); } } }; }} // namespace strategy::transform }} // namespace boost::geometry #endif // BOOST_GEOMETRY_STRATEGIES_TRANSFORM_INVERSE_TRANSFORMER_HPP ``````