1*bf2c3715SXin Li Matrix2d a, b, c; a << 1,2,3,4; b << 5,6,7,8; 2*bf2c3715SXin Li c.noalias() = a * b; // this computes the product directly to c 3*bf2c3715SXin Li cout << c << endl; 4