Chris@49: // Copyright (C) 2008-2013 NICTA (www.nicta.com.au)
Chris@49: // Copyright (C) 2008-2013 Conrad Sanderson
Chris@49: // 
Chris@49: // This Source Code Form is subject to the terms of the Mozilla Public
Chris@49: // License, v. 2.0. If a copy of the MPL was not distributed with this
Chris@49: // file, You can obtain one at http://mozilla.org/MPL/2.0/.
Chris@49: 
Chris@49: 
Chris@49: //! \addtogroup unwrap
Chris@49: //! @{
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_default(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_fixed
Chris@49:   {
Chris@49:   typedef T1 stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   unwrap_fixed(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct unwrap_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_redirect<T1, false> { typedef unwrap_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_redirect<T1, true>  { typedef unwrap_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap : public unwrap_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   inline
Chris@49:   unwrap(const T1& A)
Chris@49:     : unwrap_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap< Mat<eT> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap(const Mat<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap< Row<eT> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap(const Row<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap< Col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49: 
Chris@49:   inline
Chris@49:   unwrap(const Col<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename T2, typename glue_type>
Chris@49: struct unwrap< mtGlue<out_eT, T1, T2, glue_type> >
Chris@49:   {
Chris@49:   typedef Mat<out_eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap(const mtGlue<out_eT, T1, T2, glue_type>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<out_eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename op_type>
Chris@49: struct unwrap< mtOp<out_eT, T1, op_type> >
Chris@49:   {
Chris@49:   typedef Mat<out_eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap(const mtOp<out_eT, T1, op_type>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<out_eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //
Chris@49: //
Chris@49: //
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct quasi_unwrap_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap_default(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   // NOTE: DO NOT DIRECTLY CHECK FOR ALIASING BY TAKING THE ADDRESS OF THE "M" OBJECT IN ANY quasi_unwrap CLASS !!!
Chris@49:   const Mat<eT> M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct quasi_unwrap_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   quasi_unwrap_fixed(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&M) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct quasi_unwrap_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct quasi_unwrap_redirect<T1, false> { typedef quasi_unwrap_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct quasi_unwrap_redirect<T1, true>  { typedef quasi_unwrap_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct quasi_unwrap : public quasi_unwrap_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename quasi_unwrap_redirect<T1, is_Mat_fixed<T1>::value >::result quasi_unwrap_extra;
Chris@49:   
Chris@49:   static const bool has_subview = quasi_unwrap_extra::has_subview;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const T1& A)
Chris@49:     : quasi_unwrap_extra(A)
Chris@49:     {
Chris@49:     }
Chris@49:   
Chris@49:   using quasi_unwrap_extra::M;
Chris@49:   using quasi_unwrap_extra::is_alias;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct quasi_unwrap< Mat<eT> >
Chris@49:   {
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const Mat<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&M) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct quasi_unwrap< Row<eT> >
Chris@49:   {
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const Row<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Row<eT>& M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&M) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct quasi_unwrap< Col<eT> >
Chris@49:   {
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const Col<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Col<eT>& M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&M) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct quasi_unwrap< subview_col<eT> >
Chris@49:   {
Chris@49:   static const bool has_subview = true;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const subview_col<eT>& A)
Chris@49:     : M  ( const_cast<eT*>( A.colptr(0) ), A.n_rows, 1, false, false )
Chris@49:     , src( A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT>  M;
Chris@49:   const Mat<eT>& src;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&src) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename T2, typename glue_type>
Chris@49: struct quasi_unwrap< mtGlue<out_eT, T1, T2, glue_type> >
Chris@49:   {
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const mtGlue<out_eT, T1, T2, glue_type>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<out_eT> M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename op_type>
Chris@49: struct quasi_unwrap< mtOp<out_eT, T1, op_type> >
Chris@49:   {
Chris@49:   static const bool has_subview = false;
Chris@49:   
Chris@49:   inline
Chris@49:   quasi_unwrap(const mtOp<out_eT, T1, op_type>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<out_eT> M;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //
Chris@49: //
Chris@49: //
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check_default(const T1& A, const Mat<eT>&)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check_default(const T1& A, const bool)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check_fixed(const T1& A, const Mat<eT>& B)
Chris@49:     : M_local( (&A == &B) ? new T1(A) : 0 )
Chris@49:     , M      ( (&A == &B) ? *M_local  : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check_fixed(const T1& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new T1(A) : 0 )
Chris@49:     , M      ( is_alias ? *M_local  : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct unwrap_check_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check_redirect<T1, false> { typedef unwrap_check_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check_redirect<T1, true>  { typedef unwrap_check_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check : public unwrap_check_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   inline unwrap_check(const T1& A, const Mat<typename T1::elem_type>& B)
Chris@49:     : unwrap_check_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   
Chris@49:   inline unwrap_check(const T1& A, const bool is_alias)
Chris@49:     : unwrap_check_redirect< T1, is_Mat_fixed<T1>::value >::result(A, is_alias)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap_check< Mat<eT> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check(const Mat<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&A == &B) ? new Mat<eT>(A) : 0 )
Chris@49:     , M      ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   unwrap_check(const Mat<eT>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Mat<eT>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap_check< Row<eT> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check(const Row<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&A == &B) ? new Row<eT>(A) : 0 )
Chris@49:     , M      ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   unwrap_check(const Row<eT>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Row<eT>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct unwrap_check< Col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   unwrap_check(const Col<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&A == &B) ? new Col<eT>(A) : 0 )
Chris@49:     , M      ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   unwrap_check(const Col<eT>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Col<eT>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //
Chris@49: //
Chris@49: //
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct unwrap_check_mixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT1;
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const T1& A, const Mat<eT2>&)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const T1& A, const bool)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   const Mat<eT1> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT1>
Chris@49: struct unwrap_check_mixed< Mat<eT1> >
Chris@49:   {
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Mat<eT1>& A, const Mat<eT2>& B)
Chris@49:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Mat<eT1>(A) : 0 )
Chris@49:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Mat<eT1>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Mat<eT1>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check_mixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Mat<eT1>* M_local;
Chris@49:   const Mat<eT1>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT1>
Chris@49: struct unwrap_check_mixed< Row<eT1> >
Chris@49:   {
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Row<eT1>& A, const Mat<eT2>& B)
Chris@49:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Row<eT1>(A) : 0 )
Chris@49:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Row<eT1>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Row<eT1>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check_mixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Row<eT1>* M_local;
Chris@49:   const Row<eT1>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT1>
Chris@49: struct unwrap_check_mixed< Col<eT1> >
Chris@49:   {
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Col<eT1>& A, const Mat<eT2>& B)
Chris@49:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Col<eT1>(A) : 0 )
Chris@49:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   inline
Chris@49:   unwrap_check_mixed(const Col<eT1>& A, const bool is_alias)
Chris@49:     : M_local( is_alias ? new Col<eT1>(A) : 0 )
Chris@49:     , M      ( is_alias ? (*M_local)      : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~unwrap_check_mixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Col<eT1>* M_local;
Chris@49:   const Col<eT1>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //
Chris@49: //
Chris@49: //
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_default(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_fixed(const T1& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_redirect<T1, false> { typedef partial_unwrap_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_redirect<T1, true>  { typedef partial_unwrap_fixed<T1>   result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap : public partial_unwrap_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   inline
Chris@49:   partial_unwrap(const T1& A)
Chris@49:     : partial_unwrap_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Mat<eT> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Mat<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Row<eT> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Row<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Col<eT>& A)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap class is only used by as_scalar(),
Chris@49: // NOTE: which doesn't need to check for aliasing (ie. it doesn't need to compare the address of M with another matrix)
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< subview_col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const subview_col<eT>& A)
Chris@49:     : M  ( const_cast<eT*>( A.colptr(0) ), A.n_rows, false, false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_htrans_default(const Op<T1, op_htrans>& A)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_htrans_fixed(const Op<T1, op_htrans>& A)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_htrans_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans_redirect<T1, false> { typedef partial_unwrap_htrans_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans_redirect<T1, true>  { typedef partial_unwrap_htrans_fixed<T1>   result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap< Op<T1, op_htrans> > : public partial_unwrap_htrans_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   inline partial_unwrap(const Op<T1, op_htrans>& A)
Chris@49:     : partial_unwrap_htrans_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Mat<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Mat<eT>, op_htrans>& A)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Row<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Row<eT>, op_htrans>& A)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Col<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Col<eT>, op_htrans>& A)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap class is only used by as_scalar(),
Chris@49: // NOTE: which doesn't need to check for aliasing (ie. it doesn't need to compare the address of M with another matrix)
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< subview_col<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< subview_col<eT>, op_htrans>& A)
Chris@49:     : M  ( const_cast<eT*>( A.m.colptr(0) ), A.m.n_rows, false, false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans2_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_htrans2_default(const Op<T1, op_htrans2>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans2_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_htrans2_fixed(const Op<T1, op_htrans2>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT  val;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_htrans2_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans2_redirect<T1, false> { typedef partial_unwrap_htrans2_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_htrans2_redirect<T1, true>  { typedef partial_unwrap_htrans2_fixed<T1>   result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap< Op<T1, op_htrans2> > : public partial_unwrap_htrans2_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   inline partial_unwrap(const Op<T1, op_htrans2>& A)
Chris@49:     : partial_unwrap_htrans2_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Mat<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Mat<eT>, op_htrans2>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Row<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Row<eT>, op_htrans2>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< Col<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< Col<eT>, op_htrans2>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap class is only used by as_scalar(),
Chris@49: // NOTE: which doesn't need to check for aliasing (ie. it doesn't need to compare the address of M with another matrix)
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< Op< subview_col<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const Op< subview_col<eT>, op_htrans2>& A)
Chris@49:     : val( A.aux )
Chris@49:     , M  ( const_cast<eT*>( A.m.colptr(0) ), A.m.n_rows, false, false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_scalar_times_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_scalar_times_default(const eOp<T1, eop_scalar_times>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_scalar_times_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_scalar_times_fixed(const eOp<T1, eop_scalar_times>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT  val;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_scalar_times_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_scalar_times_redirect<T1, false> { typedef partial_unwrap_scalar_times_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_scalar_times_redirect<T1, true>  { typedef partial_unwrap_scalar_times_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap< eOp<T1, eop_scalar_times> > : public partial_unwrap_scalar_times_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<T1, eop_scalar_times>& A)
Chris@49:     : partial_unwrap_scalar_times_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Mat<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Mat<eT>,eop_scalar_times>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Row<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Row<eT>,eop_scalar_times>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Col<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Col<eT>,eop_scalar_times>& A)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // TODO: struct partial_unwrap< eOp<subview_col<eT>, eop_scalar_times> >
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_neg_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_neg_default(const eOp<T1, eop_neg>& A)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_neg_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_neg_fixed(const eOp<T1, eop_neg>& A)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_neg_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_neg_redirect<T1, false> { typedef partial_unwrap_neg_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_neg_redirect<T1, true>  { typedef partial_unwrap_neg_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap< eOp<T1, eop_neg> > : public partial_unwrap_neg_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<T1, eop_neg>& A)
Chris@49:     : partial_unwrap_neg_redirect< T1, is_Mat_fixed<T1>::value >::result(A)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Mat<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Mat<eT>,eop_neg>& A)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Row<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Row<eT>,eop_neg>& A)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap< eOp<Col<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap(const eOp<Col<eT>,eop_neg>& A)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // TODO: struct partial_unwrap< eOp<subview_col<eT>, eop_neg> >
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_check_default(const T1& A, const Mat<eT>&)
Chris@49:     : M(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_check_fixed(const T1& A, const Mat<eT>& B)
Chris@49:     : M_local( (&A == &B) ? new T1(A)  : 0 )
Chris@49:     , M      ( (&A == &B) ? (*M_local) : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_check_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_redirect<T1, false> { typedef partial_unwrap_check_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_redirect<T1, true>  { typedef partial_unwrap_check_fixed<T1>   result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check : public partial_unwrap_check_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline partial_unwrap_check(const T1& A, const Mat<eT>& B)
Chris@49:     : partial_unwrap_check_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Mat<eT> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49: 
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Mat<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A == &B) ? new Mat<eT>(A) : 0 )
Chris@49:     , M       ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Row<eT> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Row<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A == &B) ? new Row<eT>(A) : 0 )
Chris@49:     , M       ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Col<eT>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A == &B) ? new Col<eT>(A) : 0 )
Chris@49:     , M       ( (&A == &B) ? (*M_local)     : A )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap_check class is only used by the glue_times class,
Chris@49: // NOTE: which relies on partial_unwrap_check to check for aliasing
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< subview_col<eT> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const subview_col<eT>& A, const Mat<eT>& B)
Chris@49:     : M  ( const_cast<eT*>( A.colptr(0) ), A.n_rows, (&(A.m) == &B), false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_check_htrans_default(const Op<T1, op_htrans>& A, const Mat<eT>&)
Chris@49:     : M(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_check_htrans_fixed(const Op<T1, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&(A.m) == &B) ? new T1(A.m) : 0   )
Chris@49:     , M      ( (&(A.m) == &B) ? (*M_local)  : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check_htrans_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_check_htrans_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans_redirect<T1, false> { typedef partial_unwrap_check_htrans_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans_redirect<T1, true>  { typedef partial_unwrap_check_htrans_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check< Op<T1, op_htrans> > : public partial_unwrap_check_htrans_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline partial_unwrap_check(const Op<T1, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : partial_unwrap_check_htrans_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Mat<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Mat<eT>, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Row<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Row<eT>, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A.m == &B) ? new Row<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Col<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Col<eT>, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : M_local ( (&A.m == &B) ? new Col<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap_check class is only used by the glue_times class,
Chris@49: // NOTE: which relies on partial_unwrap_check to check for aliasing
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< subview_col<eT>, op_htrans> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< subview_col<eT>, op_htrans>& A, const Mat<eT>& B)
Chris@49:     : M  ( const_cast<eT*>( A.m.colptr(0) ), A.m.n_rows, (&(A.m.m) == &B), false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(1); }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = false;
Chris@49:   
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans2_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_check_htrans2_default(const Op<T1, op_htrans2>& A, const Mat<eT>&)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans2_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_check_htrans2_fixed(const Op<T1, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : val    (A.aux)
Chris@49:     , M_local( (&(A.m) == &B) ? new T1(A.m) : 0   )
Chris@49:     , M      ( (&(A.m) == &B) ? (*M_local)  : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check_htrans2_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT  val;
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_check_htrans2_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans2_redirect<T1, false> { typedef partial_unwrap_check_htrans2_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_htrans2_redirect<T1, true>  { typedef partial_unwrap_check_htrans2_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check< Op<T1, op_htrans2> > : public partial_unwrap_check_htrans2_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline partial_unwrap_check(const Op<T1, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : partial_unwrap_check_htrans2_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Mat<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Mat<eT>, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : val     (A.aux)
Chris@49:     , M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const eT       val;
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Row<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Row<eT>, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : val     (A.aux)
Chris@49:     , M_local ( (&A.m == &B) ? new Row<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const eT       val;
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< Col<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< Col<eT>, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : val     (A.aux)
Chris@49:     , M_local ( (&A.m == &B) ? new Col<eT>(A.m) : 0   )
Chris@49:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   // the order below is important
Chris@49:   const eT       val;
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap_check class is only used by the glue_times class,
Chris@49: // NOTE: which relies on partial_unwrap_check to check for aliasing
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< Op< subview_col<eT>, op_htrans2> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const Op< subview_col<eT>, op_htrans2>& A, const Mat<eT>& B)
Chris@49:     : val( A.aux )
Chris@49:     , M  ( const_cast<eT*>( A.m.colptr(0) ), A.m.n_rows, (&(A.m.m) == &B), false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = true;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_scalar_times_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_check_scalar_times_default(const eOp<T1, eop_scalar_times>& A, const Mat<eT>&)
Chris@49:     : val(A.aux)
Chris@49:     , M  (A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_scalar_times_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_check_scalar_times_fixed(const eOp<T1, eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : val    ( A.aux )
Chris@49:     , M_local( (&(A.P.Q) == &B) ? new T1(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? (*M_local)    : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check_scalar_times_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT  val;
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_check_scalar_times_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_scalar_times_redirect<T1, false> { typedef partial_unwrap_check_scalar_times_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_scalar_times_redirect<T1, true>  { typedef partial_unwrap_check_scalar_times_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check< eOp<T1, eop_scalar_times> > : public partial_unwrap_check_scalar_times_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline partial_unwrap_check(const eOp<T1, eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : partial_unwrap_check_scalar_times_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Mat<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Mat<eT>,eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : val    (A.aux)
Chris@49:     , M_local( (&(A.P.Q) == &B) ? new Mat<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Row<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Row<eT>,eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : val(A.aux)
Chris@49:     , M_local( (&(A.P.Q) == &B) ? new Row<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Col<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Col<eT>,eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : val    ( A.aux )
Chris@49:     , M_local( (&(A.P.Q) == &B) ? new Col<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT       val;
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap_check class is only used by the glue_times class,
Chris@49: // NOTE: which relies on partial_unwrap_check to check for aliasing
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<subview_col<eT>, eop_scalar_times> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<subview_col<eT>,eop_scalar_times>& A, const Mat<eT>& B)
Chris@49:     : val( A.aux )
Chris@49:     , M  ( const_cast<eT*>( A.P.Q.colptr(0) ), A.P.Q.n_rows, (&(A.P.Q.m) == &B), false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return val; }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const eT      val;
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_neg_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef Mat<eT>                stored_type;
Chris@49:   
Chris@49:   inline
Chris@49:   partial_unwrap_check_neg_default(const eOp<T1, eop_neg>& A, const Mat<eT>&)
Chris@49:     : M(A.P.Q)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Mat<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_neg_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   typedef T1                     stored_type;
Chris@49:   
Chris@49:   inline explicit
Chris@49:   partial_unwrap_check_neg_fixed(const eOp<T1, eop_neg>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&(A.P.Q) == &B) ? new T1(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? (*M_local)    : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check_neg_fixed()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const T1* M_local;
Chris@49:   const T1& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct partial_unwrap_check_neg_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_neg_redirect<T1, false> { typedef partial_unwrap_check_neg_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check_neg_redirect<T1, true>  { typedef partial_unwrap_check_neg_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct partial_unwrap_check< eOp<T1, eop_neg> > : public partial_unwrap_check_neg_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   inline partial_unwrap_check(const eOp<T1, eop_neg>& A, const Mat<eT>& B)
Chris@49:     : partial_unwrap_check_neg_redirect< T1, is_Mat_fixed<T1>::value >::result(A, B)
Chris@49:     {
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Mat<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Mat<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Mat<eT>,eop_neg>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&(A.P.Q) == &B) ? new Mat<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Mat<eT>* M_local;
Chris@49:   const Mat<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Row<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Row<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Row<eT>,eop_neg>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&(A.P.Q) == &B) ? new Row<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Row<eT>* M_local;
Chris@49:   const Row<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<Col<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<Col<eT>,eop_neg>& A, const Mat<eT>& B)
Chris@49:     : M_local( (&(A.P.Q) == &B) ? new Col<eT>(A.P.Q) : 0     )
Chris@49:     , M      ( (&(A.P.Q) == &B) ? *M_local           : A.P.Q )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   inline
Chris@49:   ~partial_unwrap_check()
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     
Chris@49:     if(M_local) { delete M_local; }
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Col<eT>* M_local;
Chris@49:   const Col<eT>& M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: // NOTE: we can get away with this shortcut as the partial_unwrap_check class is only used by the glue_times class,
Chris@49: // NOTE: which relies on partial_unwrap_check to check for aliasing
Chris@49: template<typename eT>
Chris@49: struct partial_unwrap_check< eOp<subview_col<eT>, eop_neg> >
Chris@49:   {
Chris@49:   typedef Col<eT> stored_type;
Chris@49:   
Chris@49:   arma_hot inline
Chris@49:   partial_unwrap_check(const eOp<subview_col<eT>,eop_neg>& A, const Mat<eT>& B)
Chris@49:     : M  ( const_cast<eT*>( A.P.Q.colptr(0) ), A.P.Q.n_rows, (&(A.P.Q.m) == &B), false )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_hot arma_inline eT get_val() const { return eT(-1); }
Chris@49:   
Chris@49:   static const bool do_trans = false;
Chris@49:   static const bool do_times = true;
Chris@49:   
Chris@49:   const Col<eT> M;
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //! @}