matrix storage order can now also be row-dominant (choosable for each matrix separately)

map() moves from MatrixBase to Matrix
much more documentation/examples/snippets

doc/Doxyfile.in

@@ -237,13 +237,13 @@ PERLMOD_MAKEVAR_PREFIX =
 # Configuration options related to the preprocessor   
 #---------------------------------------------------------------------------
 ENABLE_PREPROCESSING   = YES
-MACRO_EXPANSION        = NO
-EXPAND_ONLY_PREDEF     = NO
+MACRO_EXPANSION        = YES
+EXPAND_ONLY_PREDEF     = YES
 SEARCH_INCLUDES        = YES
 INCLUDE_PATH           = 
 INCLUDE_FILE_PATTERNS  = 
 PREDEFINED             = 
-EXPAND_AS_DEFINED      = 
+EXPAND_AS_DEFINED      = EIGEN_MAKE_SCALAR_OPS
 SKIP_FUNCTION_MACROS   = YES
 #---------------------------------------------------------------------------
 # Configuration::additions related to external references   

doc/examples/class_Eval.cpp

@@ -0,0 +1,28 @@
+#include <Eigen/Core.h>
+USING_PART_OF_NAMESPACE_EIGEN
+using namespace std;
+
+template<typename Scalar, typename Derived>
+const Eigen::Eval<Eigen::Transpose<Derived> >
+evaluatedTranspose(const MatrixBase<Scalar, Derived>& m)
+{
+  return m.transpose().eval();
+}
+
+int main(int, char**)
+{
+  Matrix2f M = Matrix2f::random();
+  Matrix2f m;
+  m = M;
+  cout << "Here is the matrix m:" << endl << m << endl;
+  cout << "Now we want to replace m by its own transpose." << endl;
+  cout << "If we do m = m.transpose(), then m becomes:" << endl;
+  m = m.transpose();
+  cout << m << endl << "which is wrong!" << endl;
+  cout << "Now let us instead do m = evaluatedTranspose(m). Then m becomes" << endl;
+  m = M;
+  m = evaluatedTranspose(m);
+  cout << m << endl << "which is right." << endl;
+
+  return 0;
+}

doc/snippets/MatrixBase_adjoint.cpp

@@ -0,0 +1,3 @@
+Matrix2cf m = Matrix2cf::random();
+cout << "Here is the 2x2 complex matrix m:" << endl << m << endl;
+cout << "Here is the adjoint of m:" << endl << m.adjoint() << endl;

doc/snippets/MatrixBase_block.cpp

@@ -1,3 +1,5 @@
 Matrix4d m = Matrix4d::diagonal(Vector4d(1,2,3,4));
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Here is m.block<2, 2>(2, 2):" << endl << m.block<2, 2>(2, 2) << endl;
 m.block<2, 2>(2, 0) = m.block<2, 2>(2, 2);
-cout << m << endl;
+cout << "Now the matrix m is:" << endl << m << endl;

doc/snippets/MatrixBase_diagonal.cpp

@@ -0,0 +1,4 @@
+Matrix3i m = Matrix3i::random();
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Here are the coefficients on the main diagonal of m:" << endl
+     << m.diagonal() << endl;

doc/snippets/MatrixBase_dynBlock.cpp

@@ -1,3 +1,5 @@
 Matrix3d m = Matrix3d::diagonal(Vector3d(1,2,3));
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Here is m.dynBlock(1, 1, 2, 1):" << endl << m.dynBlock(1, 1, 2, 1) << endl;
 m.dynBlock(1, 0, 2, 1) = m.dynBlock(1, 1, 2, 1);
-cout << m << endl;
+cout << "Now the matrix m is:" << endl << m << endl;

doc/snippets/MatrixBase_eval.cpp

@@ -0,0 +1,12 @@
+Matrix2f M = Matrix2f::random();
+Matrix2f m;
+m = M;
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Now we want to replace m by its own transpose." << endl;
+cout << "If we do m = m.transpose(), then m becomes:" << endl;
+m = m.transpose();
+cout << m << endl << "which is wrong!" << endl;
+cout << "Now let us instead do m = m.transpose().eval(). Then m becomes" << endl;
+m = M;
+m = m.transpose().eval();
+cout << m << endl << "which is right." << endl;

doc/snippets/MatrixBase_minor.cpp

@@ -0,0 +1,3 @@
+Matrix3i m = Matrix3i::random();
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Here is m.minor(1,1):" << endl << m.minor(1,1) << endl;

doc/snippets/MatrixBase_transpose.cpp

@@ -0,0 +1,8 @@
+Matrix2f m = Matrix2f::random();
+cout << "Here is the matrix m:" << endl << m << endl;
+cout << "Here is the transpose of m:" << endl << m.transpose() << endl;
+cout << "Here is the coefficient (1,0) in the transpose of m:" << endl
+     << m.transpose()(1,0) << endl;
+cout << "Let us overwrite this coefficient with the value 0." << endl;
+m.transpose()(1,0) = 0;
+cout << "Now the matrix m is:" << endl << m << endl;