- make MatrixBase and all expressions aware of their preferred traversal order.
  Honor this preference in operator=.
- add several methods to the API
- rework API for diagonal matrices
- add benchmarking code
diff --git a/test/miscmatrices.cpp b/test/miscmatrices.cpp
index 052d5ec..7126266 100644
--- a/test/miscmatrices.cpp
+++ b/test/miscmatrices.cpp
@@ -46,7 +46,7 @@
   VectorType v1 = VectorType::random(rows);
   v1[0];
   Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime>
-  square = MatrixType::diagonal(v1);
+  square = v1.asDiagonal();
   if(r==r2) VERIFY_IS_APPROX(square(r,r2), v1[r]);
   else VERIFY_IS_MUCH_SMALLER_THAN(square(r,r2), static_cast<Scalar>(1));
   square = MatrixType::zero(rows, rows);