max@0: // Copyright (C) 2008-2011 NICTA (www.nicta.com.au)
max@0: // Copyright (C) 2008-2011 Conrad Sanderson
max@0: // 
max@0: // This file is part of the Armadillo C++ library.
max@0: // It is provided without any warranty of fitness
max@0: // for any purpose. You can redistribute this file
max@0: // and/or modify it under the terms of the GNU
max@0: // Lesser General Public License (LGPL) as published
max@0: // by the Free Software Foundation, either version 3
max@0: // of the License or (at your option) any later version.
max@0: // (see http://www.opensource.org/licenses for more info)
max@0: 
max@0: 
max@0: //! \addtogroup unwrap
max@0: //! @{
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class unwrap
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline unwrap(const T1& A)   // TODO: change this to Base ?
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap< Mat<eT> >
max@0:   {
max@0:   public:
max@0: 
max@0:   inline unwrap(const Mat<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0: 
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap< Row<eT> >
max@0:   {
max@0:   public:
max@0: 
max@0:   inline unwrap(const Row<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0: 
max@0:   const Row<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap< Col<eT> >
max@0:   {
max@0:   public:
max@0: 
max@0:   inline unwrap(const Col<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0: 
max@0:   const Col<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename out_eT, typename T1, typename T2, typename glue_type>
max@0: class unwrap< mtGlue<out_eT, T1, T2, glue_type> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline unwrap(const mtGlue<out_eT, T1, T2, glue_type>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   const Mat<out_eT> M;
max@0:   };
max@0: 
max@0: 
max@0: template<typename out_eT, typename T1, typename op_type>
max@0: class unwrap< mtOp<out_eT, T1, op_type> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline unwrap(const mtOp<out_eT, T1, op_type>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   const Mat<out_eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: //
max@0: //
max@0: //
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class unwrap_check
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   unwrap_check(const T1& A, const Mat<eT>& B)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap_check< Mat<eT> >
max@0:   {
max@0:   public:
max@0: 
max@0:   inline
max@0:   unwrap_check(const Mat<eT>& A, const Mat<eT>& B)
max@0:     : M_local( (&A == &B) ? new Mat<eT>(A) : 0 )
max@0:     , M      ( (&A == &B) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap_check< Row<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   unwrap_check(const Row<eT>& A, const Mat<eT>& B)
max@0:     : M_local( (&A == reinterpret_cast<const Row<eT>*>(&B)) ? new Row<eT>(A) : 0 )
max@0:     , M      ( (&A == reinterpret_cast<const Row<eT>*>(&B)) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Row<eT>* M_local;
max@0:   const Row<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class unwrap_check< Col<eT> >
max@0:   {
max@0:   public:
max@0: 
max@0:   inline
max@0:   unwrap_check(const Col<eT>& A, const Mat<eT>& B)
max@0:     : M_local( (&A == reinterpret_cast<const Col<eT>*>(&B)) ? new Col<eT>(A) : 0 )
max@0:     , M      ( (&A == reinterpret_cast<const Col<eT>*>(&B)) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Col<eT>* M_local;
max@0:   const Col<eT>& M;
max@0:   
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: //
max@0: //
max@0: //
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class unwrap_check_mixed
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT1;
max@0:   
max@0:   template<typename eT2>
max@0:   inline
max@0:   unwrap_check_mixed(const T1& A, const Mat<eT2>& B)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~unwrap_check_mixed()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   const Mat<eT1> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT1>
max@0: class unwrap_check_mixed< Mat<eT1> >
max@0:   {
max@0:   public:
max@0:   
max@0:   template<typename eT2>
max@0:   inline
max@0:   unwrap_check_mixed(const Mat<eT1>& A, const Mat<eT2>& B)
max@0:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Mat<eT1>(A) : 0 )
max@0:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check_mixed()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT1>* M_local;
max@0:   const Mat<eT1>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT1>
max@0: class unwrap_check_mixed< Row<eT1> >
max@0:   {
max@0:   public:
max@0:   
max@0:   template<typename eT2>
max@0:   inline
max@0:   unwrap_check_mixed(const Row<eT1>& A, const Mat<eT2>& B)
max@0:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Row<eT1>(A) : 0 )
max@0:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check_mixed()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Row<eT1>* M_local;
max@0:   const Row<eT1>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT1>
max@0: class unwrap_check_mixed< Col<eT1> >
max@0:   {
max@0:   public:
max@0:   
max@0:   template<typename eT2>
max@0:   inline
max@0:   unwrap_check_mixed(const Col<eT1>& A, const Mat<eT2>& B)
max@0:     : M_local( (void_ptr(&A) == void_ptr(&B)) ? new Col<eT1>(A) : 0 )
max@0:     , M      ( (void_ptr(&A) == void_ptr(&B)) ? (*M_local)      : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~unwrap_check_mixed()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   // the order below is important
max@0:   const Col<eT1>* M_local;
max@0:   const Col<eT1>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: //
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline partial_unwrap(const T1& A)  // TODO: change this to Base ?
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Mat<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Mat<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Row<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Row<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Col<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Col<eT>& A)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap< Op<T1, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op<T1,op_htrans>& A)
max@0:     : M(A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Mat<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Mat<eT>, op_htrans>& A)
max@0:     : M(A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Row<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Row<eT>, op_htrans>& A)
max@0:     : M(A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Col<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Col<eT>, op_htrans>& A)
max@0:     : M(A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap< Op<T1, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op<T1,op_htrans2>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT      val;
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Mat<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Mat<eT>, op_htrans2>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Row<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Row<eT>, op_htrans2>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< Op< Col<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const Op< Col<eT>, op_htrans2>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap< eOp<T1, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap(const eOp<T1,eop_scalar_times>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT      val;
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< eOp<Mat<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const eOp<Mat<eT>,eop_scalar_times>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< eOp<Row<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const eOp<Row<eT>,eop_scalar_times>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap< eOp<Col<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap(const eOp<Col<eT>,eop_scalar_times>& A)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: //
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap_check
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline partial_unwrap_check(const T1& A, const Mat<eT>& B)
max@0:     : M(A)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Mat<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Mat<eT>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A == &B) ? new Mat<eT>(A) : 0 )
max@0:     , M       ( (&A == &B) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Row<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Row<eT>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A == &B) ? new Mat<eT>(A) : 0 )
max@0:     , M       ( (&A == &B) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Col<eT> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Col<eT>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A == &B) ? new Mat<eT>(A) : 0 )
max@0:     , M       ( (&A == &B) ? (*M_local)     : A )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap_check< Op<T1, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op<T1,op_htrans>& A, const Mat<eT>& B)
max@0:     : M(A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Mat<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Mat<eT>, op_htrans>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Row<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Row<eT>, op_htrans>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Col<eT>, op_htrans> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Col<eT>, op_htrans>& A, const Mat<eT>& B)
max@0:     : M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return eT(1); }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = false;
max@0:   
max@0:   // the order below is important
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap_check< Op<T1, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op<T1,op_htrans2>& A, const Mat<eT>& B)
max@0:     : val(A.aux)
max@0:     , M  (A.m)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT      val;
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Mat<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Mat<eT>, op_htrans2>& A, const Mat<eT>& B)
max@0:     : val     (A.aux)
max@0:     , M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   // the order below is important
max@0:   const eT       val;
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Row<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Row<eT>, op_htrans2>& A, const Mat<eT>& B)
max@0:     : val     (A.aux)
max@0:     , M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   // the order below is important
max@0:   const eT       val;
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< Op< Col<eT>, op_htrans2> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const Op< Mat<eT>, op_htrans2>& A, const Mat<eT>& B)
max@0:     : val     (A.aux)
max@0:     , M_local ( (&A.m == &B) ? new Mat<eT>(A.m) : 0   )
max@0:     , M       ( (&A.m == &B) ? (*M_local)       : A.m )
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     
max@0:     if(M_local)
max@0:       {
max@0:       delete M_local;
max@0:       }
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = true;
max@0:   static const bool do_times = true;
max@0:   
max@0:   // the order below is important
max@0:   const eT       val;
max@0:   const Mat<eT>* M_local;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename T1>
max@0: class partial_unwrap_check< eOp<T1, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const eOp<T1,eop_scalar_times>& A, const Mat<eT>& B)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT      val;
max@0:   const Mat<eT> M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< eOp<Mat<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const eOp<Mat<eT>,eop_scalar_times>& A, const Mat<eT>& B)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< eOp<Row<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const eOp<Row<eT>,eop_scalar_times>& A, const Mat<eT>& B)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: class partial_unwrap_check< eOp<Col<eT>, eop_scalar_times> >
max@0:   {
max@0:   public:
max@0:   
max@0:   inline
max@0:   partial_unwrap_check(const eOp<Col<eT>,eop_scalar_times>& A, const Mat<eT>& B)
max@0:     : val(A.aux)
max@0:     , M  (A.P.Q)
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     arma_ignore(B);
max@0:     }
max@0:   
max@0:   inline
max@0:   ~partial_unwrap_check()
max@0:     {
max@0:     arma_extra_debug_sigprint();
max@0:     }
max@0:   
max@0:   
max@0:   inline eT get_val() const { return val; }
max@0:   
max@0:   
max@0:   static const bool do_trans = false;
max@0:   static const bool do_times = true;
max@0:   
max@0:   const eT       val;
max@0:   const Mat<eT>& M;
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: //! @}