operator_cube_schur.hpp
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025 template<typename T1, typename T2>
00026 arma_inline
00027 const eGlueCube<T1, T2, eglue_cube_schur>
00028 operator%
00029 (
00030 const BaseCube<typename T1::elem_type,T1>& X,
00031 const BaseCube<typename T1::elem_type,T2>& Y
00032 )
00033 {
00034 arma_extra_debug_sigprint();
00035
00036 return eGlueCube<T1, T2, eglue_cube_schur>(X.get_ref(), Y.get_ref());
00037 }
00038
00039
00040
00041
00042 template<typename T1, typename T2>
00043 arma_inline
00044 Cube<typename promote_type<typename T1::elem_type, typename T2::elem_type>::result>
00045 operator%
00046 (
00047 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T1_result, T1>& X,
00048 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T2_result, T2>& Y
00049 )
00050 {
00051 arma_extra_debug_sigprint();
00052
00053 typedef typename T1::elem_type eT1;
00054 typedef typename T2::elem_type eT2;
00055
00056 typedef typename promote_type<eT1,eT2>::result out_eT;
00057
00058 promote_type<eT1,eT2>::check();
00059
00060 const ProxyCube<T1> A(X.get_ref());
00061 const ProxyCube<T2> B(Y.get_ref());
00062
00063 arma_debug_assert_same_size(A, B, "element-wise cube multiplication");
00064
00065 Cube<out_eT> out(A.n_rows, A.n_cols, A.n_slices);
00066
00067 out_eT* out_mem = out.memptr();
00068 const u32 n_elem = out.n_elem;
00069
00070 for(u32 i=0; i<n_elem; ++i)
00071 {
00072 out_mem[i] = upgrade_val<eT1,eT2>::apply(A[i]) * upgrade_val<eT1,eT2>::apply(B[i]);
00073 }
00074
00075 return out;
00076 }
00077
00078
00079
00080