Classes | Functions

Glue_mixed

//! More...

Classes

class  glue_mixed_times
class  glue_mixed_plus
class  glue_mixed_minus
class  glue_mixed_div
class  glue_mixed_schur

Functions

template<typename T1 , typename T2 >
static void glue_mixed_times::apply (Mat< typename eT_promoter< T1, T2 >::eT > &out, const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_times > &X)
 matrix multiplication with different element types
template<typename T1 , typename T2 >
static void glue_mixed_plus::apply (Mat< typename eT_promoter< T1, T2 >::eT > &out, const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_plus > &X)
 matrix addition with different element types
template<typename T1 , typename T2 >
static void glue_mixed_minus::apply (Mat< typename eT_promoter< T1, T2 >::eT > &out, const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_minus > &X)
 matrix subtraction with different element types
template<typename T1 , typename T2 >
static void glue_mixed_div::apply (Mat< typename eT_promoter< T1, T2 >::eT > &out, const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_div > &X)
 element-wise matrix division with different element types
template<typename T1 , typename T2 >
static void glue_mixed_schur::apply (Mat< typename eT_promoter< T1, T2 >::eT > &out, const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_schur > &X)
 element-wise matrix multiplication with different element types

Detailed Description

//!


Function Documentation

template<typename T1 , typename T2 >
void glue_mixed_times::apply ( Mat< typename eT_promoter< T1, T2 >::eT > &  out,
const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_times > &  X 
) [inline, static, inherited]

matrix multiplication with different element types

Definition at line 26 of file glue_mixed_meat.hpp.

References unwrap< T1 >::M, and Mat< eT >::n_cols.

  {
  arma_extra_debug_sigprint();
  
  typedef typename T1::elem_type eT1;
  typedef typename T2::elem_type eT2;
  
  // TODO: extend the unwrap_check framework to handle mixed matrix types
  
  const unwrap<T1> tmp1(X.A);
  const unwrap<T2> tmp2(X.B);
  
  const Mat<eT1>& A = tmp1.M;
  const Mat<eT2>& B = tmp2.M;
  
  const bool A_is_alias = ( ((void *)&out) == ((void *)&A) );
  const bool B_is_alias = ( ((void *)&out) == ((void *)&B) );
  
  const Mat<eT1>* AA_ptr = A_is_alias ? new Mat<eT1>(A) : 0;
  const Mat<eT2>* BB_ptr = B_is_alias ? new Mat<eT2>(B) : 0;
  
  const Mat<eT1>& AA = A_is_alias ? *AA_ptr : A;
  const Mat<eT2>& BB = B_is_alias ? *BB_ptr : B;
  
  arma_debug_assert_mul_size(AA, BB, "matrix multiplication");
  
  out.set_size(AA.n_rows, BB.n_cols);
  
  gemm_mixed<>::apply(out, AA, BB);
  
  if(A_is_alias == true)
    {
    delete AA_ptr;
    }
  
  if(B_is_alias == true)
    {
    delete BB_ptr;
    }
  }

template<typename T1 , typename T2 >
void glue_mixed_plus::apply ( Mat< typename eT_promoter< T1, T2 >::eT > &  out,
const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_plus > &  X 
) [inline, static, inherited]

matrix addition with different element types

Definition at line 73 of file glue_mixed_meat.hpp.

  {
  arma_extra_debug_sigprint();
  
  typedef typename T1::elem_type eT1;
  typedef typename T2::elem_type eT2;
  
  typedef typename promote_type<eT1,eT2>::result out_eT;
  
  promote_type<eT1,eT2>::check();
  
  const Proxy<T1> A(X.A);
  const Proxy<T2> B(X.B);
  
  arma_debug_assert_same_size(A, B, "matrix addition");
  
  out.set_size(A.n_rows, A.n_cols);
  
        out_eT* out_mem = out.memptr();
  const u32     n_elem  = out.n_elem;
  
  for(u32 i=0; i<n_elem; ++i)
    {
    out_mem[i] = upgrade_val<eT1,eT2>::apply(A[i]) + upgrade_val<eT1,eT2>::apply(B[i]);
    }
  }

template<typename T1 , typename T2 >
void glue_mixed_minus::apply ( Mat< typename eT_promoter< T1, T2 >::eT > &  out,
const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_minus > &  X 
) [inline, static, inherited]

matrix subtraction with different element types

Definition at line 106 of file glue_mixed_meat.hpp.

  {
  arma_extra_debug_sigprint();
  
  typedef typename T1::elem_type eT1;
  typedef typename T2::elem_type eT2;
  
  typedef typename promote_type<eT1,eT2>::result out_eT;
  
  promote_type<eT1,eT2>::check();
  
  const Proxy<T1> A(X.A);
  const Proxy<T2> B(X.B);
  
  arma_debug_assert_same_size(A, B, "matrix subtraction");
  
  out.set_size(A.n_rows, A.n_cols);
  
        out_eT* out_mem = out.memptr();
  const u32     n_elem  = out.n_elem;
  
  for(u32 i=0; i<n_elem; ++i)
    {
    out_mem[i] = upgrade_val<eT1,eT2>::apply(A[i]) - upgrade_val<eT1,eT2>::apply(B[i]);
    }
  }

template<typename T1 , typename T2 >
void glue_mixed_div::apply ( Mat< typename eT_promoter< T1, T2 >::eT > &  out,
const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_div > &  X 
) [inline, static, inherited]

element-wise matrix division with different element types

Definition at line 139 of file glue_mixed_meat.hpp.

  {
  arma_extra_debug_sigprint();
  
  typedef typename T1::elem_type eT1;
  typedef typename T2::elem_type eT2;
  
  typedef typename promote_type<eT1,eT2>::result out_eT;
  
  promote_type<eT1,eT2>::check();
  
  const Proxy<T1> A(X.A);
  const Proxy<T2> B(X.B);
  
  arma_debug_assert_same_size(A, B, "element-wise matrix division");
  
  out.set_size(A.n_rows, A.n_cols);
  
        out_eT* out_mem = out.memptr();
  const u32     n_elem  = out.n_elem;
  
  for(u32 i=0; i<n_elem; ++i)
    {
    out_mem[i] = upgrade_val<eT1,eT2>::apply(A[i]) / upgrade_val<eT1,eT2>::apply(B[i]);
    }
  }

template<typename T1 , typename T2 >
void glue_mixed_schur::apply ( Mat< typename eT_promoter< T1, T2 >::eT > &  out,
const mtGlue< typename eT_promoter< T1, T2 >::eT, T1, T2, glue_mixed_schur > &  X 
) [inline, static, inherited]

element-wise matrix multiplication with different element types

Definition at line 172 of file glue_mixed_meat.hpp.

  {
  arma_extra_debug_sigprint();
  
  typedef typename T1::elem_type eT1;
  typedef typename T2::elem_type eT2;
  
  typedef typename promote_type<eT1,eT2>::result out_eT;
  
  promote_type<eT1,eT2>::check();
  
  const Proxy<T1> A(X.A);
  const Proxy<T2> B(X.B);
  
  arma_debug_assert_same_size(A, B, "element-wise matrix multiplication");
  
  out.set_size(A.n_rows, A.n_cols);
  
        out_eT* out_mem = out.memptr();
  const u32     n_elem  = out.n_elem;
  
  for(u32 i=0; i<n_elem; ++i)
    {
    out_mem[i] = upgrade_val<eT1,eT2>::apply(A[i]) * upgrade_val<eT1,eT2>::apply(B[i]);
    }
  }