Add patch for failing prerm scripts

[ceph.git] / ceph / src / boost / boost / qvm / detail / determinant_impl.hpp

//Copyright (c) 2008-2016 Emil Dotchevski and Reverge Studios, Inc.\r
\r
//Distributed under the Boost Software License, Version 1.0. (See accompanying\r
//file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)\r
\r
#ifndef UUID_3DCF6B90AE0E11DE9A315BE555D89593\r
#define UUID_3DCF6B90AE0E11DE9A315BE555D89593\r
\r
#include <boost/qvm/inline.hpp>\r
#include <boost/qvm/mat_traits_array.hpp>\r
#include <boost/qvm/static_assert.hpp>\r
\r
namespace\r
boost\r
    {\r
    namespace\r
    qvm\r
        {\r
        namespace\r
        qvm_detail\r
            {\r
            template <int N>\r
            struct\r
            det_size\r
                {\r
                };\r
\r
            template <class M>\r
            BOOST_QVM_INLINE_TRIVIAL\r
            typename mat_traits<M>::scalar_type\r
            determinant_impl_( M const & a, det_size<2> )\r
                {\r
                return\r
                    mat_traits<M>::template read_element<0,0>(a) * mat_traits<M>::template read_element<1,1>(a) -\r
                    mat_traits<M>::template read_element<1,0>(a) * mat_traits<M>::template read_element<0,1>(a);\r
                }\r
\r
            template <class M,int N>\r
            BOOST_QVM_INLINE_RECURSION\r
            typename mat_traits<M>::scalar_type\r
            determinant_impl_( M const & a, det_size<N> )\r
                {\r
                typedef typename mat_traits<M>::scalar_type T;\r
                T m[N-1][N-1];\r
                T det=T(0);\r
                for( int j1=0; j1!=N; ++j1 )\r
                    {\r
                    for( int i=1; i!=N; ++i )\r
                        {\r
                        int j2 = 0;\r
                        for( int j=0; j!=N; ++j )\r
                            {\r
                            if( j==j1 )\r
                                continue;\r
                            m[i-1][j2] = mat_traits<M>::read_element_idx(i,j,a);\r
                            ++j2;\r
                            }\r
                        }\r
                    T d=determinant_impl_(m,det_size<N-1>());\r
                    if( j1&1 )\r
                        d=-d;\r
                    det += mat_traits<M>::read_element_idx(0,j1,a) * d;\r
                    }\r
                return det;\r
                }\r
\r
            template <class M>\r
            BOOST_QVM_INLINE_TRIVIAL\r
            typename mat_traits<M>::scalar_type\r
            determinant_impl( M const & a )\r
                {\r
                BOOST_QVM_STATIC_ASSERT(mat_traits<M>::rows==mat_traits<M>::cols);\r
                return determinant_impl_(a,det_size<mat_traits<M>::rows>());\r
                }\r
            }\r
        }\r
    }\r
\r
#endif\r