+/* -*- C++ -*- ------------------------------------------------------------
+
+Copyright (c) 2007 Jesse Anders and Demian Nave http://cmldev.net/
+
+The Configurable Math Library (CML) is distributed under the terms of the
+Boost Software License, v1.0 (see cml/LICENSE for details).
+
+ *-----------------------------------------------------------------------*/
+/** @file
+ * @brief Defines the various combinations of matrix expressions.
+ *
+ * Create unary and binary operators with macros. The available combinations
+ * are:
+ *
+ * Unary expressions:
+ *
+ * op Matrix -> Matrix
+ * op MatXpr -> Matrix
+ *
+ * Binary expressions:
+ *
+ * Matrix op Matrix -> Matrix
+ * MatXpr op Matrix -> MatXpr
+ * Matrix op MatXpr -> MatXpr
+ * MatXpr op MatXpr -> MatXpr
+ *
+ * Matrix op Scalar -> Matrix
+ * Scalar op Matrix -> Matrix
+ * MatXpr op Scalar -> MatXpr
+ * Scalar op MatXpr -> MatXpr
+ *
+ * All of the generator functions compress the expression tree by hoisting
+ * subexpressions into the containing expression. This has the effect of
+ * forcing only the root node of the expression tree to be a MatrixXpr.
+ * Every other node is a Unary or BinaryMatrixOp.
+ */
+#ifndef matop_macros_h
+#define matop_macros_h
+
+/** Declare a unary operator taking a matrix operand. */
+#define CML_MAT_UNIOP(_op_, _OpT_) \
+template<typename E, class AT, typename BO, typename L> \
+inline et::MatrixXpr< \
+ et::UnaryMatrixOp< matrix<E,AT,BO,L>, _OpT_ <E> > \
+> \
+ \
+_op_ (const matrix<E,AT,BO,L>& arg) \
+{ \
+ typedef et::UnaryMatrixOp< \
+ matrix<E,AT,BO,L>, _OpT_ <E> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(arg)); \
+}
+
+/** Declare a unary operator taking a et::MatrixXpr operand. */
+#define CML_MATXPR_UNIOP(_op_, _OpT_) \
+template<class XprT> \
+inline et::MatrixXpr< \
+ et::UnaryMatrixOp<XprT, _OpT_<typename XprT::value_type> > \
+> \
+ \
+_op_ (MATXPR_ARG_TYPE arg) \
+{ \
+ typedef et::UnaryMatrixOp< \
+ XprT, _OpT_<typename XprT::value_type> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(arg.expression())); \
+}
+
+/** Declare an operator taking two matrix operands. */
+#define CML_MAT_MAT_BINOP(_op_, _OpT_) \
+template<typename E1, class AT1, typename L1, \
+ typename E2, class AT2, typename L2, typename BO> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ matrix<E1,AT1,BO,L2>, matrix<E2,AT2,BO,L2>, _OpT_<E1,E2> > \
+> \
+ \
+_op_ ( \
+ const matrix<E1,AT1,BO,L1>& left, \
+ const matrix<E2,AT2,BO,L2>& right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ matrix<E1,AT1,BO,L1>, matrix<E2,AT2,BO,L2>, _OpT_<E1,E2> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left,right)); \
+}
+
+/** Declare an operator taking a matrix and a et::MatrixXpr. */
+#define CML_MAT_MATXPR_BINOP(_op_, _OpT_) \
+template<typename E, class AT, typename BO, typename L, class XprT> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ matrix<E,AT,BO,L>, XprT, _OpT_ <E, typename XprT::value_type> \
+ > \
+> \
+ \
+_op_ ( \
+ const matrix<E,AT,BO,L>& left, \
+ MATXPR_ARG_TYPE right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ matrix<E,AT,BO,L>, XprT, \
+ _OpT_ <E, typename XprT::value_type> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left,right.expression())); \
+}
+
+/** Declare an operator taking a et::MatrixXpr and a matrix. */
+#define CML_MATXPR_MAT_BINOP(_op_, _OpT_) \
+template<class XprT, typename E, class AT, typename BO, typename L> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ XprT, matrix<E,AT,BO,L>, _OpT_ <typename XprT::value_type, E> \
+ > \
+> \
+ \
+_op_ ( \
+ MATXPR_ARG_TYPE left, \
+ const matrix<E,AT,BO,L>& right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ XprT, matrix<E,AT,BO,L>, \
+ _OpT_ <typename XprT::value_type, E> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left.expression(),right)); \
+}
+
+/** Declare an operator taking two et::MatrixXpr operands. */
+#define CML_MATXPR_MATXPR_BINOP(_op_, _OpT_) \
+template<class XprT1, class XprT2> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ XprT1, XprT2, \
+ _OpT_ < \
+ typename XprT1::value_type, \
+ typename XprT2::value_type \
+ > \
+ > \
+> \
+ \
+_op_ ( \
+ MATXPR_ARG_TYPE_N(1) left, \
+ MATXPR_ARG_TYPE_N(2) right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ XprT1, XprT2, \
+ _OpT_ < \
+ typename XprT1::value_type, \
+ typename XprT2::value_type> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>( \
+ ExprT(left.expression(),right.expression())); \
+}
+
+
+/** Declare an operator taking a matrix and a scalar. */
+#define CML_MAT_SCALAR_BINOP(_op_, _OpT_) \
+template<typename E, class AT, typename BO, typename L, typename ScalarT>\
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ matrix<E,AT,BO,L>, ScalarT, _OpT_ <E,ScalarT> \
+ > \
+> \
+ \
+_op_ ( \
+ const matrix<E,AT,BO,L>& left, \
+ SCALAR_ARG_TYPE right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ matrix<E,AT,BO,L>, ScalarT, _OpT_ <E,ScalarT > \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left,right)); \
+}
+
+/** Declare an operator taking a scalar and a matrix. */
+#define CML_SCALAR_MAT_BINOP(_op_, _OpT_) \
+template<typename ScalarT, typename E, class AT, typename BO, typename L>\
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ ScalarT, matrix<E,AT,BO,L>, _OpT_ <ScalarT,E> \
+ > \
+> \
+ \
+_op_ ( \
+ SCALAR_ARG_TYPE left, \
+ const matrix<E,AT,BO,L>& right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ ScalarT, matrix<E,AT,BO,L>, _OpT_<ScalarT,E> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left,right)); \
+}
+
+/** Declare an operator taking a et::MatrixXpr and a scalar. */
+#define CML_MATXPR_SCALAR_BINOP(_op_, _OpT_) \
+template<class XprT, typename ScalarT> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ XprT, ScalarT, _OpT_ <typename XprT::value_type, ScalarT> \
+ > \
+> \
+ \
+_op_ ( \
+ MATXPR_ARG_TYPE left, \
+ SCALAR_ARG_TYPE right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ XprT, ScalarT, _OpT_ <typename XprT::value_type,ScalarT> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left.expression(),right)); \
+}
+
+/** Declare an operator taking a scalar and a et::MatrixXpr. */
+#define CML_SCALAR_MATXPR_BINOP(_op_, _OpT_) \
+template<typename ScalarT, class XprT> \
+inline et::MatrixXpr< \
+ et::BinaryMatrixOp< \
+ ScalarT, XprT, _OpT_ <ScalarT, typename XprT::value_type> \
+ > \
+> \
+ \
+_op_ ( \
+ SCALAR_ARG_TYPE left, \
+ MATXPR_ARG_TYPE right) \
+{ \
+ typedef et::BinaryMatrixOp< \
+ ScalarT, XprT, _OpT_ <ScalarT, typename XprT::value_type> \
+ > ExprT; \
+ return et::MatrixXpr<ExprT>(ExprT(left,right.expression())); \
+}
+
+#endif
+
+// -------------------------------------------------------------------------
+// vim:ft=cpp