mat_mat_plus_expr.hpp

矩阵运算源码最新版本

// Software License for MTL
// 
// Copyright (c) 2007 The Trustees of Indiana University. All rights reserved.
// Authors: Peter Gottschling and Andrew Lumsdaine
// 
// This file is part of the Matrix Template Library
// 
// See also license.mtl.txt in the distribution.

#ifndef MTL_MAT_MAT_PLUS_EXPR_INCLUDE
#define MTL_MAT_MAT_PLUS_EXPR_INCLUDE

#include <boost/numeric/mtl/matrix/mat_mat_op_expr.hpp>
#include <boost/numeric/mtl/operation/sfunctor.hpp>

namespace mtl { namespace matrix {

template <typename E1, typename E2>
struct mat_mat_plus_expr 
    : public mat_mat_op_expr< E1, E2, sfunctor::plus<typename E1::value_type, typename E2::value_type> >,
      public mat_expr< mat_mat_plus_expr<E1, E2> >
{
    typedef mat_mat_op_expr< E1, E2, sfunctor::plus<typename E1::value_type, typename E2::value_type> > op_base;
    typedef mat_expr< mat_mat_plus_expr<E1, E2> >                                                       crtp_base;
    typedef E1                                   first_argument_type ;
    typedef E2                                   second_argument_type ;
    
    mat_mat_plus_expr( E1 const& v1, E2 const& v2 )
	: op_base( v1, v2 ), crtp_base(*this), first(v1), second(v2)
    {}

    first_argument_type const&  first ;
    second_argument_type const& second ;
};

// Same as mat_mat_plus_expr for pair of dense matrix expressions
// Future versions will probably provide more efficient implementations for it
template <typename E1, typename E2>
struct dmat_dmat_plus_expr 
    : public mat_mat_op_expr< E1, E2, sfunctor::plus<typename E1::value_type, typename E2::value_type> >
{
    typedef mat_mat_op_expr< E1, E2, sfunctor::plus<typename E1::value_type, typename E2::value_type> > base;
    dmat_dmat_plus_expr( E1 const& v1, E2 const& v2 )
	: base( v1, v2 )
    {}
};
    


}} // Namespace mtl::matrix

#endif // MTL_MAT_MAT_PLUS_EXPR_INCLUDE