+++ /dev/null
-/* -*- 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
- *
- * @note GCC4 requires a copy constructor to elide---it won't elide a
- * compiler-generated copy constructor!
- */
-
-#ifndef matrix_class_ops_h
-#define matrix_class_ops_h
-
-#if defined(_MSC_VER) && _MSC_VER < 1400
-#pragma warning(disable:4003)
-// XXX Horrible hack to turn off warnings about "not enough actual params"
-// for the macros below.
-#endif
-
-/* This is to circumvent the problem of commas in a macro argument. It's
- * used to instantiate CML_ACCUMULATED_MATRIX_MULT:
- */
-#define TEMPLATED_MATRIX_MACRO matrix<E,AT,BO,L>
-
-/* XXX HACK!!! This is a hack to resize in the assign() functions only when
- * auto resizing is turned off.
- */
-#if !defined(CML_MATRIX_RESIZE_ON_ASSIGNMENT)
-#define _DO_MATRIX_SET_RESIZE(_R_,_C_) cml::et::detail::Resize(*this,_R_,_C_)
-#else
-#define _DO_MATRIX_SET_RESIZE(_R_,_C_)
-#endif
-
-/** Set a matrix from 2x2 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_ASSIGN_MAT_22 \
-matrix_type& \
-set( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11 \
- ) \
-{ \
- _DO_MATRIX_SET_RESIZE(2,2); \
- /* This is overkill, but simplifies size checking: */ \
- value_type v[2][2] = {{e00,e01},{e10,e11}}; \
- typedef et::OpAssign<Element,Element> OpT; \
- typedef const value_type element; \
- cml::matrix<element, external<2,2>, basis_orient, row_major> \
- src(&v[0][0]); \
- et::UnrollAssignment<OpT>(*this,src); \
- return *this; \
-}
-
-/** Create a matrix from 3x3 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_ASSIGN_MAT_33 \
-matrix_type& \
-set( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, ELEMENT_ARG_TYPE e02, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11, ELEMENT_ARG_TYPE e12, \
- ELEMENT_ARG_TYPE e20, ELEMENT_ARG_TYPE e21, ELEMENT_ARG_TYPE e22 \
- ) \
-{ \
- _DO_MATRIX_SET_RESIZE(3,3); \
- /* This is overkill, but simplifies size checking: */ \
- value_type v[3][3] = { \
- {e00,e01,e02}, \
- {e10,e11,e12}, \
- {e20,e21,e22} \
- }; \
- typedef et::OpAssign<Element,Element> OpT; \
- typedef const value_type element; \
- cml::matrix<element, external<3,3>, basis_orient, row_major> \
- src(&v[0][0]); \
- et::UnrollAssignment<OpT>(*this,src); \
- return *this; \
-}
-
-/** Create a matrix from 4x4 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_ASSIGN_MAT_44 \
-matrix_type& \
-set( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, \
- ELEMENT_ARG_TYPE e02, ELEMENT_ARG_TYPE e03, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11, \
- ELEMENT_ARG_TYPE e12, ELEMENT_ARG_TYPE e13, \
- ELEMENT_ARG_TYPE e20, ELEMENT_ARG_TYPE e21, \
- ELEMENT_ARG_TYPE e22, ELEMENT_ARG_TYPE e23, \
- ELEMENT_ARG_TYPE e30, ELEMENT_ARG_TYPE e31, \
- ELEMENT_ARG_TYPE e32, ELEMENT_ARG_TYPE e33 \
- ) \
-{ \
- _DO_MATRIX_SET_RESIZE(4,4); \
- /* This is overkill, but simplifies size checking: */ \
- value_type v[4][4] = { \
- {e00,e01,e02,e03}, \
- {e10,e11,e12,e13}, \
- {e20,e21,e22,e23}, \
- {e30,e31,e32,e33} \
- }; \
- typedef et::OpAssign<Element,Element> OpT; \
- typedef const value_type element; \
- cml::matrix<element, external<4,4>, basis_orient, row_major> \
- src(&v[0][0]); \
- et::UnrollAssignment<OpT>(*this,src); \
- return *this; \
-}
-
-/** Create a matrix from 2x2 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_CONSTRUCT_MAT_22 \
-matrix( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11 \
- ) \
-{ \
- set( \
- e00,e01, \
- e10,e11 \
- ); \
-}
-
-/** Create a matrix from 3x3 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_CONSTRUCT_MAT_33 \
-matrix( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, ELEMENT_ARG_TYPE e02, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11, ELEMENT_ARG_TYPE e12, \
- ELEMENT_ARG_TYPE e20, ELEMENT_ARG_TYPE e21, ELEMENT_ARG_TYPE e22 \
- ) \
-{ \
- set( \
- e00,e01,e02, \
- e10,e11,e12, \
- e20,e21,e22 \
- ); \
-}
-
-/** Create a matrix from 4x4 values.
- *
- * The layout assumed for the values is that of the matrix being assigned.
- */
-#define CML_CONSTRUCT_MAT_44 \
-matrix( \
- ELEMENT_ARG_TYPE e00, ELEMENT_ARG_TYPE e01, \
- ELEMENT_ARG_TYPE e02, ELEMENT_ARG_TYPE e03, \
- ELEMENT_ARG_TYPE e10, ELEMENT_ARG_TYPE e11, \
- ELEMENT_ARG_TYPE e12, ELEMENT_ARG_TYPE e13, \
- ELEMENT_ARG_TYPE e20, ELEMENT_ARG_TYPE e21, \
- ELEMENT_ARG_TYPE e22, ELEMENT_ARG_TYPE e23, \
- ELEMENT_ARG_TYPE e30, ELEMENT_ARG_TYPE e31, \
- ELEMENT_ARG_TYPE e32, ELEMENT_ARG_TYPE e33 \
- ) \
-{ \
- set( \
- e00,e01,e02,e03, \
- e10,e11,e12,e13, \
- e20,e21,e22,e23, \
- e30,e31,e32,e33 \
- ); \
-}
-
-/** Copy-construct a matrix from a fixed-size array of values. */
-#define CML_MAT_COPY_FROM_FIXED_ARRAY(_R_,_C_) \
-matrix(const value_type m[_R_][_C_]) { \
- typedef et::OpAssign<Element,Element> OpT; \
- cml::matrix<const value_type, external<_R_,_C_>, \
- basis_orient, row_major> src(&m[0][0]); \
- et::UnrollAssignment<OpT>(*this,src); \
-}
-
-/** Copy-construct a matrix from a runtime-sized array of values. */
-#define CML_MAT_COPY_FROM_ARRAY(_add_) \
-matrix(const value_type* const v, size_t R, size_t C) _add_ { \
- typedef et::OpAssign<Element,Element> OpT; \
- cml::matrix<value_type, external<>, basis_orient, \
- row_major > src(const_cast<value_type*>(v),R,C); \
- et::UnrollAssignment<OpT>(*this,src); \
-}
-
-/** Copy this matrix from another using the given elementwise op.
- *
- * @internal This is required for GCC4, since it won't elide the default
- * copy constructor.
- */
-#define CML_MAT_COPY_FROM_MATTYPE \
-matrix(const matrix_type& m) : array_type() { \
- typedef et::OpAssign <Element,Element> OpT; \
- et::UnrollAssignment<OpT>(*this,m); \
-}
-
-/** Copy this matrix from another using the given elementwise op.
- *
- * This allows copies from arbitrary matrix types.
- */
-#define CML_MAT_COPY_FROM_MAT \
-template<typename E, class AT, typename BO, typename L> \
-matrix(const TEMPLATED_MATRIX_MACRO& m) { \
- typedef et::OpAssign <Element,E> OpT; \
- et::UnrollAssignment<OpT>(*this,m); \
-}
-
-/** Declare a function to copy this matrix from a matrix expression. */
-#define CML_MAT_COPY_FROM_MATXPR \
-template<class XprT> \
-matrix(MATXPR_ARG_TYPE e) { \
- /* Verify that a promotion exists at compile time: */ \
- typedef typename et::MatrixPromote< \
- matrix_type, typename XprT::result_type>::type result_type; \
- typedef typename XprT::value_type src_value_type; \
- typedef et::OpAssign <Element,src_value_type> OpT; \
- et::UnrollAssignment<OpT>(*this,e); \
-}
-
-#if defined(CML_USE_GENERATED_MATRIX_ASSIGN_OP)
-#define CML_MAT_ASSIGN_FROM_MATTYPE
-#else
-/** Assign from the same matrix type.
- *
- * @param m the matrix to copy from.
- *
- * @note This is required for GCC4, otherwise it generates a slow
- * assignment operator using memcpy.
- *
- * @note ICC9/Linux-x86 seems to prefer its own assignment operator (need
- * to figure out why).
- */
-#define CML_MAT_ASSIGN_FROM_MATTYPE \
-matrix_type& operator=(const matrix_type& m) { \
- typedef et::OpAssign<Element,Element> OpT; \
- et::UnrollAssignment<OpT>(*this,m); \
- return *this; \
-}
-#endif
-
-
-/** Assign this matrix from another using the given elementwise op.
- *
- * This allows assignment from arbitrary matrix types.
- *
- * @param _op_ the operator (e.g. +=)
- * @param _op_name_ the op functor (e.g. et::OpAssign)
- */
-#define CML_MAT_ASSIGN_FROM_MAT(_op_, _op_name_) \
-template<typename E, class AT, typename BO, typename L> matrix_type& \
-operator _op_ (const TEMPLATED_MATRIX_MACRO& m) { \
- typedef _op_name_ <Element,E> OpT; \
- et::UnrollAssignment<OpT>(*this,m); \
- return *this; \
-}
-
-/** Declare a function to assign this matrix from a matrix expression.
- *
- * @param _op_ the operator (e.g. +=)
- * @param _op_name_ the op functor (e.g. et::OpAssign)
- */
-#define CML_MAT_ASSIGN_FROM_MATXPR(_op_, _op_name_) \
-template<class XprT> matrix_type& \
-operator _op_ (MATXPR_ARG_TYPE e) { \
- /* Verify that a promotion exists at compile time: */ \
- typedef typename et::MatrixPromote< \
- matrix_type, typename XprT::result_type>::type result_type; \
- typedef typename XprT::value_type src_value_type; \
- typedef _op_name_ <Element,src_value_type> OpT; \
- et::UnrollAssignment<OpT>(*this,e); \
- return *this; \
-}
-
-/** Declare a function to assign this matrix from a scalar.
- *
- * @param _op_ the operator (e.g. +=)
- * @param _op_name_ the op functor (e.g. et::OpAssign)
- *
- * @internal This shouldn't be used for ops, like +=, which aren't
- * defined in vector algebra.
- */
-#define CML_MAT_ASSIGN_FROM_SCALAR(_op_, _op_name_) \
-matrix_type& operator _op_ (ELEMENT_ARG_TYPE s) { \
- typedef _op_name_ <Element,value_type> OpT; \
- et::UnrollAssignment<OpT>(*this,s); \
- return *this; \
-}
-
-
-/** Accumulated matrix multiplication.
- *
- * @throws std::invalid_argument if the matrices are not square.
- */
-#define CML_ACCUMULATED_MATRIX_MULT(_arg_type_) \
-matrix_type& operator*=(_arg_type_ m) { \
- typedef typename et::MatrixPromote< \
- matrix_type, _arg_type_>::type result_type; \
- cml::et::CheckedSquare(*this, typename result_type::size_tag()); \
- return (*this = (*this)*m); \
-}
-
-
-/* These should only be used for testing: */
-#define CML_MATRIX_BRACE_OPERATORS \
-template<class Matrix> struct row_ref { \
- typedef typename Matrix::reference reference; \
- reference operator[](size_t col) { return m(row,col); } \
- Matrix& m; \
- size_t row; \
-}; \
- \
-template<class Matrix> struct const_row_ref { \
- typedef typename Matrix::const_reference const_reference; \
- const_reference operator[](size_t col) const { return m(row,col); } \
- const Matrix& m; \
- size_t row; \
-}; \
- \
-row_ref<matrix_type> operator[](size_t row) { \
- row_ref<matrix_type> ref = { *this, row }; return ref; \
-} \
- \
-const_row_ref<matrix_type> operator[](size_t row) const { \
- const_row_ref<matrix_type> ref = { *this, row }; return ref; \
-}
-
-#endif
-
-// -------------------------------------------------------------------------
-// vim:ft=cpp