RMUL2025/lib/eigen/doc/snippets/SelfAdjointView_operatorNorm.cpp

MatrixXd ones = MatrixXd::Ones(3, 3);
cout << "The operator norm of the 3x3 matrix of ones is " << ones.selfadjointView<Lower>().operatorNorm() << endl;