max@0: // Copyright (C) 2009-2012 NICTA (www.nicta.com.au)
max@0: // Copyright (C) 2009-2012 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 glue_relational
max@0: //! @{
max@0: 
max@0: 
max@0: 
max@0: #undef operator_rel
max@0: #undef operator_str
max@0: 
max@0: #undef arma_applier_mat
max@0: #undef arma_applier_cube
max@0: 
max@0: 
max@0: #define arma_applier_mat(operator_rel, operator_str) \
max@0:   {\
max@0:   const Proxy<T1> P1(X.A);\
max@0:   const Proxy<T2> P2(X.B);\
max@0:   \
max@0:   arma_debug_assert_same_size(P1, P2, operator_str);\
max@0:   \
max@0:   const bool bad_alias = (Proxy<T1>::has_subview && P1.is_alias(out)) || (Proxy<T2>::has_subview && P2.is_alias(out));\
max@0:   \
max@0:   if(bad_alias == false)\
max@0:     {\
max@0:     \
max@0:     const uword n_rows = P1.get_n_rows();\
max@0:     const uword n_cols = P1.get_n_cols();\
max@0:     \
max@0:     out.set_size(n_rows, n_cols);\
max@0:     \
max@0:     uword* 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:       typename Proxy<T1>::ea_type A = P1.get_ea();\
max@0:       typename Proxy<T2>::ea_type B = P2.get_ea();\
max@0:       \
max@0:       const uword n_elem = out.n_elem;\
max@0:       \
max@0:       for(uword i=0; i<n_elem; ++i)\
max@0:         {\
max@0:         out_mem[i] = (A[i] operator_rel B[i]) ? uword(1) : uword(0);\
max@0:         }\
max@0:       }\
max@0:     else\
max@0:       {\
max@0:       uword count = 0;\
max@0:       \
max@0:       for(uword col=0; col<n_cols; ++col)\
max@0:       for(uword row=0; row<n_rows; ++row, ++count)\
max@0:         {\
max@0:         out_mem[count] = (P1.at(row,col) operator_rel P2.at(row,col)) ? uword(1) : uword(0);\
max@0:         }\
max@0:       }\
max@0:     }\
max@0:   else\
max@0:     {\
max@0:     const unwrap<typename Proxy<T1>::stored_type> tmp1(P1.Q);\
max@0:     const unwrap<typename Proxy<T2>::stored_type> tmp2(P2.Q);\
max@0:     \
max@0:     out = (tmp1.M) operator_rel (tmp2.M);\
max@0:     }\
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: 
max@0: #define arma_applier_cube(operator_rel, operator_str) \
max@0:   {\
max@0:   const ProxyCube<T1> P1(X.A);\
max@0:   const ProxyCube<T2> P2(X.B);\
max@0:   \
max@0:   arma_debug_assert_same_size(P1, P2, operator_str);\
max@0:   \
max@0:   const bool bad_alias = (ProxyCube<T1>::has_subview && P1.is_alias(out)) || (ProxyCube<T2>::has_subview && P2.is_alias(out));\
max@0:   \
max@0:   if(bad_alias == false)\
max@0:     {\
max@0:     \
max@0:     const uword n_rows   = P1.get_n_rows();\
max@0:     const uword n_cols   = P1.get_n_cols();\
max@0:     const uword n_slices = P1.get_n_slices();\
max@0:     \
max@0:     out.set_size(n_rows, n_cols, n_slices);\
max@0:     \
max@0:     uword* 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:       typename ProxyCube<T1>::ea_type A = P1.get_ea();\
max@0:       typename ProxyCube<T2>::ea_type B = P2.get_ea();\
max@0:       \
max@0:       const uword n_elem = out.n_elem;\
max@0:       \
max@0:       for(uword i=0; i<n_elem; ++i)\
max@0:         {\
max@0:         out_mem[i] = (A[i] operator_rel B[i]) ? uword(1) : uword(0);\
max@0:         }\
max@0:       }\
max@0:     else\
max@0:       {\
max@0:       uword count = 0;\
max@0:       \
max@0:       for(uword slice = 0; slice < n_slices; ++slice)\
max@0:       for(uword col   = 0; col   < n_cols;   ++col)\
max@0:       for(uword row   = 0; row   < n_rows;   ++row, ++count)\
max@0:         {\
max@0:         out_mem[count] = (P1.at(row,col,slice) operator_rel P2.at(row,col,slice)) ? uword(1) : uword(0);\
max@0:         }\
max@0:       }\
max@0:     }\
max@0:   else\
max@0:     {\
max@0:     const unwrap_cube<typename ProxyCube<T1>::stored_type> tmp1(P1.Q);\
max@0:     const unwrap_cube<typename ProxyCube<T2>::stored_type> tmp2(P2.Q);\
max@0:     \
max@0:     out = (tmp1.M) operator_rel (tmp2.M);\
max@0:     }\
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_lt::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_lt>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(<, "operator<");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_gt::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_gt>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(>, "operator>");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_lteq::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_lteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(<=, "operator<=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_gteq::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_gteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(>=, "operator>=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_eq::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_eq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(==, "operator==");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_noteq::apply
max@0:   (
max@0:         Mat   <uword>& out,
max@0:   const mtGlue<uword, T1, T2, glue_rel_noteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_mat(!=, "operator!=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: //
max@0: //
max@0: //
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_lt::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_lt>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(<, "operator<");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_gt::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_gt>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(>, "operator>");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_lteq::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_lteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(<=, "operator<=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_gteq::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_gteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(>=, "operator>=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_eq::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_eq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(==, "operator==");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_rel_noteq::apply
max@0:   (
max@0:         Cube      <uword>& out,
max@0:   const mtGlueCube<uword, T1, T2, glue_rel_noteq>& X
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   arma_applier_cube(!=, "operator!=");
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: #undef arma_applier_mat
max@0: #undef arma_applier_cube
max@0: 
max@0: 
max@0: 
max@0: //! @}