max@0: // Copyright (C) 2008-2010 NICTA (www.nicta.com.au)
max@0: // Copyright (C) 2008-2010 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 op_dotext
max@0: //! @{
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: inline
max@0: eT
max@0: op_dotext::direct_rowvec_mat_colvec
max@0:   (
max@0:   const eT*      A_mem,
max@0:   const Mat<eT>& B,
max@0:   const eT*      C_mem
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword cost_AB = B.n_cols;
max@0:   const uword cost_BC = B.n_rows;
max@0:   
max@0:   if(cost_AB <= cost_BC)
max@0:     {
max@0:     podarray<eT> tmp(B.n_cols);
max@0:     
max@0:     for(uword col=0; col<B.n_cols; ++col)
max@0:       {
max@0:       const eT* B_coldata = B.colptr(col);
max@0:       
max@0:       eT val = eT(0);
max@0:       for(uword i=0; i<B.n_rows; ++i)
max@0:         {
max@0:         val += A_mem[i] * B_coldata[i];
max@0:         }
max@0:         
max@0:       tmp[col] = val;
max@0:       }
max@0:     
max@0:     return op_dot::direct_dot(B.n_cols, tmp.mem, C_mem);
max@0:     }
max@0:   else
max@0:     {
max@0:     podarray<eT> tmp(B.n_rows);
max@0:     
max@0:     for(uword row=0; row<B.n_rows; ++row)
max@0:       {
max@0:       eT val = eT(0);
max@0:       for(uword col=0; col<B.n_cols; ++col)
max@0:         {
max@0:         val += B.at(row,col) * C_mem[col];
max@0:         }
max@0:       
max@0:       tmp[row] = val;
max@0:       }
max@0:     
max@0:     return op_dot::direct_dot(B.n_rows, A_mem, tmp.mem);
max@0:     }
max@0:   
max@0:   
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: inline
max@0: eT
max@0: op_dotext::direct_rowvec_transmat_colvec
max@0:   (
max@0:   const eT*      A_mem,
max@0:   const Mat<eT>& B,
max@0:   const eT*      C_mem
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const uword cost_AB = B.n_rows;
max@0:   const uword cost_BC = B.n_cols;
max@0:   
max@0:   if(cost_AB <= cost_BC)
max@0:     {
max@0:     podarray<eT> tmp(B.n_rows);
max@0:     
max@0:     for(uword row=0; row<B.n_rows; ++row)
max@0:       {
max@0:       eT val = eT(0);
max@0:       
max@0:       for(uword i=0; i<B.n_cols; ++i)
max@0:         {
max@0:         val += A_mem[i] * B.at(row,i);
max@0:         }
max@0:         
max@0:       tmp[row] = val;
max@0:       }
max@0:     
max@0:     return op_dot::direct_dot(B.n_rows, tmp.mem, C_mem);
max@0:     }
max@0:   else
max@0:     {
max@0:     podarray<eT> tmp(B.n_cols);
max@0:     
max@0:     for(uword col=0; col<B.n_cols; ++col)
max@0:       {
max@0:       const eT* B_coldata = B.colptr(col);
max@0:       
max@0:       eT val = eT(0);
max@0:       
max@0:       for(uword i=0; i<B.n_rows; ++i)
max@0:         {
max@0:         val += B_coldata[i] * C_mem[i];
max@0:         }
max@0:       
max@0:       tmp[col] = val;
max@0:       }
max@0:     
max@0:     return op_dot::direct_dot(B.n_cols, A_mem, tmp.mem);
max@0:     }
max@0:   
max@0:   
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: inline
max@0: eT
max@0: op_dotext::direct_rowvec_diagmat_colvec
max@0:   (
max@0:   const eT*      A_mem,
max@0:   const Mat<eT>& B,
max@0:   const eT*      C_mem
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   eT val = eT(0);
max@0: 
max@0:   for(uword i=0; i<B.n_rows; ++i)
max@0:     {
max@0:     val += A_mem[i] * B.at(i,i) * C_mem[i];
max@0:     }
max@0: 
max@0:   return val;
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: inline
max@0: eT
max@0: op_dotext::direct_rowvec_invdiagmat_colvec
max@0:   (
max@0:   const eT*      A_mem,
max@0:   const Mat<eT>& B,
max@0:   const eT*      C_mem
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   eT val = eT(0);
max@0: 
max@0:   for(uword i=0; i<B.n_rows; ++i)
max@0:     {
max@0:     val += (A_mem[i] * C_mem[i]) / B.at(i,i);
max@0:     }
max@0: 
max@0:   return val;
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: template<typename eT>
max@0: inline
max@0: eT
max@0: op_dotext::direct_rowvec_invdiagvec_colvec
max@0:   (
max@0:   const eT*      A_mem,
max@0:   const Mat<eT>& B,
max@0:   const eT*      C_mem
max@0:   )
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   const eT* B_mem = B.mem;
max@0:   
max@0:   eT val = eT(0);
max@0: 
max@0:   for(uword i=0; i<B.n_elem; ++i)
max@0:     {
max@0:     val += (A_mem[i] * C_mem[i]) / B_mem[i];
max@0:     }
max@0: 
max@0:   return val;
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: //! @}