Chris@49: // Copyright (C) 2010-2013 NICTA (www.nicta.com.au)
Chris@49: // Copyright (C) 2010-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 Proxy
Chris@49: //! @{
Chris@49: 
Chris@49: 
Chris@49: // ea_type is the "element accessor" type,
Chris@49: // which can provide access to elements via operator[]
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_default
Chris@49:   {
Chris@49:   inline Proxy_default(const T1&)
Chris@49:     {
Chris@49:     arma_type_check(( is_arma_type<T1>::value == false ));
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_fixed
Chris@49:   {
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef T1                                       stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const T1&                                aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = true;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = T1::is_row;
Chris@49:   static const bool is_col = T1::is_col;
Chris@49:   
Chris@49:   arma_aligned const T1& Q;
Chris@49:   
Chris@49:   inline explicit Proxy_fixed(const T1& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline static uword get_n_rows() { return T1::n_rows; }
Chris@49:   arma_inline static uword get_n_cols() { return T1::n_cols; }
Chris@49:   arma_inline static uword get_n_elem() { return T1::n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&Q) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const
Chris@49:     {
Chris@49:     #if defined(ARMA_HAVE_ALIGNED_ATTRIBUTE)
Chris@49:       return true;
Chris@49:     #else
Chris@49:       return memory::is_aligned(Q.memptr());
Chris@49:     #endif
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct Proxy_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_redirect<T1, false> { typedef Proxy_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_redirect<T1, true>  { typedef Proxy_fixed<T1>   result; };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy : public Proxy_redirect<T1, is_Mat_fixed<T1>::value >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   inline Proxy(const T1& A)
Chris@49:     : Proxy_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: class Proxy< Mat<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Mat<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const Mat<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Mat<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&Q) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< Col<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Col<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Col<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const Col<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Col<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];        }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row];      }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q.at_alt(i); }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&Q) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< Row<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Row<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Row<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = true;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const Row<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Row<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return 1;        }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];        }
Chris@49:   arma_inline elem_type at         (const uword, const uword col) const { return Q[col];      }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q.at_alt(i); }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&Q) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, typename gen_type>
Chris@49: class Proxy< Gen<T1, gen_type > >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Gen<T1, gen_type>                        stored_type;
Chris@49:   typedef const Gen<T1, gen_type>&                 ea_type;
Chris@49:   typedef const Gen<T1, gen_type>&                 aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = Gen<T1, gen_type>::prefer_at_accessor;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = Gen<T1, gen_type>::is_row;
Chris@49:   static const bool is_col = Gen<T1, gen_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const Gen<T1, gen_type>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Gen<T1, gen_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return (is_row ? 1 : Q.n_rows);                           }
Chris@49:   arma_inline uword get_n_cols() const { return (is_col ? 1 : Q.n_cols);                           }
Chris@49:   arma_inline uword get_n_elem() const { return (is_row ? 1 : Q.n_rows) * (is_col ? 1 : Q.n_cols); }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q[i];           }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return Gen<T1, gen_type>::is_simple; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, typename op_type>
Chris@49: class Proxy< Op<T1, op_type> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = Op<T1, op_type>::is_row;
Chris@49:   static const bool is_col = Op<T1, op_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const Mat<elem_type> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<T1, op_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy_diagvec_mat
Chris@49:   {
Chris@49:   inline Proxy_diagvec_mat(const T1&) {}
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy_diagvec_mat< Op<T1, op_diagvec> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef diagview<elem_type>                      stored_type;
Chris@49:   typedef const diagview<elem_type>&               ea_type;
Chris@49:   typedef const diagview<elem_type>&               aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const Mat<elem_type>&     R;
Chris@49:   arma_aligned const diagview<elem_type> Q;
Chris@49:   
Chris@49:   inline explicit Proxy_diagvec_mat(const Op<T1, op_diagvec>& A)
Chris@49:     : R(A.m), Q( R.diag( (A.aux_uword_b > 0) ? -sword(A.aux_uword_a) : sword(A.aux_uword_a) ) )
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];         }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q.at(row, 0); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];         }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&R) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy_diagvec_expr
Chris@49:   {
Chris@49:   inline Proxy_diagvec_expr(const T1&) {}
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy_diagvec_expr< Op<T1, op_diagvec> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const Mat<elem_type> Q;
Chris@49:   
Chris@49:   inline explicit Proxy_diagvec_expr(const Op<T1, op_diagvec>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];         }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q.at(row, 0); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q.at_alt(i);  }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct Proxy_diagvec_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_diagvec_redirect< Op<T1, op_diagvec>, true > { typedef Proxy_diagvec_mat < Op<T1, op_diagvec> > result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_diagvec_redirect< Op<T1, op_diagvec>, false> { typedef Proxy_diagvec_expr< Op<T1, op_diagvec> > result; };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy< Op<T1, op_diagvec> >
Chris@49:   : public Proxy_diagvec_redirect< Op<T1, op_diagvec>, is_Mat<T1>::value >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename Proxy_diagvec_redirect< Op<T1, op_diagvec>, is_Mat<T1>::value >::result Proxy_diagvec;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<T1, op_diagvec>& A)
Chris@49:     : Proxy_diagvec(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_xtrans_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   arma_aligned const Mat<eT> Q;
Chris@49:   
Chris@49:   arma_hot
Chris@49:   inline Proxy_xtrans_default(const T1& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
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 Proxy_xtrans_vector
Chris@49:   {
Chris@49:   inline Proxy_xtrans_vector(const T1&) {}
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_xtrans_vector< Op<T1, op_htrans> >
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = quasi_unwrap<T1>::has_subview;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = true;
Chris@49:   
Chris@49:   arma_aligned const quasi_unwrap<T1> U; // avoid copy if T1 is a Row, Col or subview_col
Chris@49:   arma_aligned const Mat<eT>          Q;
Chris@49:   
Chris@49:   inline Proxy_xtrans_vector(const Op<T1, op_htrans>& A)
Chris@49:     : U(A.m)
Chris@49:     , Q(const_cast<eT*>(U.M.memptr()), U.M.n_cols, U.M.n_rows, false, false)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return U.is_alias(X); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_xtrans_vector< Op<T1, op_strans> >
Chris@49:   {
Chris@49:   typedef typename T1::elem_type eT;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = quasi_unwrap<T1>::has_subview;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = true;
Chris@49:   
Chris@49:   arma_aligned const quasi_unwrap<T1> U; // avoid copy if T1 is a Row, Col or subview_col
Chris@49:   arma_aligned const Mat<eT>          Q;
Chris@49:   
Chris@49:   inline Proxy_xtrans_vector(const Op<T1, op_strans>& A)
Chris@49:     : U(A.m)
Chris@49:     , Q(const_cast<eT*>(U.M.memptr()), U.M.n_cols, U.M.n_rows, false, false)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return U.is_alias(X); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct Proxy_xtrans_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_xtrans_redirect<T1, false> { typedef Proxy_xtrans_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_xtrans_redirect<T1, true>  { typedef Proxy_xtrans_vector<T1>  result; };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy< Op<T1, op_htrans> >
Chris@49:   : public
Chris@49:     Proxy_xtrans_redirect
Chris@49:       <
Chris@49:       Op<T1, op_htrans>,
Chris@49:       ((is_complex<typename T1::elem_type>::value == false) && ((Op<T1, op_htrans>::is_row) || (Op<T1, op_htrans>::is_col)) )
Chris@49:       >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef
Chris@49:   typename
Chris@49:   Proxy_xtrans_redirect
Chris@49:     <
Chris@49:     Op<T1, op_htrans>,
Chris@49:     ((is_complex<typename T1::elem_type>::value == false) && ((Op<T1, op_htrans>::is_row) || (Op<T1, op_htrans>::is_col)) )
Chris@49:     >::result
Chris@49:   Proxy_xtrans;
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = Proxy_xtrans::prefer_at_accessor;
Chris@49:   static const bool has_subview        = Proxy_xtrans::has_subview;
Chris@49:   static const bool is_fixed           = Proxy_xtrans::is_fixed;
Chris@49:   static const bool fake_mat           = Proxy_xtrans::fake_mat;
Chris@49:   
Chris@49:   // NOTE: the Op class takes care of swapping row and col for op_htrans
Chris@49:   static const bool is_row = Op<T1, op_htrans>::is_row;
Chris@49:   static const bool is_col = Op<T1, op_htrans>::is_col;
Chris@49:   
Chris@49:   using Proxy_xtrans::Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<T1, op_htrans>& A)
Chris@49:     : Proxy_xtrans(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Proxy_xtrans::is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy< Op<T1, op_strans> >
Chris@49:   : public
Chris@49:     Proxy_xtrans_redirect
Chris@49:       <
Chris@49:       Op<T1, op_strans>,
Chris@49:       ( (Op<T1, op_strans>::is_row) || (Op<T1, op_strans>::is_col) )
Chris@49:       >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef
Chris@49:   typename
Chris@49:   Proxy_xtrans_redirect
Chris@49:     <
Chris@49:     Op<T1, op_strans>,
Chris@49:     ( (Op<T1, op_strans>::is_row) || (Op<T1, op_strans>::is_col) )
Chris@49:     >::result
Chris@49:   Proxy_xtrans;
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = Proxy_xtrans::prefer_at_accessor;
Chris@49:   static const bool has_subview        = Proxy_xtrans::has_subview;
Chris@49:   static const bool is_fixed           = Proxy_xtrans::is_fixed;
Chris@49:   static const bool fake_mat           = Proxy_xtrans::fake_mat;
Chris@49:   
Chris@49:   // NOTE: the Op class takes care of swapping row and col for op_strans
Chris@49:   static const bool is_row = Op<T1, op_strans>::is_row;
Chris@49:   static const bool is_col = Op<T1, op_strans>::is_col;
Chris@49:   
Chris@49:   using Proxy_xtrans::Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<T1, op_strans>& A)
Chris@49:     : Proxy_xtrans(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Proxy_xtrans::is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct Proxy_subview_row_htrans_cx
Chris@49:   {
Chris@49:   typedef eT                                elem_type;
Chris@49:   typedef typename get_pod_type<eT>::result pod_type;
Chris@49:   typedef subview_row_htrans<eT>            stored_type;
Chris@49:   typedef const subview_row_htrans<eT>&     ea_type;
Chris@49:   typedef const subview_row_htrans<eT>&     aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_row_htrans<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy_subview_row_htrans_cx(const Op<subview_row<eT>, op_htrans>& A)
Chris@49:     : Q(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.sv_row.m)) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct Proxy_subview_row_htrans_non_cx
Chris@49:   {
Chris@49:   typedef eT                                elem_type;
Chris@49:   typedef typename get_pod_type<eT>::result pod_type;
Chris@49:   typedef subview_row_strans<eT>            stored_type;
Chris@49:   typedef const subview_row_strans<eT>&     ea_type;
Chris@49:   typedef const subview_row_strans<eT>&     aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_row_strans<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy_subview_row_htrans_non_cx(const Op<subview_row<eT>, op_htrans>& A)
Chris@49:     : Q(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.sv_row.m)) == void_ptr(&X)); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT, bool condition>
Chris@49: struct Proxy_subview_row_htrans_redirect {};
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct Proxy_subview_row_htrans_redirect<eT, true>  { typedef Proxy_subview_row_htrans_cx<eT>      result; };
Chris@49: 
Chris@49: template<typename eT>
Chris@49: struct Proxy_subview_row_htrans_redirect<eT, false> { typedef Proxy_subview_row_htrans_non_cx<eT>  result; };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< Op<subview_row<eT>, op_htrans> >
Chris@49:   : public
Chris@49:     Proxy_subview_row_htrans_redirect
Chris@49:       <
Chris@49:       eT,
Chris@49:       is_complex<eT>::value
Chris@49:       >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef
Chris@49:   typename
Chris@49:   Proxy_subview_row_htrans_redirect
Chris@49:       <
Chris@49:       eT,
Chris@49:       is_complex<eT>::value
Chris@49:       >::result
Chris@49:   Proxy_sv_row_ht;
Chris@49:   
Chris@49:   typedef typename Proxy_sv_row_ht::elem_type   elem_type;
Chris@49:   typedef typename Proxy_sv_row_ht::pod_type    pod_type;
Chris@49:   typedef typename Proxy_sv_row_ht::stored_type stored_type;
Chris@49:   typedef typename Proxy_sv_row_ht::ea_type     ea_type;
Chris@49:   typedef typename Proxy_sv_row_ht::ea_type     aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = Proxy_sv_row_ht::prefer_at_accessor;
Chris@49:   static const bool has_subview        = Proxy_sv_row_ht::has_subview;
Chris@49:   static const bool is_fixed           = Proxy_sv_row_ht::is_fixed;
Chris@49:   static const bool fake_mat           = Proxy_sv_row_ht::fake_mat;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   using Proxy_sv_row_ht::Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<subview_row<eT>, op_htrans>& A)
Chris@49:     : Proxy_sv_row_ht(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Proxy_sv_row_ht::is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< Op<subview_row<eT>, op_strans> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                elem_type;
Chris@49:   typedef typename get_pod_type<eT>::result pod_type;
Chris@49:   typedef subview_row_strans<eT>            stored_type;
Chris@49:   typedef const subview_row_strans<eT>&     ea_type;
Chris@49:   typedef const subview_row_strans<eT>&     aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_row_strans<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<subview_row<eT>, op_strans>& A)
Chris@49:     : Q(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.sv_row.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49:   
Chris@49: template<typename T>
Chris@49: class Proxy< Op< Row< std::complex<T> >, op_htrans> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename std::complex<T>  eT;
Chris@49: 
Chris@49:   typedef typename std::complex<T>  elem_type;
Chris@49:   typedef T                         pod_type;
Chris@49:   typedef xvec_htrans<eT>           stored_type;
Chris@49:   typedef const xvec_htrans<eT>&    ea_type;
Chris@49:   typedef const xvec_htrans<eT>&    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   const xvec_htrans<eT> Q;
Chris@49:   const Row<eT>&        src;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op< Row< std::complex<T> >, op_htrans>& A)
Chris@49:     : Q  (A.m.memptr(), A.m.n_rows, A.m.n_cols)
Chris@49:     , src(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
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:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T>
Chris@49: class Proxy< Op< Col< std::complex<T> >, op_htrans> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename std::complex<T>  eT;
Chris@49: 
Chris@49:   typedef typename std::complex<T>  elem_type;
Chris@49:   typedef T                         pod_type;
Chris@49:   typedef xvec_htrans<eT>           stored_type;
Chris@49:   typedef const xvec_htrans<eT>&    ea_type;
Chris@49:   typedef const xvec_htrans<eT>&    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = true;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   const xvec_htrans<eT> Q;
Chris@49:   const Col<eT>&        src;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op< Col< std::complex<T> >, op_htrans>& A)
Chris@49:     : Q  (A.m.memptr(), A.m.n_rows, A.m.n_cols)
Chris@49:     , src(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return 1;        }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword, const uword col) const { return Q[col]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
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:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T>
Chris@49: class Proxy< Op< subview_col< std::complex<T> >, op_htrans> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename std::complex<T>  eT;
Chris@49:   
Chris@49:   typedef typename std::complex<T>  elem_type;
Chris@49:   typedef T                         pod_type;
Chris@49:   typedef xvec_htrans<eT>           stored_type;
Chris@49:   typedef const xvec_htrans<eT>&    ea_type;
Chris@49:   typedef const xvec_htrans<eT>&    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = true;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   const xvec_htrans<eT>  Q;
Chris@49:   const subview_col<eT>& src;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op< subview_col< std::complex<T> >, op_htrans>& A)
Chris@49:     : Q  (A.m.colptr(0), A.m.n_rows, A.m.n_cols)
Chris@49:     , src(A.m)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return 1;        }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword, const uword col) const { return Q[col]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return void_ptr(&src.m) == void_ptr(&X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_htrans2_default
Chris@49:   {
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const Mat<elem_type> Q;
Chris@49:   
Chris@49:   arma_hot
Chris@49:   inline Proxy_htrans2_default(const T1& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_htrans2_vector
Chris@49:   {
Chris@49:   inline Proxy_htrans2_vector(const T1&) {}
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_htrans2_vector< Op<T1, op_htrans2> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                           elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result         pod_type;
Chris@49:   typedef       eOp< Op<T1, op_htrans>, eop_scalar_times>  stored_type;
Chris@49:   typedef const eOp< Op<T1, op_htrans>, eop_scalar_times>& ea_type;
Chris@49:   typedef const eOp< Op<T1, op_htrans>, eop_scalar_times>& aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = eOp< Op<T1, op_htrans>, eop_scalar_times>::prefer_at_accessor;
Chris@49:   static const bool has_subview        = eOp< Op<T1, op_htrans>, eop_scalar_times>::has_subview;
Chris@49:   static const bool is_fixed           = eOp< Op<T1, op_htrans>, eop_scalar_times>::is_fixed;
Chris@49:   static const bool fake_mat           = eOp< Op<T1, op_htrans>, eop_scalar_times>::fake_mat;
Chris@49:   
Chris@49:   // NOTE: the Op class takes care of swapping row and col for op_htrans
Chris@49:   static const bool is_row = eOp< Op<T1, op_htrans>, eop_scalar_times>::is_row;
Chris@49:   static const bool is_col = eOp< Op<T1, op_htrans>, eop_scalar_times>::is_col;
Chris@49:   
Chris@49:   arma_aligned const      Op<T1, op_htrans>                     R;
Chris@49:   arma_aligned const eOp< Op<T1, op_htrans>, eop_scalar_times > Q;
Chris@49:   
Chris@49:   inline explicit Proxy_htrans2_vector(const Op<T1, op_htrans2>& A)
Chris@49:     : R(A.m)
Chris@49:     , Q(R, A.aux)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.get_n_rows(); }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.get_n_cols(); }
Chris@49:   arma_inline uword get_n_elem() const { return Q.get_n_elem();              }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Q.P.is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return Q.P.is_aligned(); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, bool condition>
Chris@49: struct Proxy_htrans2_redirect {};
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_htrans2_redirect<T1, false> { typedef Proxy_htrans2_default<T1> result; };
Chris@49: 
Chris@49: template<typename T1>
Chris@49: struct Proxy_htrans2_redirect<T1, true>  { typedef Proxy_htrans2_vector<T1>  result; };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1>
Chris@49: class Proxy< Op<T1, op_htrans2> >
Chris@49:   : public
Chris@49:     Proxy_htrans2_redirect
Chris@49:       <
Chris@49:       Op<T1, op_htrans2>,
Chris@49:       ( (Op<T1, op_htrans2>::is_row) || (Op<T1, op_htrans2>::is_col) )
Chris@49:       >::result
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef
Chris@49:   typename
Chris@49:   Proxy_htrans2_redirect
Chris@49:     <
Chris@49:     Op<T1, op_htrans2>,
Chris@49:     ( (Op<T1, op_htrans2>::is_row) || (Op<T1, op_htrans2>::is_col) )
Chris@49:     >::result
Chris@49:   Proxy_htrans2;
Chris@49:   
Chris@49:   typedef typename Proxy_htrans2::elem_type       elem_type;
Chris@49:   typedef typename Proxy_htrans2::pod_type        pod_type;
Chris@49:   typedef typename Proxy_htrans2::stored_type     stored_type;
Chris@49:   typedef typename Proxy_htrans2::ea_type         ea_type;
Chris@49:   typedef typename Proxy_htrans2::aligned_ea_type aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = Proxy_htrans2::prefer_at_accessor;
Chris@49:   static const bool has_subview        = Proxy_htrans2::has_subview;
Chris@49:   static const bool is_fixed           = Proxy_htrans2::is_fixed;
Chris@49:   static const bool fake_mat           = Proxy_htrans2::fake_mat;
Chris@49:   
Chris@49:   static const bool is_row = Proxy_htrans2::is_row;
Chris@49:   static const bool is_col = Proxy_htrans2::is_col;
Chris@49:   
Chris@49:   using Proxy_htrans2::Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Op<T1, op_htrans2>& A)
Chris@49:     : Proxy_htrans2(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Proxy_htrans2::get_n_rows(); }
Chris@49:   arma_inline uword get_n_cols() const { return Proxy_htrans2::get_n_cols(); }
Chris@49:   arma_inline uword get_n_elem() const { return Proxy_htrans2::get_n_elem(); }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Proxy_htrans2::get_ea();         }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Proxy_htrans2::get_aligned_ea(); }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Proxy_htrans2::is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return Proxy_htrans2::is_aligned(); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, typename T2, typename glue_type>
Chris@49: class Proxy< Glue<T1, T2, glue_type> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<elem_type>                           stored_type;
Chris@49:   typedef const elem_type*                         ea_type;
Chris@49:   typedef const Mat<elem_type>&                    aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = Glue<T1, T2, glue_type>::is_row;
Chris@49:   static const bool is_col = Glue<T1, T2, glue_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const Mat<elem_type> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const Glue<T1, T2, glue_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49: 
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< subview<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef subview<eT>                              stored_type;
Chris@49:   typedef const subview<eT>&                       ea_type;
Chris@49:   typedef const subview<eT>&                       aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = true;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const subview<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q[i];           }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< subview_col<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef subview_col<eT>                          stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const subview_col<eT>&                   aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_col<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_col<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];        }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row];      }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q.at_alt(i); }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.colmem; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;        }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.colmem); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< subview_row<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef subview_row<eT>                          stored_type;
Chris@49:   typedef const subview_row<eT>&                   ea_type;
Chris@49:   typedef const subview_row<eT>&                   aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = true;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const subview_row<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_row<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return 1;        }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword, const uword col) const { return Q[col]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< subview_row_strans<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef subview_row_strans<eT>                   stored_type;
Chris@49:   typedef const subview_row_strans<eT>&            ea_type;
Chris@49:   typedef const subview_row_strans<eT>&            aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_row_strans<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_row_strans<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.sv_row.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< subview_row_htrans<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef subview_row_htrans<eT>                   stored_type;
Chris@49:   typedef const subview_row_htrans<eT>&            ea_type;
Chris@49:   typedef const subview_row_htrans<eT>&            aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const subview_row_htrans<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_row_htrans<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];   }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row]; }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.sv_row.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< xvec_htrans<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Mat<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const Mat<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const xvec_htrans<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];          }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row,col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT, typename T1>
Chris@49: class Proxy< subview_elem1<eT,T1> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Mat<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const Mat<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_elem1<eT,T1>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];        }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q[row];      }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q.at_alt(i); }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT, typename T1, typename T2>
Chris@49: class Proxy< subview_elem2<eT,T1,T2> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef Mat<eT>                                  stored_type;
Chris@49:   typedef const eT*                                ea_type;
Chris@49:   typedef const Mat<eT>&                           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = false;
Chris@49:   
Chris@49:   arma_aligned const Mat<eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const subview_elem2<eT,T1,T2>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename eT>
Chris@49: class Proxy< diagview<eT> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef eT                                       elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef diagview<eT>                             stored_type;
Chris@49:   typedef const diagview<eT>&                      ea_type;
Chris@49:   typedef const diagview<eT>&                      aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = true;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = false;
Chris@49:   static const bool is_col = true;
Chris@49:   
Chris@49:   arma_aligned const diagview<eT>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const diagview<eT>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return 1;        }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem; }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                const { return Q[i];         }
Chris@49:   arma_inline elem_type at         (const uword row, const uword) const { return Q.at(row, 0); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                const { return Q[i];         }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (void_ptr(&(Q.m)) == void_ptr(&X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return false; }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, typename eop_type>
Chris@49: class Proxy< eOp<T1, eop_type > >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef eOp<T1, eop_type>                        stored_type;
Chris@49:   typedef const eOp<T1, eop_type>&                 ea_type;
Chris@49:   typedef const eOp<T1, eop_type>&                 aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = eOp<T1, eop_type>::prefer_at_accessor;
Chris@49:   static const bool has_subview        = eOp<T1, eop_type>::has_subview;
Chris@49:   static const bool is_fixed           = eOp<T1, eop_type>::is_fixed;
Chris@49:   static const bool fake_mat           = eOp<T1, eop_type>::fake_mat;
Chris@49:   
Chris@49:   static const bool is_row = eOp<T1, eop_type>::is_row;
Chris@49:   static const bool is_col = eOp<T1, eop_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const eOp<T1, eop_type>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const eOp<T1, eop_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.get_n_rows(); }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.get_n_cols(); }
Chris@49:   arma_inline uword get_n_elem() const { return Q.get_n_elem();              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return Q.P.is_alias(X); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return Q.P.is_aligned(); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename T1, typename T2, typename eglue_type>
Chris@49: class Proxy< eGlue<T1, T2, eglue_type > >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef typename T1::elem_type                   elem_type;
Chris@49:   typedef typename get_pod_type<elem_type>::result pod_type;
Chris@49:   typedef eGlue<T1, T2, eglue_type>                stored_type;
Chris@49:   typedef const eGlue<T1, T2, eglue_type>&         ea_type;
Chris@49:   typedef const eGlue<T1, T2, eglue_type>&         aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = eGlue<T1, T2, eglue_type>::prefer_at_accessor;
Chris@49:   static const bool has_subview        = eGlue<T1, T2, eglue_type>::has_subview;
Chris@49:   static const bool is_fixed           = eGlue<T1, T2, eglue_type>::is_fixed;
Chris@49:   static const bool fake_mat           = eGlue<T1, T2, eglue_type>::fake_mat;
Chris@49:   
Chris@49:   static const bool is_row = eGlue<T1, T2, eglue_type>::is_row;
Chris@49:   static const bool is_col = eGlue<T1, T2, eglue_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const eGlue<T1, T2, eglue_type>& Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const eGlue<T1, T2, eglue_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.get_n_rows(); }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.get_n_cols(); }
Chris@49:   arma_inline uword get_n_elem() const { return Q.get_n_elem();              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];           }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row, col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);    }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q; }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q; }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>& X) const { return (Q.P1.is_alias(X) || Q.P2.is_alias(X)); }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return (Q.P1.is_aligned() && Q.P2.is_aligned()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename op_type>
Chris@49: class Proxy< mtOp<out_eT, T1, op_type> >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef          out_eT                       elem_type;
Chris@49:   typedef typename get_pod_type<out_eT>::result pod_type;
Chris@49:   typedef          Mat<out_eT>                  stored_type;
Chris@49:   typedef          const elem_type*             ea_type;
Chris@49:   typedef          const Mat<out_eT>&           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = mtOp<out_eT, T1, op_type>::is_row;
Chris@49:   static const bool is_col = mtOp<out_eT, T1, op_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const Mat<out_eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const mtOp<out_eT, T1, op_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];          }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row,col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: template<typename out_eT, typename T1, typename T2, typename glue_type>
Chris@49: class Proxy< mtGlue<out_eT, T1, T2, glue_type > >
Chris@49:   {
Chris@49:   public:
Chris@49:   
Chris@49:   typedef          out_eT                       elem_type;
Chris@49:   typedef typename get_pod_type<out_eT>::result pod_type;
Chris@49:   typedef          Mat<out_eT>                  stored_type;
Chris@49:   typedef          const elem_type*             ea_type;
Chris@49:   typedef          const Mat<out_eT>&           aligned_ea_type;
Chris@49:   
Chris@49:   static const bool prefer_at_accessor = false;
Chris@49:   static const bool has_subview        = false;
Chris@49:   static const bool is_fixed           = false;
Chris@49:   static const bool fake_mat           = false;
Chris@49:   
Chris@49:   static const bool is_row = mtGlue<out_eT, T1, T2, glue_type>::is_row;
Chris@49:   static const bool is_col = mtGlue<out_eT, T1, T2, glue_type>::is_col;
Chris@49:   
Chris@49:   arma_aligned const Mat<out_eT> Q;
Chris@49:   
Chris@49:   inline explicit Proxy(const mtGlue<out_eT, T1, T2, glue_type>& A)
Chris@49:     : Q(A)
Chris@49:     {
Chris@49:     arma_extra_debug_sigprint();
Chris@49:     }
Chris@49:   
Chris@49:   arma_inline uword get_n_rows() const { return is_row ? 1 : Q.n_rows; }
Chris@49:   arma_inline uword get_n_cols() const { return is_col ? 1 : Q.n_cols; }
Chris@49:   arma_inline uword get_n_elem() const { return Q.n_elem;              }
Chris@49:   
Chris@49:   arma_inline elem_type operator[] (const uword i)                    const { return Q[i];          }
Chris@49:   arma_inline elem_type at         (const uword row, const uword col) const { return Q.at(row,col); }
Chris@49:   arma_inline elem_type at_alt     (const uword i)                    const { return Q.at_alt(i);   }
Chris@49:   
Chris@49:   arma_inline         ea_type         get_ea() const { return Q.memptr(); }
Chris@49:   arma_inline aligned_ea_type get_aligned_ea() const { return Q;          }
Chris@49:   
Chris@49:   template<typename eT2>
Chris@49:   arma_inline bool is_alias(const Mat<eT2>&) const { return false; }
Chris@49:   
Chris@49:   arma_inline bool is_aligned() const { return memory::is_aligned(Q.memptr()); }
Chris@49:   };
Chris@49: 
Chris@49: 
Chris@49: 
Chris@49: //! @}