1 typedef Matrix<double,3,Dynamic> Matrix3Xd; 2 Matrix3Xd M = Matrix3Xd::Random(3,5); 3 Projective3d P(Matrix4d::Random()); 4 cout << "The matrix M is:" << endl << M << endl << endl; 5 cout << "M.colwise().homogeneous():" << endl << M.colwise().homogeneous() << endl << endl; 6 cout << "P * M.colwise().homogeneous():" << endl << P * M.colwise().homogeneous() << endl << endl; 7 cout << "P * M.colwise().homogeneous().hnormalized(): " << endl << (P * M.colwise().homogeneous()).colwise().hnormalized() << endl << endl;