1 MatrixXf matA(2,2);
2 matA << 2, 0,  0, 2;
3 matA.noalias() = matA * matA;
4 cout << matA;
5