- 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/doc/snippets/MatrixBase_asDiagonal.cpp b/doc/snippets/MatrixBase_asDiagonal.cpp
new file mode 100644
index 0000000..4cbff82
--- /dev/null
+++ b/doc/snippets/MatrixBase_asDiagonal.cpp
@@ -0,0 +1 @@
+cout << Vector3i(2,5,6).asDiagonal() << endl;