Jaesung Oh / Eigen

Eigne Matrix Class Library

Dependents:   MPC_current_control HydraulicControlBoard_SW AHRS Test_ekf ... more

src/Core/Block.h@0:13a5d365ba16, 2016-10-13 (annotated)

Committer:
ykuroda
Date:
Thu Oct 13 04:07:23 2016 +0000
Revision:
0:13a5d365ba16
First commint, Eigne Matrix Class Library

Who changed what in which revision?

UserRevisionLine numberNew contents of line
ykuroda 0:13a5d365ba16 1 // This file is part of Eigen, a lightweight C++ template library
ykuroda 0:13a5d365ba16 2 // for linear algebra.
ykuroda 0:13a5d365ba16 3 //
ykuroda 0:13a5d365ba16 4 // Copyright (C) 2008 Gael Guennebaud <gael.guennebaud@inria.fr>
ykuroda 0:13a5d365ba16 5 // Copyright (C) 2006-2010 Benoit Jacob <jacob.benoit.1@gmail.com>
ykuroda 0:13a5d365ba16 6 //
ykuroda 0:13a5d365ba16 7 // This Source Code Form is subject to the terms of the Mozilla
ykuroda 0:13a5d365ba16 8 // Public License v. 2.0. If a copy of the MPL was not distributed
ykuroda 0:13a5d365ba16 9 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
ykuroda 0:13a5d365ba16 10
ykuroda 0:13a5d365ba16 11 #ifndef EIGEN_BLOCK_H
ykuroda 0:13a5d365ba16 12 #define EIGEN_BLOCK_H
ykuroda 0:13a5d365ba16 13
ykuroda 0:13a5d365ba16 14 namespace Eigen {
ykuroda 0:13a5d365ba16 15
ykuroda 0:13a5d365ba16 16 /** \class Block
ykuroda 0:13a5d365ba16 17 * \ingroup Core_Module
ykuroda 0:13a5d365ba16 18 *
ykuroda 0:13a5d365ba16 19 * \brief Expression of a fixed-size or dynamic-size block
ykuroda 0:13a5d365ba16 20 *
ykuroda 0:13a5d365ba16 21 * \param XprType the type of the expression in which we are taking a block
ykuroda 0:13a5d365ba16 22 * \param BlockRows the number of rows of the block we are taking at compile time (optional)
ykuroda 0:13a5d365ba16 23 * \param BlockCols the number of columns of the block we are taking at compile time (optional)
ykuroda 0:13a5d365ba16 24 *
ykuroda 0:13a5d365ba16 25 * This class represents an expression of either a fixed-size or dynamic-size block. It is the return
ykuroda 0:13a5d365ba16 26 * type of DenseBase::block(Index,Index,Index,Index) and DenseBase::block<int,int>(Index,Index) and
ykuroda 0:13a5d365ba16 27 * most of the time this is the only way it is used.
ykuroda 0:13a5d365ba16 28 *
ykuroda 0:13a5d365ba16 29 * However, if you want to directly maniputate block expressions,
ykuroda 0:13a5d365ba16 30 * for instance if you want to write a function returning such an expression, you
ykuroda 0:13a5d365ba16 31 * will need to use this class.
ykuroda 0:13a5d365ba16 32 *
ykuroda 0:13a5d365ba16 33 * Here is an example illustrating the dynamic case:
ykuroda 0:13a5d365ba16 34 * \include class_Block.cpp
ykuroda 0:13a5d365ba16 35 * Output: \verbinclude class_Block.out
ykuroda 0:13a5d365ba16 36 *
ykuroda 0:13a5d365ba16 37 * \note Even though this expression has dynamic size, in the case where \a XprType
ykuroda 0:13a5d365ba16 38 * has fixed size, this expression inherits a fixed maximal size which means that evaluating
ykuroda 0:13a5d365ba16 39 * it does not cause a dynamic memory allocation.
ykuroda 0:13a5d365ba16 40 *
ykuroda 0:13a5d365ba16 41 * Here is an example illustrating the fixed-size case:
ykuroda 0:13a5d365ba16 42 * \include class_FixedBlock.cpp
ykuroda 0:13a5d365ba16 43 * Output: \verbinclude class_FixedBlock.out
ykuroda 0:13a5d365ba16 44 *
ykuroda 0:13a5d365ba16 45 * \sa DenseBase::block(Index,Index,Index,Index), DenseBase::block(Index,Index), class VectorBlock
ykuroda 0:13a5d365ba16 46 */
ykuroda 0:13a5d365ba16 47
ykuroda 0:13a5d365ba16 48 namespace internal {
ykuroda 0:13a5d365ba16 49 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel>
ykuroda 0:13a5d365ba16 50 struct traits<Block<XprType, BlockRows, BlockCols, InnerPanel> > : traits<XprType>
ykuroda 0:13a5d365ba16 51 {
ykuroda 0:13a5d365ba16 52 typedef typename traits<XprType>::Scalar Scalar;
ykuroda 0:13a5d365ba16 53 typedef typename traits<XprType>::StorageKind StorageKind;
ykuroda 0:13a5d365ba16 54 typedef typename traits<XprType>::XprKind XprKind;
ykuroda 0:13a5d365ba16 55 typedef typename nested<XprType>::type XprTypeNested;
ykuroda 0:13a5d365ba16 56 typedef typename remove_reference<XprTypeNested>::type _XprTypeNested;
ykuroda 0:13a5d365ba16 57 enum{
ykuroda 0:13a5d365ba16 58 MatrixRows = traits<XprType>::RowsAtCompileTime,
ykuroda 0:13a5d365ba16 59 MatrixCols = traits<XprType>::ColsAtCompileTime,
ykuroda 0:13a5d365ba16 60 RowsAtCompileTime = MatrixRows == 0 ? 0 : BlockRows,
ykuroda 0:13a5d365ba16 61 ColsAtCompileTime = MatrixCols == 0 ? 0 : BlockCols,
ykuroda 0:13a5d365ba16 62 MaxRowsAtCompileTime = BlockRows==0 ? 0
ykuroda 0:13a5d365ba16 63 : RowsAtCompileTime != Dynamic ? int(RowsAtCompileTime)
ykuroda 0:13a5d365ba16 64 : int(traits<XprType>::MaxRowsAtCompileTime),
ykuroda 0:13a5d365ba16 65 MaxColsAtCompileTime = BlockCols==0 ? 0
ykuroda 0:13a5d365ba16 66 : ColsAtCompileTime != Dynamic ? int(ColsAtCompileTime)
ykuroda 0:13a5d365ba16 67 : int(traits<XprType>::MaxColsAtCompileTime),
ykuroda 0:13a5d365ba16 68 XprTypeIsRowMajor = (int(traits<XprType>::Flags)&RowMajorBit) != 0,
ykuroda 0:13a5d365ba16 69 IsDense = is_same<StorageKind,Dense>::value,
ykuroda 0:13a5d365ba16 70 IsRowMajor = (IsDense&&MaxRowsAtCompileTime==1&&MaxColsAtCompileTime!=1) ? 1
ykuroda 0:13a5d365ba16 71 : (IsDense&&MaxColsAtCompileTime==1&&MaxRowsAtCompileTime!=1) ? 0
ykuroda 0:13a5d365ba16 72 : XprTypeIsRowMajor,
ykuroda 0:13a5d365ba16 73 HasSameStorageOrderAsXprType = (IsRowMajor == XprTypeIsRowMajor),
ykuroda 0:13a5d365ba16 74 InnerSize = IsRowMajor ? int(ColsAtCompileTime) : int(RowsAtCompileTime),
ykuroda 0:13a5d365ba16 75 InnerStrideAtCompileTime = HasSameStorageOrderAsXprType
ykuroda 0:13a5d365ba16 76 ? int(inner_stride_at_compile_time<XprType>::ret)
ykuroda 0:13a5d365ba16 77 : int(outer_stride_at_compile_time<XprType>::ret),
ykuroda 0:13a5d365ba16 78 OuterStrideAtCompileTime = HasSameStorageOrderAsXprType
ykuroda 0:13a5d365ba16 79 ? int(outer_stride_at_compile_time<XprType>::ret)
ykuroda 0:13a5d365ba16 80 : int(inner_stride_at_compile_time<XprType>::ret),
ykuroda 0:13a5d365ba16 81 MaskPacketAccessBit = (InnerSize == Dynamic || (InnerSize % packet_traits<Scalar>::size) == 0)
ykuroda 0:13a5d365ba16 82 && (InnerStrideAtCompileTime == 1)
ykuroda 0:13a5d365ba16 83 ? PacketAccessBit : 0,
ykuroda 0:13a5d365ba16 84 MaskAlignedBit = (InnerPanel && (OuterStrideAtCompileTime!=Dynamic) && (((OuterStrideAtCompileTime * int(sizeof(Scalar))) % 16) == 0)) ? AlignedBit : 0,
ykuroda 0:13a5d365ba16 85 FlagsLinearAccessBit = (RowsAtCompileTime == 1 || ColsAtCompileTime == 1 || (InnerPanel && (traits<XprType>::Flags&LinearAccessBit))) ? LinearAccessBit : 0,
ykuroda 0:13a5d365ba16 86 FlagsLvalueBit = is_lvalue<XprType>::value ? LvalueBit : 0,
ykuroda 0:13a5d365ba16 87 FlagsRowMajorBit = IsRowMajor ? RowMajorBit : 0,
ykuroda 0:13a5d365ba16 88 Flags0 = traits<XprType>::Flags & ( (HereditaryBits & ~RowMajorBit) |
ykuroda 0:13a5d365ba16 89 DirectAccessBit |
ykuroda 0:13a5d365ba16 90 MaskPacketAccessBit |
ykuroda 0:13a5d365ba16 91 MaskAlignedBit),
ykuroda 0:13a5d365ba16 92 Flags = Flags0 | FlagsLinearAccessBit | FlagsLvalueBit | FlagsRowMajorBit
ykuroda 0:13a5d365ba16 93 };
ykuroda 0:13a5d365ba16 94 };
ykuroda 0:13a5d365ba16 95
ykuroda 0:13a5d365ba16 96 template<typename XprType, int BlockRows=Dynamic, int BlockCols=Dynamic, bool InnerPanel = false,
ykuroda 0:13a5d365ba16 97 bool HasDirectAccess = internal::has_direct_access<XprType>::ret> class BlockImpl_dense;
ykuroda 0:13a5d365ba16 98
ykuroda 0:13a5d365ba16 99 } // end namespace internal
ykuroda 0:13a5d365ba16 100
ykuroda 0:13a5d365ba16 101 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel, typename StorageKind> class BlockImpl;
ykuroda 0:13a5d365ba16 102
ykuroda 0:13a5d365ba16 103 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel> class Block
ykuroda 0:13a5d365ba16 104 : public BlockImpl<XprType, BlockRows, BlockCols, InnerPanel, typename internal::traits<XprType>::StorageKind>
ykuroda 0:13a5d365ba16 105 {
ykuroda 0:13a5d365ba16 106 typedef BlockImpl<XprType, BlockRows, BlockCols, InnerPanel, typename internal::traits<XprType>::StorageKind> Impl;
ykuroda 0:13a5d365ba16 107 public:
ykuroda 0:13a5d365ba16 108 //typedef typename Impl::Base Base;
ykuroda 0:13a5d365ba16 109 typedef Impl Base;
ykuroda 0:13a5d365ba16 110 EIGEN_GENERIC_PUBLIC_INTERFACE(Block)
ykuroda 0:13a5d365ba16 111 EIGEN_INHERIT_ASSIGNMENT_OPERATORS(Block)
ykuroda 0:13a5d365ba16 112
ykuroda 0:13a5d365ba16 113 /** Column or Row constructor
ykuroda 0:13a5d365ba16 114 */
ykuroda 0:13a5d365ba16 115 inline Block(XprType& xpr, Index i) : Impl(xpr,i)
ykuroda 0:13a5d365ba16 116 {
ykuroda 0:13a5d365ba16 117 eigen_assert( (i>=0) && (
ykuroda 0:13a5d365ba16 118 ((BlockRows==1) && (BlockCols==XprType::ColsAtCompileTime) && i<xpr.rows())
ykuroda 0:13a5d365ba16 119 ||((BlockRows==XprType::RowsAtCompileTime) && (BlockCols==1) && i<xpr.cols())));
ykuroda 0:13a5d365ba16 120 }
ykuroda 0:13a5d365ba16 121
ykuroda 0:13a5d365ba16 122 /** Fixed-size constructor
ykuroda 0:13a5d365ba16 123 */
ykuroda 0:13a5d365ba16 124 inline Block(XprType& xpr, Index a_startRow, Index a_startCol)
ykuroda 0:13a5d365ba16 125 : Impl(xpr, a_startRow, a_startCol)
ykuroda 0:13a5d365ba16 126 {
ykuroda 0:13a5d365ba16 127 EIGEN_STATIC_ASSERT(RowsAtCompileTime!=Dynamic && ColsAtCompileTime!=Dynamic,THIS_METHOD_IS_ONLY_FOR_FIXED_SIZE)
ykuroda 0:13a5d365ba16 128 eigen_assert(a_startRow >= 0 && BlockRows >= 1 && a_startRow + BlockRows <= xpr.rows()
ykuroda 0:13a5d365ba16 129 && a_startCol >= 0 && BlockCols >= 1 && a_startCol + BlockCols <= xpr.cols());
ykuroda 0:13a5d365ba16 130 }
ykuroda 0:13a5d365ba16 131
ykuroda 0:13a5d365ba16 132 /** Dynamic-size constructor
ykuroda 0:13a5d365ba16 133 */
ykuroda 0:13a5d365ba16 134 inline Block(XprType& xpr,
ykuroda 0:13a5d365ba16 135 Index a_startRow, Index a_startCol,
ykuroda 0:13a5d365ba16 136 Index blockRows, Index blockCols)
ykuroda 0:13a5d365ba16 137 : Impl(xpr, a_startRow, a_startCol, blockRows, blockCols)
ykuroda 0:13a5d365ba16 138 {
ykuroda 0:13a5d365ba16 139 eigen_assert((RowsAtCompileTime==Dynamic || RowsAtCompileTime==blockRows)
ykuroda 0:13a5d365ba16 140 && (ColsAtCompileTime==Dynamic || ColsAtCompileTime==blockCols));
ykuroda 0:13a5d365ba16 141 eigen_assert(a_startRow >= 0 && blockRows >= 0 && a_startRow <= xpr.rows() - blockRows
ykuroda 0:13a5d365ba16 142 && a_startCol >= 0 && blockCols >= 0 && a_startCol <= xpr.cols() - blockCols);
ykuroda 0:13a5d365ba16 143 }
ykuroda 0:13a5d365ba16 144 };
ykuroda 0:13a5d365ba16 145
ykuroda 0:13a5d365ba16 146 // The generic default implementation for dense block simplu forward to the internal::BlockImpl_dense
ykuroda 0:13a5d365ba16 147 // that must be specialized for direct and non-direct access...
ykuroda 0:13a5d365ba16 148 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel>
ykuroda 0:13a5d365ba16 149 class BlockImpl<XprType, BlockRows, BlockCols, InnerPanel, Dense>
ykuroda 0:13a5d365ba16 150 : public internal::BlockImpl_dense<XprType, BlockRows, BlockCols, InnerPanel>
ykuroda 0:13a5d365ba16 151 {
ykuroda 0:13a5d365ba16 152 typedef internal::BlockImpl_dense<XprType, BlockRows, BlockCols, InnerPanel> Impl;
ykuroda 0:13a5d365ba16 153 typedef typename XprType::Index Index;
ykuroda 0:13a5d365ba16 154 public:
ykuroda 0:13a5d365ba16 155 typedef Impl Base;
ykuroda 0:13a5d365ba16 156 EIGEN_INHERIT_ASSIGNMENT_OPERATORS(BlockImpl)
ykuroda 0:13a5d365ba16 157 inline BlockImpl(XprType& xpr, Index i) : Impl(xpr,i) {}
ykuroda 0:13a5d365ba16 158 inline BlockImpl(XprType& xpr, Index a_startRow, Index a_startCol) : Impl(xpr, a_startRow, a_startCol) {}
ykuroda 0:13a5d365ba16 159 inline BlockImpl(XprType& xpr, Index a_startRow, Index a_startCol, Index blockRows, Index blockCols)
ykuroda 0:13a5d365ba16 160 : Impl(xpr, a_startRow, a_startCol, blockRows, blockCols) {}
ykuroda 0:13a5d365ba16 161 };
ykuroda 0:13a5d365ba16 162
ykuroda 0:13a5d365ba16 163 namespace internal {
ykuroda 0:13a5d365ba16 164
ykuroda 0:13a5d365ba16 165 /** \internal Internal implementation of dense Blocks in the general case. */
ykuroda 0:13a5d365ba16 166 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel, bool HasDirectAccess> class BlockImpl_dense
ykuroda 0:13a5d365ba16 167 : public internal::dense_xpr_base<Block<XprType, BlockRows, BlockCols, InnerPanel> >::type
ykuroda 0:13a5d365ba16 168 {
ykuroda 0:13a5d365ba16 169 typedef Block<XprType, BlockRows, BlockCols, InnerPanel> BlockType;
ykuroda 0:13a5d365ba16 170 public:
ykuroda 0:13a5d365ba16 171
ykuroda 0:13a5d365ba16 172 typedef typename internal::dense_xpr_base<BlockType>::type Base;
ykuroda 0:13a5d365ba16 173 EIGEN_DENSE_PUBLIC_INTERFACE(BlockType)
ykuroda 0:13a5d365ba16 174 EIGEN_INHERIT_ASSIGNMENT_OPERATORS(BlockImpl_dense)
ykuroda 0:13a5d365ba16 175
ykuroda 0:13a5d365ba16 176 class InnerIterator;
ykuroda 0:13a5d365ba16 177
ykuroda 0:13a5d365ba16 178 /** Column or Row constructor
ykuroda 0:13a5d365ba16 179 */
ykuroda 0:13a5d365ba16 180 inline BlockImpl_dense(XprType& xpr, Index i)
ykuroda 0:13a5d365ba16 181 : m_xpr(xpr),
ykuroda 0:13a5d365ba16 182 // It is a row if and only if BlockRows==1 and BlockCols==XprType::ColsAtCompileTime,
ykuroda 0:13a5d365ba16 183 // and it is a column if and only if BlockRows==XprType::RowsAtCompileTime and BlockCols==1,
ykuroda 0:13a5d365ba16 184 // all other cases are invalid.
ykuroda 0:13a5d365ba16 185 // The case a 1x1 matrix seems ambiguous, but the result is the same anyway.
ykuroda 0:13a5d365ba16 186 m_startRow( (BlockRows==1) && (BlockCols==XprType::ColsAtCompileTime) ? i : 0),
ykuroda 0:13a5d365ba16 187 m_startCol( (BlockRows==XprType::RowsAtCompileTime) && (BlockCols==1) ? i : 0),
ykuroda 0:13a5d365ba16 188 m_blockRows(BlockRows==1 ? 1 : xpr.rows()),
ykuroda 0:13a5d365ba16 189 m_blockCols(BlockCols==1 ? 1 : xpr.cols())
ykuroda 0:13a5d365ba16 190 {}
ykuroda 0:13a5d365ba16 191
ykuroda 0:13a5d365ba16 192 /** Fixed-size constructor
ykuroda 0:13a5d365ba16 193 */
ykuroda 0:13a5d365ba16 194 inline BlockImpl_dense(XprType& xpr, Index a_startRow, Index a_startCol)
ykuroda 0:13a5d365ba16 195 : m_xpr(xpr), m_startRow(a_startRow), m_startCol(a_startCol),
ykuroda 0:13a5d365ba16 196 m_blockRows(BlockRows), m_blockCols(BlockCols)
ykuroda 0:13a5d365ba16 197 {}
ykuroda 0:13a5d365ba16 198
ykuroda 0:13a5d365ba16 199 /** Dynamic-size constructor
ykuroda 0:13a5d365ba16 200 */
ykuroda 0:13a5d365ba16 201 inline BlockImpl_dense(XprType& xpr,
ykuroda 0:13a5d365ba16 202 Index a_startRow, Index a_startCol,
ykuroda 0:13a5d365ba16 203 Index blockRows, Index blockCols)
ykuroda 0:13a5d365ba16 204 : m_xpr(xpr), m_startRow(a_startRow), m_startCol(a_startCol),
ykuroda 0:13a5d365ba16 205 m_blockRows(blockRows), m_blockCols(blockCols)
ykuroda 0:13a5d365ba16 206 {}
ykuroda 0:13a5d365ba16 207
ykuroda 0:13a5d365ba16 208 inline Index rows() const { return m_blockRows.value(); }
ykuroda 0:13a5d365ba16 209 inline Index cols() const { return m_blockCols.value(); }
ykuroda 0:13a5d365ba16 210
ykuroda 0:13a5d365ba16 211 inline Scalar& coeffRef(Index rowId, Index colId)
ykuroda 0:13a5d365ba16 212 {
ykuroda 0:13a5d365ba16 213 EIGEN_STATIC_ASSERT_LVALUE(XprType)
ykuroda 0:13a5d365ba16 214 return m_xpr.const_cast_derived()
ykuroda 0:13a5d365ba16 215 .coeffRef(rowId + m_startRow.value(), colId + m_startCol.value());
ykuroda 0:13a5d365ba16 216 }
ykuroda 0:13a5d365ba16 217
ykuroda 0:13a5d365ba16 218 inline const Scalar& coeffRef(Index rowId, Index colId) const
ykuroda 0:13a5d365ba16 219 {
ykuroda 0:13a5d365ba16 220 return m_xpr.derived()
ykuroda 0:13a5d365ba16 221 .coeffRef(rowId + m_startRow.value(), colId + m_startCol.value());
ykuroda 0:13a5d365ba16 222 }
ykuroda 0:13a5d365ba16 223
ykuroda 0:13a5d365ba16 224 EIGEN_STRONG_INLINE const CoeffReturnType coeff(Index rowId, Index colId) const
ykuroda 0:13a5d365ba16 225 {
ykuroda 0:13a5d365ba16 226 return m_xpr.coeff(rowId + m_startRow.value(), colId + m_startCol.value());
ykuroda 0:13a5d365ba16 227 }
ykuroda 0:13a5d365ba16 228
ykuroda 0:13a5d365ba16 229 inline Scalar& coeffRef(Index index)
ykuroda 0:13a5d365ba16 230 {
ykuroda 0:13a5d365ba16 231 EIGEN_STATIC_ASSERT_LVALUE(XprType)
ykuroda 0:13a5d365ba16 232 return m_xpr.const_cast_derived()
ykuroda 0:13a5d365ba16 233 .coeffRef(m_startRow.value() + (RowsAtCompileTime == 1 ? 0 : index),
ykuroda 0:13a5d365ba16 234 m_startCol.value() + (RowsAtCompileTime == 1 ? index : 0));
ykuroda 0:13a5d365ba16 235 }
ykuroda 0:13a5d365ba16 236
ykuroda 0:13a5d365ba16 237 inline const Scalar& coeffRef(Index index) const
ykuroda 0:13a5d365ba16 238 {
ykuroda 0:13a5d365ba16 239 return m_xpr.const_cast_derived()
ykuroda 0:13a5d365ba16 240 .coeffRef(m_startRow.value() + (RowsAtCompileTime == 1 ? 0 : index),
ykuroda 0:13a5d365ba16 241 m_startCol.value() + (RowsAtCompileTime == 1 ? index : 0));
ykuroda 0:13a5d365ba16 242 }
ykuroda 0:13a5d365ba16 243
ykuroda 0:13a5d365ba16 244 inline const CoeffReturnType coeff(Index index) const
ykuroda 0:13a5d365ba16 245 {
ykuroda 0:13a5d365ba16 246 return m_xpr
ykuroda 0:13a5d365ba16 247 .coeff(m_startRow.value() + (RowsAtCompileTime == 1 ? 0 : index),
ykuroda 0:13a5d365ba16 248 m_startCol.value() + (RowsAtCompileTime == 1 ? index : 0));
ykuroda 0:13a5d365ba16 249 }
ykuroda 0:13a5d365ba16 250
ykuroda 0:13a5d365ba16 251 template<int LoadMode>
ykuroda 0:13a5d365ba16 252 inline PacketScalar packet(Index rowId, Index colId) const
ykuroda 0:13a5d365ba16 253 {
ykuroda 0:13a5d365ba16 254 return m_xpr.template packet<Unaligned>
ykuroda 0:13a5d365ba16 255 (rowId + m_startRow.value(), colId + m_startCol.value());
ykuroda 0:13a5d365ba16 256 }
ykuroda 0:13a5d365ba16 257
ykuroda 0:13a5d365ba16 258 template<int LoadMode>
ykuroda 0:13a5d365ba16 259 inline void writePacket(Index rowId, Index colId, const PacketScalar& val)
ykuroda 0:13a5d365ba16 260 {
ykuroda 0:13a5d365ba16 261 m_xpr.const_cast_derived().template writePacket<Unaligned>
ykuroda 0:13a5d365ba16 262 (rowId + m_startRow.value(), colId + m_startCol.value(), val);
ykuroda 0:13a5d365ba16 263 }
ykuroda 0:13a5d365ba16 264
ykuroda 0:13a5d365ba16 265 template<int LoadMode>
ykuroda 0:13a5d365ba16 266 inline PacketScalar packet(Index index) const
ykuroda 0:13a5d365ba16 267 {
ykuroda 0:13a5d365ba16 268 return m_xpr.template packet<Unaligned>
ykuroda 0:13a5d365ba16 269 (m_startRow.value() + (RowsAtCompileTime == 1 ? 0 : index),
ykuroda 0:13a5d365ba16 270 m_startCol.value() + (RowsAtCompileTime == 1 ? index : 0));
ykuroda 0:13a5d365ba16 271 }
ykuroda 0:13a5d365ba16 272
ykuroda 0:13a5d365ba16 273 template<int LoadMode>
ykuroda 0:13a5d365ba16 274 inline void writePacket(Index index, const PacketScalar& val)
ykuroda 0:13a5d365ba16 275 {
ykuroda 0:13a5d365ba16 276 m_xpr.const_cast_derived().template writePacket<Unaligned>
ykuroda 0:13a5d365ba16 277 (m_startRow.value() + (RowsAtCompileTime == 1 ? 0 : index),
ykuroda 0:13a5d365ba16 278 m_startCol.value() + (RowsAtCompileTime == 1 ? index : 0), val);
ykuroda 0:13a5d365ba16 279 }
ykuroda 0:13a5d365ba16 280
ykuroda 0:13a5d365ba16 281 #ifdef EIGEN_PARSED_BY_DOXYGEN
ykuroda 0:13a5d365ba16 282 /** \sa MapBase::data() */
ykuroda 0:13a5d365ba16 283 inline const Scalar* data() const;
ykuroda 0:13a5d365ba16 284 inline Index innerStride() const;
ykuroda 0:13a5d365ba16 285 inline Index outerStride() const;
ykuroda 0:13a5d365ba16 286 #endif
ykuroda 0:13a5d365ba16 287
ykuroda 0:13a5d365ba16 288 const typename internal::remove_all<typename XprType::Nested>::type& nestedExpression() const
ykuroda 0:13a5d365ba16 289 {
ykuroda 0:13a5d365ba16 290 return m_xpr;
ykuroda 0:13a5d365ba16 291 }
ykuroda 0:13a5d365ba16 292
ykuroda 0:13a5d365ba16 293 Index startRow() const
ykuroda 0:13a5d365ba16 294 {
ykuroda 0:13a5d365ba16 295 return m_startRow.value();
ykuroda 0:13a5d365ba16 296 }
ykuroda 0:13a5d365ba16 297
ykuroda 0:13a5d365ba16 298 Index startCol() const
ykuroda 0:13a5d365ba16 299 {
ykuroda 0:13a5d365ba16 300 return m_startCol.value();
ykuroda 0:13a5d365ba16 301 }
ykuroda 0:13a5d365ba16 302
ykuroda 0:13a5d365ba16 303 protected:
ykuroda 0:13a5d365ba16 304
ykuroda 0:13a5d365ba16 305 const typename XprType::Nested m_xpr;
ykuroda 0:13a5d365ba16 306 const internal::variable_if_dynamic<Index, XprType::RowsAtCompileTime == 1 ? 0 : Dynamic> m_startRow;
ykuroda 0:13a5d365ba16 307 const internal::variable_if_dynamic<Index, XprType::ColsAtCompileTime == 1 ? 0 : Dynamic> m_startCol;
ykuroda 0:13a5d365ba16 308 const internal::variable_if_dynamic<Index, RowsAtCompileTime> m_blockRows;
ykuroda 0:13a5d365ba16 309 const internal::variable_if_dynamic<Index, ColsAtCompileTime> m_blockCols;
ykuroda 0:13a5d365ba16 310 };
ykuroda 0:13a5d365ba16 311
ykuroda 0:13a5d365ba16 312 /** \internal Internal implementation of dense Blocks in the direct access case.*/
ykuroda 0:13a5d365ba16 313 template<typename XprType, int BlockRows, int BlockCols, bool InnerPanel>
ykuroda 0:13a5d365ba16 314 class BlockImpl_dense<XprType,BlockRows,BlockCols, InnerPanel,true>
ykuroda 0:13a5d365ba16 315 : public MapBase<Block<XprType, BlockRows, BlockCols, InnerPanel> >
ykuroda 0:13a5d365ba16 316 {
ykuroda 0:13a5d365ba16 317 typedef Block<XprType, BlockRows, BlockCols, InnerPanel> BlockType;
ykuroda 0:13a5d365ba16 318 public:
ykuroda 0:13a5d365ba16 319
ykuroda 0:13a5d365ba16 320 typedef MapBase<BlockType> Base;
ykuroda 0:13a5d365ba16 321 EIGEN_DENSE_PUBLIC_INTERFACE(BlockType)
ykuroda 0:13a5d365ba16 322 EIGEN_INHERIT_ASSIGNMENT_OPERATORS(BlockImpl_dense)
ykuroda 0:13a5d365ba16 323
ykuroda 0:13a5d365ba16 324 /** Column or Row constructor
ykuroda 0:13a5d365ba16 325 */
ykuroda 0:13a5d365ba16 326 inline BlockImpl_dense(XprType& xpr, Index i)
ykuroda 0:13a5d365ba16 327 : Base(internal::const_cast_ptr(&xpr.coeffRef(
ykuroda 0:13a5d365ba16 328 (BlockRows==1) && (BlockCols==XprType::ColsAtCompileTime) ? i : 0,
ykuroda 0:13a5d365ba16 329 (BlockRows==XprType::RowsAtCompileTime) && (BlockCols==1) ? i : 0)),
ykuroda 0:13a5d365ba16 330 BlockRows==1 ? 1 : xpr.rows(),
ykuroda 0:13a5d365ba16 331 BlockCols==1 ? 1 : xpr.cols()),
ykuroda 0:13a5d365ba16 332 m_xpr(xpr)
ykuroda 0:13a5d365ba16 333 {
ykuroda 0:13a5d365ba16 334 init();
ykuroda 0:13a5d365ba16 335 }
ykuroda 0:13a5d365ba16 336
ykuroda 0:13a5d365ba16 337 /** Fixed-size constructor
ykuroda 0:13a5d365ba16 338 */
ykuroda 0:13a5d365ba16 339 inline BlockImpl_dense(XprType& xpr, Index startRow, Index startCol)
ykuroda 0:13a5d365ba16 340 : Base(internal::const_cast_ptr(&xpr.coeffRef(startRow,startCol))), m_xpr(xpr)
ykuroda 0:13a5d365ba16 341 {
ykuroda 0:13a5d365ba16 342 init();
ykuroda 0:13a5d365ba16 343 }
ykuroda 0:13a5d365ba16 344
ykuroda 0:13a5d365ba16 345 /** Dynamic-size constructor
ykuroda 0:13a5d365ba16 346 */
ykuroda 0:13a5d365ba16 347 inline BlockImpl_dense(XprType& xpr,
ykuroda 0:13a5d365ba16 348 Index startRow, Index startCol,
ykuroda 0:13a5d365ba16 349 Index blockRows, Index blockCols)
ykuroda 0:13a5d365ba16 350 : Base(internal::const_cast_ptr(&xpr.coeffRef(startRow,startCol)), blockRows, blockCols),
ykuroda 0:13a5d365ba16 351 m_xpr(xpr)
ykuroda 0:13a5d365ba16 352 {
ykuroda 0:13a5d365ba16 353 init();
ykuroda 0:13a5d365ba16 354 }
ykuroda 0:13a5d365ba16 355
ykuroda 0:13a5d365ba16 356 const typename internal::remove_all<typename XprType::Nested>::type& nestedExpression() const
ykuroda 0:13a5d365ba16 357 {
ykuroda 0:13a5d365ba16 358 return m_xpr;
ykuroda 0:13a5d365ba16 359 }
ykuroda 0:13a5d365ba16 360
ykuroda 0:13a5d365ba16 361 /** \sa MapBase::innerStride() */
ykuroda 0:13a5d365ba16 362 inline Index innerStride() const
ykuroda 0:13a5d365ba16 363 {
ykuroda 0:13a5d365ba16 364 return internal::traits<BlockType>::HasSameStorageOrderAsXprType
ykuroda 0:13a5d365ba16 365 ? m_xpr.innerStride()
ykuroda 0:13a5d365ba16 366 : m_xpr.outerStride();
ykuroda 0:13a5d365ba16 367 }
ykuroda 0:13a5d365ba16 368
ykuroda 0:13a5d365ba16 369 /** \sa MapBase::outerStride() */
ykuroda 0:13a5d365ba16 370 inline Index outerStride() const
ykuroda 0:13a5d365ba16 371 {
ykuroda 0:13a5d365ba16 372 return m_outerStride;
ykuroda 0:13a5d365ba16 373 }
ykuroda 0:13a5d365ba16 374
ykuroda 0:13a5d365ba16 375 #ifndef __SUNPRO_CC
ykuroda 0:13a5d365ba16 376 // FIXME sunstudio is not friendly with the above friend...
ykuroda 0:13a5d365ba16 377 // META-FIXME there is no 'friend' keyword around here. Is this obsolete?
ykuroda 0:13a5d365ba16 378 protected:
ykuroda 0:13a5d365ba16 379 #endif
ykuroda 0:13a5d365ba16 380
ykuroda 0:13a5d365ba16 381 #ifndef EIGEN_PARSED_BY_DOXYGEN
ykuroda 0:13a5d365ba16 382 /** \internal used by allowAligned() */
ykuroda 0:13a5d365ba16 383 inline BlockImpl_dense(XprType& xpr, const Scalar* data, Index blockRows, Index blockCols)
ykuroda 0:13a5d365ba16 384 : Base(data, blockRows, blockCols), m_xpr(xpr)
ykuroda 0:13a5d365ba16 385 {
ykuroda 0:13a5d365ba16 386 init();
ykuroda 0:13a5d365ba16 387 }
ykuroda 0:13a5d365ba16 388 #endif
ykuroda 0:13a5d365ba16 389
ykuroda 0:13a5d365ba16 390 protected:
ykuroda 0:13a5d365ba16 391 void init()
ykuroda 0:13a5d365ba16 392 {
ykuroda 0:13a5d365ba16 393 m_outerStride = internal::traits<BlockType>::HasSameStorageOrderAsXprType
ykuroda 0:13a5d365ba16 394 ? m_xpr.outerStride()
ykuroda 0:13a5d365ba16 395 : m_xpr.innerStride();
ykuroda 0:13a5d365ba16 396 }
ykuroda 0:13a5d365ba16 397
ykuroda 0:13a5d365ba16 398 typename XprType::Nested m_xpr;
ykuroda 0:13a5d365ba16 399 Index m_outerStride;
ykuroda 0:13a5d365ba16 400 };
ykuroda 0:13a5d365ba16 401
ykuroda 0:13a5d365ba16 402 } // end namespace internal
ykuroda 0:13a5d365ba16 403
ykuroda 0:13a5d365ba16 404 } // end namespace Eigen
ykuroda 0:13a5d365ba16 405
ykuroda 0:13a5d365ba16 406 #endif // EIGEN_BLOCK_H