Eigen  3.4.90 (git rev a4098ac676528a83cfb73d4d26ce1b42ec05f47c)
Interfacing with raw buffers: the Map class

This page explains how to work with "raw" C/C++ arrays. This can be useful in a variety of contexts, particularly when "importing" vectors and matrices from other libraries into Eigen.


Occasionally you may have a pre-defined array of numbers that you want to use within Eigen as a vector or matrix. While one option is to make a copy of the data, most commonly you probably want to re-use this memory as an Eigen type. Fortunately, this is very easy with the Map class.

Map types and declaring Map variables

A Map object has a type defined by its Eigen equivalent:

A matrix or vector expression mapping an existing array of data.
Definition: Map.h:98

Note that, in this default case, a Map requires just a single template parameter.

To construct a Map variable, you need two other pieces of information: a pointer to the region of memory defining the array of coefficients, and the desired shape of the matrix or vector. For example, to define a matrix of float with sizes determined at compile time, you might do the following:

Map<MatrixXf> mf(pf,rows,columns);

where pf is a float * pointing to the array of memory. A fixed-size read-only vector of integers might be declared as

where pi is an int *. In this case the size does not have to be passed to the constructor, because it is already specified by the Matrix/Array type.

Note that Map does not have a default constructor; you must pass a pointer to initialize the object. However, you can work around this requirement (see Changing the mapped array).

Map is flexible enough to accommodate a variety of different data representations. There are two other (optional) template parameters:

Map<typename MatrixType,
int MapOptions,
typename StrideType>

Using Map variables

You can use a Map object just like any other Eigen type:

typedef Matrix<float,1,Dynamic> MatrixType;
typedef Map<MatrixType> MapType;
typedef Map<const MatrixType> MapTypeConst; // a read-only map
const int n_dims = 5;
MatrixType m1(n_dims), m2(n_dims);
float *p = &m2(0); // get the address storing the data for m2
MapType m2map(p,m2.size()); // m2map shares data with m2
MapTypeConst m2mapconst(p,m2.size()); // a read-only accessor for m2
cout << "m1: " << m1 << endl;
cout << "m2: " << m2 << endl;
cout << "Squared euclidean distance: " << (m1-m2).squaredNorm() << endl;
cout << "Squared euclidean distance, using map: " <<
(m1-m2map).squaredNorm() << endl;
m2map(3) = 7; // this will change m2, since they share the same array
cout << "Updated m2: " << m2 << endl;
cout << "m2 coefficient 2, constant accessor: " << m2mapconst(2) << endl;
/* m2mapconst(2) = 5; */ // this yields a compile-time error
The matrix class, also used for vectors and row-vectors.
Definition: Matrix.h:182
m1:      -1  -0.737   0.511 -0.0827  0.0655
m2: -0.562 -0.906  0.358  0.359  0.869
Squared euclidean distance: 1.08
Squared euclidean distance, using map: 1.08
Updated m2: -0.562 -0.906  0.358      7  0.869
m2 coefficient 2, constant accessor: 0.358

All Eigen functions are written to accept Map objects just like other Eigen types. However, when writing your own functions taking Eigen types, this does not happen automatically: a Map type is not identical to its Dense equivalent. See Writing Functions Taking Eigen Types as Parameters for details.

Changing the mapped array

It is possible to change the array of a Map object after declaration, using the C++ "placement new" syntax:

int data[] = {1,2,3,4,5,6,7,8,9};
Map<RowVectorXi> v(data,4);
cout << "The mapped vector v is: " << v << "\n";
new (&v) Map<RowVectorXi>(data+4,5);
cout << "Now v is: " << v << "\n";
The mapped vector v is: 1 2 3 4
Now v is: 5 6 7 8 9

Despite appearances, this does not invoke the memory allocator, because the syntax specifies the location for storing the result.

This syntax makes it possible to declare a Map object without first knowing the mapped array's location in memory:

Map<Matrix3f> A(NULL); // don't try to use this matrix yet!
VectorXf b(n_matrices);
for (int i = 0; i < n_matrices; i++)
new (&A) Map<Matrix3f>(get_matrix_pointer(i));
b(i) = A.trace();