// // Copyright (c) 2002-2003 // Toon Knapen, Kresimir Fresl, Joerg Walter // // Distributed under 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_UBLAS_RAW_ #define _BOOST_UBLAS_RAW_ namespace boost { namespace numeric { namespace ublas { namespace raw { // We need data_const() mostly due to MSVC 6.0. // But how shall we write portable code otherwise? template < typename V > BOOST_UBLAS_INLINE int size( const V &v ) ; template < typename V > BOOST_UBLAS_INLINE int size( const vector_reference &v ) ; template < typename M > BOOST_UBLAS_INLINE int size1( const M &m ) ; template < typename M > BOOST_UBLAS_INLINE int size2( const M &m ) ; template < typename M > BOOST_UBLAS_INLINE int size1( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE int size2( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m, row_major_tag ) ; template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m, column_major_tag ) ; template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m ) ; template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const matrix_reference &m ) ; template < typename V > BOOST_UBLAS_INLINE int stride( const V &v ) ; template < typename V > BOOST_UBLAS_INLINE int stride( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE int stride( const vector_slice &v ) ; template < typename M > BOOST_UBLAS_INLINE int stride( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE int stride( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE int stride1( const M &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride2( const M &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_reference &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE int stride1( const c_matrix &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE int stride2( const c_matrix &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_slice &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_slice &m ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::const_pointer data( const MV &mv ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::const_pointer data_const( const MV &mv ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::pointer data( MV &mv ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_reference &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data_const( const vector_reference &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_reference &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer data( const c_vector &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer data_const( const c_vector &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::pointer data( c_vector &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_slice &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data_const( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data_const( const vector_slice &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_slice &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data_const( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_reference &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::array_type::const_pointer data( const c_matrix &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::array_type::const_pointer data_const( const c_matrix &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::pointer data( c_matrix &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data_const( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data_const( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data( const matrix_slice &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data_const( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer data_const( const matrix_slice &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_slice &m ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::const_pointer base( const MV &mv ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::const_pointer base_const( const MV &mv ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::pointer base( MV &mv ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base( const vector_reference &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base_const( const vector_reference &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_reference &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer base( const c_vector &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer base_const( const c_vector &v ) ; template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::pointer base( c_vector &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base( const vector_slice &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base_const( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer base_const( const vector_slice &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_slice &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base_const( const matrix_reference &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_reference &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::array_type::const_pointer base( const c_matrix &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::array_type::const_pointer base_const( const c_matrix &m ) ; template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::pointer base( c_matrix &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base_const( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base_const( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base( const matrix_slice &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base_const( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::array_type::const_pointer base_const( const matrix_slice &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_slice &m ) ; template < typename MV > BOOST_UBLAS_INLINE typename MV::size_type start( const MV &mv ) ; template < typename V > BOOST_UBLAS_INLINE typename V::size_type start( const vector_range &v ) ; template < typename V > BOOST_UBLAS_INLINE typename V::size_type start( const vector_slice &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_row &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_column &v ) ; template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_range &m ) ; template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_slice &m ) ; template < typename V > BOOST_UBLAS_INLINE int size( const V &v ) { return v.size() ; } template < typename V > BOOST_UBLAS_INLINE int size( const vector_reference &v ) { return size( v ) ; } template < typename M > BOOST_UBLAS_INLINE int size1( const M &m ) { return m.size1() ; } template < typename M > BOOST_UBLAS_INLINE int size2( const M &m ) { return m.size2() ; } template < typename M > BOOST_UBLAS_INLINE int size1( const matrix_reference &m ) { return size1( m.expression() ) ; } template < typename M > BOOST_UBLAS_INLINE int size2( const matrix_reference &m ) { return size2( m.expression() ) ; } template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m, row_major_tag ) { return m.size2() ; } template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m, column_major_tag ) { return m.size1() ; } template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const M &m ) { return leading_dimension( m, typename M::orientation_category() ) ; } template < typename M > BOOST_UBLAS_INLINE int leading_dimension( const matrix_reference &m ) { return leading_dimension( m.expression() ) ; } template < typename V > BOOST_UBLAS_INLINE int stride( const V &v ) { return 1 ; } template < typename V > BOOST_UBLAS_INLINE int stride( const vector_range &v ) { return stride( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE int stride( const vector_slice &v ) { return v.stride() * stride( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride( const matrix_row &v ) { return stride2( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride( const matrix_column &v ) { return stride1( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride1( const M &m ) { typedef typename M::functor_type functor_type; return functor_type::one1( m.size1(), m.size2() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride2( const M &m ) { typedef typename M::functor_type functor_type; return functor_type::one2( m.size1(), m.size2() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_reference &m ) { return stride1( m.expression() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_reference &m ) { return stride2( m.expression() ) ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE int stride1( const c_matrix &m ) { return N ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE int stride2( const c_matrix &m ) { return 1 ; } template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_range &m ) { return stride1( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride1( const matrix_slice &m ) { return m.stride1() * stride1( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_range &m ) { return stride2( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE int stride2( const matrix_slice &m ) { return m.stride2() * stride2( m.data() ) ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::array_type::const_pointer data( const MV &mv ) { return &mv.data().begin()[0] ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::array_type::const_pointer data_const( const MV &mv ) { return &mv.data().begin()[0] ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::pointer data( MV &mv ) { return &mv.data().begin()[0] ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_reference &v ) { return data( v.expression () ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data_const( const vector_reference &v ) { return data_const( v.expression () ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_reference &v ) { return data( v.expression () ) ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer data( const c_vector &v ) { return v.data() ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::array_type::const_pointer data_const( const c_vector &v ) { return v.data() ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::pointer data( c_vector &v ) { return v.data() ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_range &v ) { return data( v.data() ) + v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data( const vector_slice &v ) { return data( v.data() ) + v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::array_type::const_pointer data_const( const vector_range &v ) { return data_const( v.data() ) + v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer data_const( const vector_slice &v ) { return data_const( v.data() ) + v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_range &v ) { return data( v.data() ) + v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer data( vector_slice &v ) { return data( v.data() ) + v.start() * stride (v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data( const matrix_reference &m ) { return data( m.expression () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data_const( const matrix_reference &m ) { return data_const( m.expression () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_reference &m ) { return data( m.expression () ) ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::const_pointer data( const c_matrix &m ) { return m.data() ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::const_pointer data_const( const c_matrix &m ) { return m.data() ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::pointer data( c_matrix &m ) { return m.data() ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data( const matrix_row &v ) { return data( v.data() ) + v.index() * stride1( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data( const matrix_column &v ) { return data( v.data() ) + v.index() * stride2( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data_const( const matrix_row &v ) { return data_const( v.data() ) + v.index() * stride1( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data_const( const matrix_column &v ) { return data_const( v.data() ) + v.index() * stride2( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_row &v ) { return data( v.data() ) + v.index() * stride1( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_column &v ) { return data( v.data() ) + v.index() * stride2( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data( const matrix_range &m ) { return data( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data( const matrix_slice &m ) { return data( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data_const( const matrix_range &m ) { return data_const( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer data_const( const matrix_slice &m ) { return data_const( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_range &m ) { return data( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer data( matrix_slice &m ) { return data( m.data() ) + m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::const_pointer base( const MV &mv ) { return &mv.data().begin()[0] ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::const_pointer base_const( const MV &mv ) { return &mv.data().begin()[0] ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::array_type::pointer base( MV &mv ) { return &mv.data().begin()[0] ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base( const vector_reference &v ) { return base( v.expression () ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base_const( const vector_reference &v ) { return base_const( v.expression () ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_reference &v ) { return base( v.expression () ) ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::const_pointer base( const c_vector &v ) { return v.data() ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::array_type::const_pointer base_const( const c_vector &v ) { return v.data() ; } template < typename T, std::size_t N > BOOST_UBLAS_INLINE typename c_vector::pointer base( c_vector &v ) { return v.data() ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base( const vector_range &v ) { return base( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base( const vector_slice &v ) { return base( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base_const( const vector_range &v ) { return base_const( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::const_pointer base_const( const vector_slice &v ) { return base_const( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_range &v ) { return base( v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::array_type::pointer base( vector_slice &v ) { return base( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base( const matrix_reference &m ) { return base( m.expression () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base_const( const matrix_reference &m ) { return base_const( m.expression () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_reference &m ) { return base( m.expression () ) ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::const_pointer base( const c_matrix &m ) { return m.data() ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::array_type::const_pointer base_const( const c_matrix &m ) { return m.data() ; } template < typename T, std::size_t M, std::size_t N > BOOST_UBLAS_INLINE typename c_matrix::pointer base( c_matrix &m ) { return m.data() ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base( const matrix_row &v ) { return base( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base( const matrix_column &v ) { return base( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base_const( const matrix_row &v ) { return base_const( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base_const( const matrix_column &v ) { return base_const( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_row &v ) { return base( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_column &v ) { return base( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base( const matrix_range &m ) { return base( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base( const matrix_slice &m ) { return base( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base_const( const matrix_range &m ) { return base_const( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::const_pointer base_const( const matrix_slice &m ) { return base_const( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_range &m ) { return base( m.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::array_type::pointer base( matrix_slice &m ) { return base( m.data() ) ; } template < typename MV > BOOST_UBLAS_INLINE typename MV::size_type start( const MV &mv ) { return 0 ; } template < typename V > BOOST_UBLAS_INLINE typename V::size_type start( const vector_range &v ) { return v.start() * stride (v.data() ) ; } template < typename V > BOOST_UBLAS_INLINE typename V::size_type start( const vector_slice &v ) { return v.start() * stride (v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_row &v ) { return v.index() * stride1( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_column &v ) { return v.index() * stride2( v.data() ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_range &m ) { return m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } template < typename M > BOOST_UBLAS_INLINE typename M::size_type start( const matrix_slice &m ) { return m.start1() * stride1( m.data () ) + m.start2() * stride2( m.data () ) ; } }}}} #endif