CommaInitializer.h
Go to the documentation of this file.
00001 // This file is part of Eigen, a lightweight C++ template library
00002 // for linear algebra.
00003 //
00004 // Copyright (C) 2008 Gael Guennebaud <gael.guennebaud@inria.fr>
00005 // Copyright (C) 2006-2008 Benoit Jacob <jacob.benoit.1@gmail.com>
00006 //
00007 // This Source Code Form is subject to the terms of the Mozilla
00008 // Public License v. 2.0. If a copy of the MPL was not distributed
00009 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
00010 
00011 #ifndef EIGEN_COMMAINITIALIZER_H
00012 #define EIGEN_COMMAINITIALIZER_H
00013 
00014 namespace Eigen { 
00015 
00027 template<typename XprType>
00028 struct CommaInitializer
00029 {
00030   typedef typename XprType::Scalar Scalar;
00031   typedef typename XprType::Index Index;
00032 
00033   inline CommaInitializer(XprType& xpr, const Scalar& s)
00034     : m_xpr(xpr), m_row(0), m_col(1), m_currentBlockRows(1)
00035   {
00036     m_xpr.coeffRef(0,0) = s;
00037   }
00038 
00039   template<typename OtherDerived>
00040   inline CommaInitializer(XprType& xpr, const DenseBase<OtherDerived>& other)
00041     : m_xpr(xpr), m_row(0), m_col(other.cols()), m_currentBlockRows(other.rows())
00042   {
00043     m_xpr.block(0, 0, other.rows(), other.cols()) = other;
00044   }
00045 
00046   /* inserts a scalar value in the target matrix */
00047   CommaInitializer& operator,(const Scalar& s)
00048   {
00049     if (m_col==m_xpr.cols())
00050     {
00051       m_row+=m_currentBlockRows;
00052       m_col = 0;
00053       m_currentBlockRows = 1;
00054       eigen_assert(m_row<m_xpr.rows()
00055         && "Too many rows passed to comma initializer (operator<<)");
00056     }
00057     eigen_assert(m_col<m_xpr.cols()
00058       && "Too many coefficients passed to comma initializer (operator<<)");
00059     eigen_assert(m_currentBlockRows==1);
00060     m_xpr.coeffRef(m_row, m_col++) = s;
00061     return *this;
00062   }
00063 
00064   /* inserts a matrix expression in the target matrix */
00065   template<typename OtherDerived>
00066   CommaInitializer& operator,(const DenseBase<OtherDerived>& other)
00067   {
00068     if(other.cols()==0 || other.rows()==0)
00069       return *this;
00070     if (m_col==m_xpr.cols())
00071     {
00072       m_row+=m_currentBlockRows;
00073       m_col = 0;
00074       m_currentBlockRows = other.rows();
00075       eigen_assert(m_row+m_currentBlockRows<=m_xpr.rows()
00076         && "Too many rows passed to comma initializer (operator<<)");
00077     }
00078     eigen_assert(m_col<m_xpr.cols()
00079       && "Too many coefficients passed to comma initializer (operator<<)");
00080     eigen_assert(m_currentBlockRows==other.rows());
00081     if (OtherDerived::SizeAtCompileTime != Dynamic)
00082       m_xpr.template block<OtherDerived::RowsAtCompileTime != Dynamic ? OtherDerived::RowsAtCompileTime : 1,
00083                               OtherDerived::ColsAtCompileTime != Dynamic ? OtherDerived::ColsAtCompileTime : 1>
00084                     (m_row, m_col) = other;
00085     else
00086       m_xpr.block(m_row, m_col, other.rows(), other.cols()) = other;
00087     m_col += other.cols();
00088     return *this;
00089   }
00090 
00091   inline ~CommaInitializer()
00092   {
00093     eigen_assert((m_row+m_currentBlockRows) == m_xpr.rows()
00094          && m_col == m_xpr.cols()
00095          && "Too few coefficients passed to comma initializer (operator<<)");
00096   }
00097 
00105   inline XprType& finished() { return m_xpr; }
00106 
00107   XprType& m_xpr;   // target expression
00108   Index m_row;              // current row id
00109   Index m_col;              // current col id
00110   Index m_currentBlockRows; // current block height
00111 };
00112 
00124 template<typename Derived>
00125 inline CommaInitializer<Derived> DenseBase<Derived>::operator<< (const Scalar& s)
00126 {
00127   return CommaInitializer<Derived>(*static_cast<Derived*>(this), s);
00128 }
00129 
00131 template<typename Derived>
00132 template<typename OtherDerived>
00133 inline CommaInitializer<Derived>
00134 DenseBase<Derived>::operator<<(const DenseBase<OtherDerived>& other)
00135 {
00136   return CommaInitializer<Derived>(*static_cast<Derived *>(this), other);
00137 }
00138 
00139 } // end namespace Eigen
00140 
00141 #endif // EIGEN_COMMAINITIALIZER_H


win_eigen
Author(s): Daniel Stonier
autogenerated on Wed Sep 16 2015 07:10:26