1namespace Eigen { 2 3/** \page TopicStorageOrders Storage orders 4 5There are two different storage orders for matrices and two-dimensional arrays: column-major and row-major. 6This page explains these storage orders and how to specify which one should be used. 7 8<b>Table of contents</b> 9 - \ref TopicStorageOrdersIntro 10 - \ref TopicStorageOrdersInEigen 11 - \ref TopicStorageOrdersWhich 12 13 14\section TopicStorageOrdersIntro Column-major and row-major storage 15 16The entries of a matrix form a two-dimensional grid. However, when the matrix is stored in memory, the entries 17have to somehow be laid out linearly. There are two main ways to do this, by row and by column. 18 19We say that a matrix is stored in \b row-major order if it is stored row by row. The entire first row is 20stored first, followed by the entire second row, and so on. Consider for example the matrix 21 22\f[ 23A = \begin{bmatrix} 248 & 2 & 2 & 9 \\ 259 & 1 & 4 & 4 \\ 263 & 5 & 4 & 5 27\end{bmatrix}. 28\f] 29 30If this matrix is stored in row-major order, then the entries are laid out in memory as follows: 31 32\code 8 2 2 9 9 1 4 4 3 5 4 5 \endcode 33 34On the other hand, a matrix is stored in \b column-major order if it is stored column by column, starting with 35the entire first column, followed by the entire second column, and so on. If the above matrix is stored in 36column-major order, it is laid out as follows: 37 38\code 8 9 3 2 1 5 2 4 4 9 4 5 \endcode 39 40This example is illustrated by the following Eigen code. It uses the PlainObjectBase::data() function, which 41returns a pointer to the memory location of the first entry of the matrix. 42 43<table class="example"> 44<tr><th>Example</th><th>Output</th></tr> 45<tr><td> 46\include TopicStorageOrders_example.cpp 47</td> 48<td> 49\verbinclude TopicStorageOrders_example.out 50</td></tr></table> 51 52 53\section TopicStorageOrdersInEigen Storage orders in Eigen 54 55The storage order of a matrix or a two-dimensional array can be set by specifying the \c Options template 56parameter for Matrix or Array. As \ref TutorialMatrixClass explains, the %Matrix class template has six 57template parameters, of which three are compulsory (\c Scalar, \c RowsAtCompileTime and \c ColsAtCompileTime) 58and three are optional (\c Options, \c MaxRowsAtCompileTime and \c MaxColsAtCompileTime). If the \c Options 59parameter is set to \c RowMajor, then the matrix or array is stored in row-major order; if it is set to 60\c ColMajor, then it is stored in column-major order. This mechanism is used in the above Eigen program to 61specify the storage order. 62 63If the storage order is not specified, then Eigen defaults to storing the entry in column-major. This is also 64the case if one of the convenience typedefs (\c Matrix3f, \c ArrayXXd, etc.) is used. 65 66Matrices and arrays using one storage order can be assigned to matrices and arrays using the other storage 67order, as happens in the above program when \c Arowmajor is initialized using \c Acolmajor. Eigen will reorder 68the entries automatically. More generally, row-major and column-major matrices can be mixed in an expression 69as we want. 70 71 72\section TopicStorageOrdersWhich Which storage order to choose? 73 74So, which storage order should you use in your program? There is no simple answer to this question; it depends 75on your application. Here are some points to keep in mind: 76 77 - Your users may expect you to use a specific storage order. Alternatively, you may use other libraries than 78 Eigen, and these other libraries may expect a certain storage order. In these cases it may be easiest and 79 fastest to use this storage order in your whole program. 80 - Algorithms that traverse a matrix row by row will go faster when the matrix is stored in row-major order 81 because of better data locality. Similarly, column-by-column traversal is faster for column-major 82 matrices. It may be worthwhile to experiment a bit to find out what is faster for your particular 83 application. 84 - The default in Eigen is column-major. Naturally, most of the development and testing of the Eigen library 85 is thus done with column-major matrices. This means that, even though we aim to support column-major and 86 row-major storage orders transparently, the Eigen library may well work best with column-major matrices. 87 88*/ 89} 90