max@0: // Copyright (C) 2010-2011 NICTA (www.nicta.com.au)
max@0: // Copyright (C) 2010-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 eglue_core
max@0: //! @{
max@0: 
max@0: 
max@0: 
max@0: class eglue_plus : public eglue_core<eglue_plus>
max@0:   {
max@0:   public:
max@0:   
max@0:   inline static const char* text() { return "addition"; }
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: class eglue_minus : public eglue_core<eglue_minus>
max@0:   {
max@0:   public:
max@0:   
max@0:   inline static const char* text() { return "subtraction"; }
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: class eglue_div : public eglue_core<eglue_div>
max@0:   {
max@0:   public:
max@0:   
max@0:   inline static const char* text() { return "element-wise division"; }
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: class eglue_schur : public eglue_core<eglue_schur>
max@0:   {
max@0:   public:
max@0:   
max@0:   inline static const char* text() { return "element-wise multiplication"; }
max@0:   };
max@0: 
max@0: 
max@0: 
max@0: #undef arma_applier_1
max@0: #undef arma_applier_2
max@0: #undef arma_applier_3
max@0: #undef operatorA
max@0: #undef operatorB
max@0: 
max@0: #define arma_applier_1(operatorA, operatorB) \
max@0:   {\
max@0:   uword i,j;\
max@0:   \
max@0:   for(i=0, j=1; j<n_elem; i+=2, j+=2)\
max@0:     {\
max@0:     eT tmp_i = P1[i];\
max@0:     eT tmp_j = P1[j];\
max@0:     \
max@0:     tmp_i operatorB##= P2[i];\
max@0:     tmp_j operatorB##= P2[j];\
max@0:     \
max@0:     out_mem[i] operatorA tmp_i;\
max@0:     out_mem[j] operatorA tmp_j;\
max@0:     }\
max@0:   \
max@0:   if(i < n_elem)\
max@0:     {\
max@0:     out_mem[i] operatorA P1[i] operatorB P2[i];\
max@0:     }\
max@0:   }
max@0:   
max@0: 
max@0: 
max@0: #define arma_applier_2(operatorA, operatorB) \
max@0:   {\
max@0:   uword count = 0;\
max@0:   \
max@0:   for(uword col=0; col<n_cols; ++col)\
max@0:     {\
max@0:     uword i,j;\
max@0:     \
max@0:     for(i=0, j=1; j<n_rows; i+=2, j+=2, count+=2)\
max@0:       {\
max@0:       eT tmp_i = P1.at(i,col);\
max@0:       eT tmp_j = P1.at(j,col);\
max@0:       \
max@0:       tmp_i operatorB##= P2.at(i,col);\
max@0:       tmp_j operatorB##= P2.at(j,col);\
max@0:       \
max@0:       out_mem[count  ] operatorA tmp_i;\
max@0:       out_mem[count+1] operatorA tmp_j;\
max@0:       }\
max@0:     \
max@0:     if(i < n_rows)\
max@0:       {\
max@0:       out_mem[count] operatorA P1.at(i,col) operatorB P2.at(i,col);\
max@0:       ++count;\
max@0:       }\
max@0:     }\
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: #define arma_applier_3(operatorA, operatorB) \
max@0:   {\
max@0:   uword count = 0;\
max@0:   \
max@0:   for(uword slice=0; slice<n_slices; ++slice)\
max@0:     {\
max@0:     for(uword col=0; col<n_cols; ++col)\
max@0:       {\
max@0:       uword i,j;\
max@0:       \
max@0:       for(i=0, j=1; j<n_rows; i+=2, j+=2, count+=2)\
max@0:         {\
max@0:         eT tmp_i = P1.at(i,col,slice);\
max@0:         eT tmp_j = P1.at(j,col,slice);\
max@0:         \
max@0:         tmp_i operatorB##= P2.at(i,col,slice);\
max@0:         tmp_j operatorB##= P2.at(j,col,slice);\
max@0:         \
max@0:         out_mem[count  ] operatorA tmp_i;\
max@0:         out_mem[count+1] operatorA tmp_j;\
max@0:         }\
max@0:       \
max@0:       if(i < n_rows)\
max@0:         {\
max@0:         out_mem[count] operatorA P1.at(i,col,slice) operatorB P2.at(i,col,slice);\
max@0:         ++count;\
max@0:         }\
max@0:       }\
max@0:     }\
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: //
max@0: // matrices
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply(Mat<typename T1::elem_type>& out, const eGlue<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   const bool prefer_at_accessor = (Proxy<T1>::prefer_at_accessor || Proxy<T2>::prefer_at_accessor);
max@0:   
max@0:   // NOTE: we're assuming that the matrix has already been set to the correct size and there is no aliasing;
max@0:   // size setting and alias checking is done by either the Mat contructor or operator=()
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename Proxy<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename Proxy<T2>::ea_type P2 = x.P2.get_ea();
max@0:   
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows = out.n_rows;
max@0:     const uword n_cols = out.n_cols;
max@0:   
max@0:     const Proxy<T1>& P1 = x.P1;
max@0:     const Proxy<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_2(=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_2(=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_2(=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_2(=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_plus(Mat<typename T1::elem_type>& out, const eGlue<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_debug_assert_same_size(out, x.P1, "addition");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (Proxy<T1>::prefer_at_accessor || Proxy<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename Proxy<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename Proxy<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(+=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(+=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(+=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(+=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows = out.n_rows;
max@0:     const uword n_cols = out.n_cols;
max@0:     
max@0:     const Proxy<T1>& P1 = x.P1;
max@0:     const Proxy<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_2(+=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_2(+=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_2(+=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_2(+=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_minus(Mat<typename T1::elem_type>& out, const eGlue<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_debug_assert_same_size(out, x.P1, "subtraction");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (Proxy<T1>::prefer_at_accessor || Proxy<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename Proxy<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename Proxy<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(-=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(-=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(-=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(-=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows = out.n_rows;
max@0:     const uword n_cols = out.n_cols;
max@0:     
max@0:     const Proxy<T1>& P1 = x.P1;
max@0:     const Proxy<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_2(-=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_2(-=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_2(-=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_2(-=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_schur(Mat<typename T1::elem_type>& out, const eGlue<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_debug_assert_same_size(out, x.P1, "element-wise multiplication");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (Proxy<T1>::prefer_at_accessor || Proxy<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename Proxy<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename Proxy<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(*=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(*=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(*=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(*=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows = out.n_rows;
max@0:     const uword n_cols = out.n_cols;
max@0:     
max@0:     const Proxy<T1>& P1 = x.P1;
max@0:     const Proxy<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_2(*=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_2(*=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_2(*=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_2(*=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_div(Mat<typename T1::elem_type>& out, const eGlue<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_debug_assert_same_size(out, x.P1, "element-wise division");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (Proxy<T1>::prefer_at_accessor || Proxy<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename Proxy<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename Proxy<T2>::ea_type P2 = x.P2.get_ea();
max@0:   
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(/=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(/=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(/=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(/=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows = out.n_rows;
max@0:     const uword n_cols = out.n_cols;
max@0:     
max@0:     const Proxy<T1>& P1 = x.P1;
max@0:     const Proxy<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_2(*=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_2(*=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_2(*=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_2(*=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: //
max@0: // cubes
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply(Cube<typename T1::elem_type>& out, const eGlueCube<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   const bool prefer_at_accessor = (ProxyCube<T1>::prefer_at_accessor || ProxyCube<T2>::prefer_at_accessor);
max@0:   
max@0:   // NOTE: we're assuming that the cube has already been set to the correct size and there is no aliasing;
max@0:   // size setting and alias checking is done by either the Cube contructor or operator=()
max@0:   
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem = out.n_elem;
max@0:     
max@0:     typename ProxyCube<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename ProxyCube<T2>::ea_type P2 = x.P2.get_ea();
max@0:   
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const uword n_rows   = x.get_n_rows();
max@0:     const uword n_cols   = x.get_n_cols();
max@0:     const uword n_slices = x.get_n_slices();
max@0:   
max@0:     const ProxyCube<T1>& P1 = x.P1;
max@0:     const ProxyCube<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_3(=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_3(=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_3(=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_3(=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_plus(Cube<typename T1::elem_type>& out, const eGlueCube<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword n_rows   = x.get_n_rows();
max@0:   const uword n_cols   = x.get_n_cols();
max@0:   const uword n_slices = x.get_n_slices();
max@0:   
max@0:   arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "addition");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (ProxyCube<T1>::prefer_at_accessor || ProxyCube<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem  = out.n_elem;
max@0:     
max@0:     typename ProxyCube<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename ProxyCube<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(+=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(+=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(+=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(+=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const ProxyCube<T1>& P1 = x.P1;
max@0:     const ProxyCube<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_3(+=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_3(+=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_3(+=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_3(+=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_minus(Cube<typename T1::elem_type>& out, const eGlueCube<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword n_rows   = x.get_n_rows();
max@0:   const uword n_cols   = x.get_n_cols();
max@0:   const uword n_slices = x.get_n_slices();
max@0:   
max@0:   arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "subtraction");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (ProxyCube<T1>::prefer_at_accessor || ProxyCube<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem  = out.n_elem;
max@0:     
max@0:     typename ProxyCube<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename ProxyCube<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(-=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(-=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(-=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(-=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const ProxyCube<T1>& P1 = x.P1;
max@0:     const ProxyCube<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_3(-=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_3(-=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_3(-=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_3(-=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_schur(Cube<typename T1::elem_type>& out, const eGlueCube<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword n_rows   = x.get_n_rows();
max@0:   const uword n_cols   = x.get_n_cols();
max@0:   const uword n_slices = x.get_n_slices();
max@0:   
max@0:   arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "element-wise multiplication");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (ProxyCube<T1>::prefer_at_accessor || ProxyCube<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem  = out.n_elem;
max@0:     
max@0:     typename ProxyCube<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename ProxyCube<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(*=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(*=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(*=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(*=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const ProxyCube<T1>& P1 = x.P1;
max@0:     const ProxyCube<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_3(*=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_3(*=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_3(*=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_3(*=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eglue_type>
max@0: template<typename T1, typename T2>
max@0: arma_hot
max@0: inline
max@0: void
max@0: eglue_core<eglue_type>::apply_inplace_div(Cube<typename T1::elem_type>& out, const eGlueCube<T1, T2, eglue_type>& x)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword n_rows   = x.get_n_rows();
max@0:   const uword n_cols   = x.get_n_cols();
max@0:   const uword n_slices = x.get_n_slices();
max@0:   
max@0:   arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "element-wise division");
max@0:   
max@0:   typedef typename T1::elem_type eT;
max@0:   
max@0:   eT* out_mem = out.memptr();
max@0:   
max@0:   const bool prefer_at_accessor = (ProxyCube<T1>::prefer_at_accessor || ProxyCube<T2>::prefer_at_accessor);
max@0:   
max@0:   if(prefer_at_accessor == false)
max@0:     {
max@0:     const uword n_elem  = out.n_elem;
max@0:     
max@0:     typename ProxyCube<T1>::ea_type P1 = x.P1.get_ea();
max@0:     typename ProxyCube<T2>::ea_type P2 = x.P2.get_ea();
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_1(/=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_1(/=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_1(/=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_1(/=, *); }
max@0:     }
max@0:   else
max@0:     {
max@0:     const ProxyCube<T1>& P1 = x.P1;
max@0:     const ProxyCube<T2>& P2 = x.P2;
max@0:     
max@0:          if(is_same_type<eglue_type, eglue_plus >::value == true) { arma_applier_3(/=, +); }
max@0:     else if(is_same_type<eglue_type, eglue_minus>::value == true) { arma_applier_3(/=, -); }
max@0:     else if(is_same_type<eglue_type, eglue_div  >::value == true) { arma_applier_3(/=, /); }
max@0:     else if(is_same_type<eglue_type, eglue_schur>::value == true) { arma_applier_3(/=, *); }
max@0:     }
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: #undef arma_applier_1
max@0: #undef arma_applier_2
max@0: #undef arma_applier_3
max@0: 
max@0: 
max@0: 
max@0: //! @}